Temporal analysis of blazar flux is a powerful tool to draw inferences about the emission processes and physics of these sources. In the most general case, the available light curves are irregularly sampled and influenced by gaps, and in addition are also affected by correlated noise, making their analysis complicated. Gaussian processes may offer a viable tool to assess the statistical significance of proposed periods in light curves characterized by any sampling and noise pattern. We infer the significance of the periods proposed in the literature for two well known blazars with multiple claims of possible year-long periodicity: PG1553+113 and PKS2155-304, in the high-energy and optical bands. Adding a periodic component to the modeling gives a better statistical description of the analyzed light curves. The improvement is rather solid for PG1553+113, both at high energies and in the optical, while for PKS2155-304 at high energies the improvement is not yet strong enough to allow cogent claims, and no evidence for periodicity emerged by the analysis in the optical. Modeling a light curve by means of Gaussian processes, in spite of being relatively computationally demanding, allows us to derive a wealth of information about the data under study and suggests an original analysis framework for light curves of astrophysical interest.

We investigate the environments of galaxies around BL Lac objects and FR I radiogalaxies, the alleged parent populations of misaligned sources. We compare the environment of a sample of 50 BL Lac objects at 0.1<z<0.33 with that of a sample of 90 FR I galaxies at 0.1<z< 0.15. The galaxy environment is estimated using the SDSS images in the i-band. We find that the galaxy excess density within 0.5 Mpc around of FR I radiogalaxies is a factor ~ 2 larger than that around BL Lacs. This implies a reconsideration of the parent population of BL Lac objects.

The availability of about a decade of uninterrupted sky monitoring by the Fermi satellite has made possible to study long-term quasi-periodicities for high-energy sources. It is therefore not a surprise that for several blazars in the recent literature claims for such periodicities, with various level of confidence, have been reported. The confirmation of these findings could be of tremendous importance for the physical description of this category of sources and have consequences for the gravitational wave background interpretation. In this work we carry out a temporal analysis of the Fermi light curves for several of the sources mentioned in recent literature by means of a homogeneous procedure and find that, globally, no strong cases for blazar year-long quasi-periodicities can be confirmed. The computed power spectral densities are all essentially consistent with being generated by red-noise only. We further discuss the meaning and the limitations of the present analysis.

PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5{\sigma}) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broad-band emission is compatible with the EC scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR.

Blazars are a sub-population of quasars whose jets are nearly aligned with the line-of-sight, which tend to exhibit multi-wavelength variability on a variety of timescales. Quasi-periodic variability on year-like timescales has been detected in a number of bright sources, and has been connected to the orbital motion of a putative massive black hole binary. If this were indeed the case, those blazar binaries would contribute to the nanohertz gravitational-wave stochastic background. We test the binary hypothesis for the blazar population observed by the \textit{Fermi} Gamma-Ray Space Telescope, which consists of BL Lacertae objects and flat-spectrum radio quasars. Using mock populations informed by the luminosity functions for BL Lacertae objects and flat-spectrum radio quasars with redshifts $z \le 2$, we calculate the expected gravitational wave background and compare it to recent pulsar timing array upper limits. The two are consistent only if a fraction $\lesssim 10^{-3}$ of blazars hosts a binary with orbital periods $<5$ years. We therefore conclude that binarity cannot significantly explain year-like quasi-periodicity in blazars.

We review the reports of possible year-long quasi-periodicities of BL Lac objects in the $\gamma$-ray and optical bands, and present a homogeneous time analysis of the light curves of PKS2155$-$304, PG1553+113, and BL Lac. Based on results from a survey covering the entire Fermi $\gamma$-ray sky we have estimated the fraction of possible quasi-periodic BL Lac objects. We compared the cyclical behaviour in BL Lac objects with that derived from the search of possible optical periodicities in quasars, and find that at z$\lesssim$1 the cosmic density of quasi-periodic BL Lac objects is larger than that of quasi-periodic quasars. If the BL Lac quasi-periodicities were due to a supermassive binary black hole (SBBH) scenario, there could be a tension with the upper limits on the gravitational wave background measured by the pulsar timing array. The argument clearly indicates the difficulties of generally associating quasi-periodicities of BL Lac objects with SBBHs.

We report on a study of the galaxy environments of low redshift physical quasars pairs. We selected 20 pairs having projected separation $<$ 0.5 Mpc and difference of systemic velocity $<$ 800 km s$^{-1}$. Using Sloan Digital Sky Survey images we evaluated the galaxy overdensity around these quasars in pairs and then compare it with that of a sample of isolated quasars with same redshift and luminosity. It is found that on average there is a systematic larger overdensity of galaxies around quasars in pairs with respect to that of isolated quasars. This may represent a significant link between nuclear activity and galaxy environment. However, at odds with that, the closest quasar pairs seem to inhabit poorer environments. Implications of present results and perspectives for future work are briefly discussed.

We examine the 2008-2016 gamma-ray and optical light curves of a number of bright Fermi blazars. In a fraction of them, the periodograms show possible evidence of quasi-periodicities related in the two bands. This coincidence strengthens their physical meaning. Comparing with results from the periodicity search of quasars, the presence of quasi-periodicities in blazars suggests that the basic condition for its observability is related to the relativistic jet in the observer direction, but the overall picture remains uncertain.

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 have considered the \textit{Fermi} $\gamma$-ray light curve of the blazar PKS 0537-441 during a high state extending from 2008/08/10 to 2011/08/27. The periodogram exhibits a peak at T $\sim$ 280 d, with a significance of $\sim$ 99.7 \%. A peak of similar relevance at $\frac{1}{2}$ T is found in the optical light curves. Considering the entire duration of the \textit{ Fermi} light curve 2008-2015, no significant peak is revealed, while the optical one remains meaningful. Comparing with recent observations of PKS 2155-304 and PG 1553+113 it seems that month-year oscillations can characterize some blazars. Month-scale oscillations can also show up only during phases of enhanced or bursting emission like in the case of PKS 0537-441.

We searched for quasi-periodicities on year-like timescales in the light curves of 6 blazars in the optical - near infrared bands and we made a comparison with the high energy emission. We obtained optical/NIR light curves from REM photometry plus archival SMARTS data and we accessed the Fermi light curves for the $\gamma$-ray data. The periodograms often show strong peaks in the optical and gamma-ray bands, which in some cases may be inter-related. The significance of the revealed peaks is then discussed, taking into account that the noise is frequency dependent. Quasi-periodicities on a year-like timescale appear to occur often in blazars. No straightforward model describing these possible periodicities is yet available, but some plausible interpretations for the physical mechanisms causing periodic variabilities of these sources are examined.

Blazars are astrophysical sources whose emission is dominated by non-thermal processes, typically interpreted as synchrotron and inverse Compton emission. Although the general picture is rather robust and consistent with observations, many aspects are still unexplored. Polarimetric monitoring can offer a wealth of information about the physical processes in blazars. Models with largely different physical ingredients can often provide almost indistinguishable predictions for the total flux, but usually are characterized by markedly different polarization properties. We explore, with a pilot study, the possibility to derive structural information about the emitting regions of blazars by means of a joint analysis of rapid variability of the total and polarized flux at optical wavelengths. Short timescale (from tens of seconds to a couple of minutes) optical linear polarimetry and photometry for two blazars, BL Lacertae and PKS 1424+240, was carried out with the PAOLO polarimeter at the 3.6m Telescopio Nazionale Galileo. Several hours of almost continuous observations were obtained for both sources. Our intense monitoring allowed us to draw strongly different scenarios for BL Lacertae and PKS 1424+240, with the former characterized by intense variability on time-scales from hours to a few minutes and the latter practically constant in total flux. Essentially the same behavior is observed for the polarized flux and the position angle. The variability time-scales turned out to be as short as a few minutes, although involving only a few percent variation of the flux. The polarization variability time-scale is generally consistent with the total flux variability. Total and polarized flux appear to be essentially uncorrelated. However, even during our relatively short monitoring, different regimes can be singled out. (abridged)

We have searched for periodicities in our VRIJHK photometry of PKS 2155-304, which covers the years 2005-2012. A peak of the Fourier spectrum with high significance is found at T$\sim$315 days, confirming the recent findings by \cite{Zhang2014}. The examination of the gamma-ray light curves from the\textit{ Fermi} archives yields a significant signal at $\sim$ 2T, which, while nominally significant, involves data spanning only $\sim$ 6T. Assuming a black hole mass of $10^{9}\ M_{\odot}$ the Keplerian distance corresponding to the quasi-period T is $\sim10^{16}$ cm, about 50 Schwarzschild radii.

We investigate the properties of the galaxy environment of a sample of 14 low redshift (z $<$ 0.85) quasar physical pairs extracted from SDSS DR10 archives. The pairs have a systemic radial velocity difference $\Delta V_\parallel \leqslant$ 600 $km \ s^{-1}$ (based on [OIII]5007 \AA \ line) and projected distance $ R_\bot \leqslant$ 600 kpc. The physical association of the pairs is statistically confirmed at a level of $\sim$ 90 %. For most of the images of these quasars we are able to resolve their host galaxies that turn out to be on average similar to those of quasars not in pairs. We also found that quasars in a pair are on average in region of modest galaxy overdensity extending up 0.5 Mpc from the QSO. This galaxy overdensity is indistinguishable from that of a homogeneous sample of isolated quasars at the same redshift and with similar host galaxy luminosity. These results, albeit derived from a small (but homogeneous) sample of objects, suggest that the rare activation of two quasars with small physical separation does not require any extraordinary environment.

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

We present Very Large Telescope optical spectroscopy of nine BL Lac objects of unknown redshift belonging to the list of optically selected radio loud BL Lacs candidates. We explore their spectroscopic properties and the possible link with gamma ray emission. From the new observations we determine the redshift of four objects from faint emission lines or from absorption features of the host galaxy. In three cases we find narrow intervening absorptions from which a lower limit to the redshift is inferred. For the remaining two featureless sources, lower limits to the redshift are deduced from the very absence of spectral lines. A search for gamma counterpart emission shows that six out of nine are Fermi gamma-ray emitters with two new detections. Our analysis suggests that most of the BL Lac still lacking of redshift information are most probably located at high redshift.

Context. We present the light curves of six BL Lac objects, PKS 0537-441, PKS 0735+17, OJ 287, PKS 2005-489, PKS 2155-304, W Comae, and of the flat spectrum radio quasar PKS 1510-089, as a part of a photometric monitoring program in the near-infrared/optical bands started in 2004. All sources are Fermi blazars. Aims. Our purpose is to investigate flux and spectral variability on short and long time scales. Systematic monitoring, independent of the activity of the source, guarantees large sample size statistics, and allows an unbiased view of different activity states on weekly or daily time scales for the whole timeframe and on nightly timescales for some epochs. Methods. Data were obtained with the REM telescope located at the ESO premises of La Silla (Chile). Light curves were gathered in the optical/near-infrared VRIJHK bands from April 2005 to June 2012. Results. Variability >= 3 mag is observed in PKS 0537-441, PKS 1510-089 and PKS 2155-304, the largest ranges spanned in the near-infrared. The color intensity plots show rather different morphologies. The spectral energy distributions in general are well fitted by a power law, with some deviations that are more apparent in low states. Some variability episodes during a night interval are well documented for PKS 0537-441 and PKS 2155-304. For the latter source the variability time scale implies a large relativistic beaming factor.