In this work, we study the probability density functions that best describe the $\gamma$-ray flux distribution of decade-long light curves of a sample of blazars. For averaged behavior over this period, there was maximum likelihood estimated log-stable distribution; for most sources leading to standard log-normal distribution ($\alpha=2$); however, other sources clearly displayed heavy tail distributions ($\alpha<2$), suggesting underlying multiplicative process of infinite variance. For sequences normalized using the found log-stable distributions, there was performed proposed novel non-stationarity and autocorrelation analysis. The former allowed to quantitatively evaluate non-stationarity of each source as maximizing log-likelihood rate of forgetting in modeled evolution of PDFs, also for evaluation of local variability allowing e.g. for anomaly detection suggesting changes of behavior. Discussed autocorrelation analysis looked at lag $l$ dependence of statistical behavior of all $\{(y_t,y_{t+l})\}$ points, described by various mixed moments - allowing to quantitatively point multiple characteristic time scales of the objects, for example, suggesting hidden periodic processes, with statistical interpretations of their contributions.

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.

The search for periodic signals from blazars has become an actively pursued field of research in recent years. This is because periodic signals bring us information about the processes occurring near the innermost regions of blazars, which are mostly inaccessible to our direct view. Such signals provide insights into some of the extreme conditions that take place in the vicinity of supermassive black holes that lead to the launch of the relativistic jets. In addition, studies of characteristic timescales in blazar light curves shed light on some of the challenging issues in blazar physics that includes disk-jet connection, strong gravity near fast rotating supermassive black holes and release of gravitational waves from binary supermassive black hole systems. However, a number of issues associated with the search quasi-periodic oscillations (QPOs) in blazar e.g., red-noise dominance, modest significance of the detection, periodic modulation lasting for only a couple of cycles and their transient nature, make it difficult to estimate the true significance of the detection. Consequently, it also becomes difficult to make meaningful inferences about the nature of the on-going processes. In this proceedings, results of study focused on searching QPOs in a number of blazar multi-frequency light curves are summarized. The time series analyses of long term observations of the blazars revealed presence of year-timescale QPOs in the sources including OJ 287 (optical), Mrk 501 (gamma-ray), J1043+2408 (radio) and PKS 0219-164 (radio). As likely explanations, we discuss a number of scenarios including binary supermassive black hole systems, Lense-Thirring precession, and jet precession.

We present an in-depth and systematic variability study of a sample of 20 powerful blazars, including 12 BL Lacs and 8 flat spectrum radio quasars, applying various analysis tools such as flux distribution, symmetry analysis, and time series analysis on the decade-long Fermi/LAT observations. The results show that blazars with steeper \gama-ray spectral indexes are found to be more variable; and the \gama-ray flux distribution closely resembles lognormal probability distribution function. The statistical variability properties of the sources as studied by power spectral density analysis are consistent with \emph{flicker noise} ($P(\nu)\propto1/\nu$) -- an indication of long-memory processes at work. Statistical analysis of the distribution of flux rise and decay rates in the light curves of the sources, aimed at distinguishing between particle acceleration and energy dissipation timescales, counter-intuitively suggests that both kinds of rates follow a similar distribution and the derived mean variability timescales are in the order of a few weeks. The corresponding emission region size is used to constrain location of \gama-ray production sites in the sources to be a few parsecs. Additionally, using Lomb-Scargle periodogram and weighted wavelet z-transform methods and extensive Monte Carlo simulations, we detected year timescale quasi-periodic oscillations in the sources S5 0716+714, Mrk 421, ON +325, PKS 1424-418 and PKS 2155-304. The detection significance was computed taking proper account of the red-noise and other artifacts inherent in the observations. We explain the results in the light of current blazar models with relativistic shocks propagating down the jet viewed close to the line of sight.

The typical blazar S5 0716$+$714 is very interesting due to its rapid and large amplitude variability and high duty cycle of micro-variability in optical band. We analyze the observations in I, R and V bands obtained with the $1.0m$ telescope at Weihai observatory of Shandong University from 2011 to 2018. The model of synchrotron radiation from turbulent cells in a jet has been proposed as a mechanism for explaining micro-variability seen in blazar light curves. Parameters such as the sizes of turbulent cells, the enhanced particle densities, and the location of the turbulent cells in the jet can be studied using this model. The model predicts a time lag between variations as observed in different frequency bands. Automatic model fitting method for micro-variability is developed, and the fitting results of our multi-frequency micro-variability observations support the model. The results show that both the amplitude and duration of flares decomposed from the micro-variability light curves confirm to the log-normal distribution. The turbulent cell size is within the range of about 5 to 55 AU, and the time lags of the micro-variability flares between the I-R and R-V bands should be several minutes. The time lags obtained from the turbulence model are consistent with the fitting statistical results, and the time lags of flares are correlated with the time lags of the whole light curve.

Propagation of ultra-high energy photons in the solar magnetosphere gives rise to cascades comprising thousands of photons. We study the cascade development using Monte Carlo simulations and find that the photons in the cascades are spatially extended over hundreds of kilometers as they arrive at the top of the Earth's atmosphere. We compare results from simulations which use two models of the solar magnetic field, and show that although signatures of such cascades are different for the models used, for practical detection purpose in the ground-based detectors, they are similar.

Quasi-periodic oscillations (QPO) originating from innermost regions of blazars can provide unique perspective of some of the burning issues in blazar studies including disk-jet connection, launch of relativistic jets from the central engine, and other extreme conditions near the fast rotating supermassive black holes. However, a number of hurdles associated with searching QPOs in blazars e.g., red-noise dominance, modest significance of the detection and periodic modulation lasting for only a couple of cycles, make it difficult to estimate the true significance of the detection. In this work, we report a $\sim$ 330-day QPO in the Fermi/LAT observations of the blazar Mrk 501 spanning nearly a decade. To establish consistency of the result, we adopted multiple approaches to the time series analysis and employed four widely known methods. Among these, Lomb-Scargle periodogram and weighted wavelet z-transform represent frequency domain based methods whereas epoch folding and z-transformed discrete auto-correlation function are time-domain based analysis. Power spectrum response method was followed to properly account for the red-noise, largely inherent in blazar light curves. Both local and global significance of the signal were found to be above 99\% over possible spurious detection. In the context where not many $\gamma$-ray QPOs have been reported to last more than 5 cycles, this might be one of the few instances where we witness a sub-year timescale $\gamma$-ray QPO persisting nearly 7 cycles. A number of possible scenarios linked with binary supermassive black hole, relativistic jets, and accretion disks can be invoked to explain the transient QPO.

Search for periodic signals from blazars has become widely discussed topic in recent years. In the scenario that such periodic changes originate from the innermost regions of blazars, the signals bear imprints of the processes occurring near the central engine, which is mostly inaccessible to our direct view. Such signals provide insights into various aspect of blazar studies including disk-jet connection, magnetic field configuration and, more importantly, strong gravity near the supermassive black holes and release of gravitational waves from the binary supermassive black hole systems. In this work, we report detection of a periodic signal in the radio light curve of the blazar J1043+2408 spanning $\sim$10.5 years. We performed multiple methods of time series analysis, namely, epoch folding, Lomb-Scargle periodogram, and discrete auto-correlation function. All three methods consistently reveal a repeating signal with a periodicity of $\sim560$ days. To robustly account for the red-noise processes usually dominant in the blazar variability and other possible artifacts, a large number of Monte Carlo simulations were performed. This allowed us to estimate a high significance (99.9\% local and 99.4\% global) against possible spurious detection. As possible explanations, we discuss a number of scenarios including binary supermassive black hole system, Lense-Thirring precession and jet precession.

In this presentation, we report the results of intraday variability in the optical (BVRI bands) and hard X-ray band (3-79 keV) in a number of blazars. In the optical microvariability studies of the blazars S5 0716+714 and BL Lac, we observed many interesting features such as rapid variability, large variability amplitude, presence of characteristic timescales, bluer-when-brighter achromatic behavior, and single power-law power spectral density. In {\it NuSTAR} observations of several blazars, using spectral and timing analysis, we found similar features consistent with the optical studies. In addition, in BL Lacs we estimated the Lorentz factor of the population of highest energy electrons emitting synchrotron emission, and whereas in flat-spectrum radio quasars, using external Compton models, we estimated the energy of the lower end of the injected electrons to be a few tens of Lorentz factors. In addition, we find that the low flux state exhibit more rapid variability in contrast to the previously reported results showing high flux states displaying rapid variability. In both the studies, the size of the emission regions estimated using variability timescales turn out to be an order magnitude smaller than the gravitational radius of a typical black-hole masses between $\sim 10^8-10^9$ solar masses, believed to be at the center of the radio-loud AGN. The results of the studies suggest that these low-amplitude rapid variability might originate as a result of magnetohydrodynamical instabilities near the base of the jets triggered by the processes modulated by the magnetic field at accretion disc.

Here we explore the disk-jet connection in the broad-line radio quasar 4C+74.26, utilizing the results of the multiwavelength monitoring of the source. The target is unique in that its radiative output at radio wavelengths is dominated by a moderately-beamed nuclear jet, at optical frequencies by the accretion disk, and in the hard X-ray range by the disk corona. Our analysis reveals a correlation (local and global significance of 96\% and 98\%, respectively) between the optical and radio bands, with the disk lagging behind the jet by $250 \pm 42$ days. We discuss the possible explanation for this, speculating that the observed disk and the jet flux changes are generated by magnetic fluctuations originating within the innermost parts of a truncated disk, and that the lag is related to a delayed radiative response of the disk when compared with the propagation timescale of magnetic perturbations along relativistic outflow. This scenario is supported by the re-analysis of the NuSTAR data, modelled in terms of a relativistic reflection from the disk illuminated by the coronal emission, which returns the inner disk radius $R_{\rm in}/R_{\rm ISCO} =35^{+40}_{-16}$. We discuss the global energetics in the system, arguing that while the accretion proceeds at the Eddington rate, with the accretion-related bolometric luminosity $L_{\rm bol} \sim 9 \times 10^{46}$ erg s$^{-1}$ $\sim 0.2 L_{\rm Edd}$, the jet total kinetic energy $L_\textrm{j} \sim 4 \times 10^{44}$ erg s$^{-1}$, inferred from the dynamical modelling of the giant radio lobes in the source, constitutes only a small fraction of the available accretion power.

Context. Investigation of the hard X-ray emission properties of blazars is key to the understanding of the central engine of the sources and associated jet process. In particular, simultaneous spectral and timing analyses of the intra-day hard X-ray observations provide us a means to peer into the compact innermost blazar regions, not accessible to our current instruments. Aims. The primary objective of the work is to associate the observed hard X-ray variability properties in the blazars to their flux and spectral states, thereby, based on the correlation among them, extract the details about the emission regions and the processes occurring near the central engine. Methods. We carried out timing, spectral and cross-correlation analysis of 31 NuSTAR observations of 13 blazars. We investigated the spectral shapes of the sources using single power-law, broken power-law and log-parabola models. We also studies the co-relation between the soft and the hard emission using z-transformed discrete correlation function. Results. We found that for most of the sources the hard X-ray emission can be well represented by log-parabola model; and that the spectral slopes for different blazar sub-classes are consistent with so called blazar sequence. We noted a close connection between the flux and spectral slope within the source sub-class in the sense that high flux and/or flux states tend to be harder in spectra. In BL Lacertae objects, assuming particle acceleration by diffusive shocks and synchrotron cooling as the dominant processes governing the observed flux variability, we constrain the magnetic field of the emission region to be a few gauss; whereas in flat-spectrum radio quasars, using external Compton models, we estimate the energy of the lower end of the injected electrons to be a few Lorentz factors.

Although the photon structure is most efficiently studied with the accelerator instruments, there is also a scientifically complementary potential in investigations on photons produced in the outer space. This potential is already being explored with gamma ray telescopes, ultra-high energy cosmic ray observatories and, since very recently, by the Cosmic-Ray Extremely Distributed Observatory (CREDO). Unlike the former instruments focused on detection of single photons, CREDO aims at the detection of cascades (ensembles) of photons originating even at astrophysical distances. If at least a part of such a cascade reaches Earth, it might produce a unique pattern composed of a number of air showers observable by an appropriately dense array of standard detectors. If the energies of air showers constituting the pattern are relatively low and if the typical distances between the neighbors are large, the ensemble character of the whole phenomenon might remain uncovered, unless the CREDO strategy is implemented.

In this work, we present the results of our multi-band microvariability study of the famous blazar BL Lac. We performed microvariablity observations of the source in the optical VRI bands for 4 nights in the year 2016. We studied the intranight flux and spectral variability of the source in detail with an objective to characterize microvariability in the blazars, a frequently observed phenomenon in blazars. The results show that the source often displays a fast flux variability with an amplitude as large as ~ 0.2 magnitude within a few hours, and that the color variability in the similar time scales can be characterized as "bluer-when-brighter" trend. We also observed markedly curved optical spectrum during one of the nights. Furthermore, the correlation between multi-band emission shows that in general the emission in all the bands are highly correlated; and in one of the nights V band emission was found to lead the I band emission by ~ 20 minutes. The search for characteristic timescale using auto-correlation function and the structure function analyses reveals characteristic timescale of ~ 48 minutes in one of the R band observations. We try to explain the observed results in the context of the passage of shock waves through the relativistic outflows in blazars.

In this work, we explore the long-term variability properties of the blazar PKS 0219-164 in the radio and the $\gamma$-ray regime, utilizing the OVRO 15 GHz and the \fermi/LAT observations from the period 2008--2017. We found that $\gamma$-ray emission is more variable than the radio emission implying that $\gamma$ ray emission possibly originated in more compact regions while the radio emission represented continuum emission from the large scale jets. Also, in $\gamma$-ray the source exhibited spectral variability characterized by the \emph{ softer-when-brighter} trend, a less frequently observed feature in the high energy emission by BL Lacs. In radio, using Lomb-Scargle periodogram and weighted wavelet z-transform, we detected a strong signal of quasi-periodic oscillation (QPO) with a periodicity of 270 $\pm$ 26 days with possible harmonics of 550 $\pm$ 42 and 1150 $\pm$ 157 days periods. At a time when the issue of the detection of QPOs in blazars is still under debate, the observed QPO with high statistical significance ( $\sim$ 97\% -- 99\% global significance over underlying red-noise processes) and persistent over nearly 10 oscillations could make one of the strongest cases for the detection of QPOs in blazar light curves. We discuss various blazar models that might lead to the $\gamma$-ray and radio variability, QPO, and the achromatic behavior seen in the high energy emission from the source.

Detection of periodicity in the broad-band non-thermal emission of blazars has so far been proven to be elusive. However, there are a number of scenarios which could lead to quasi-periodic variations in blazar light curves. For example, orbital or thermal/viscous period of accreting matter around central supermassive black holes could, in principle, be imprinted in the multi-wavelength emission of small-scale blazar jets, carrying as such crucial information about plasma conditions within the jet launching regions. In this paper, we present the results of our time series analysis of $\sim 9.2$ year-long, and exceptionally well-sampled optical light curve of the BL Lac OJ 287. The study primarily uses the data from our own observations performed at the Mt. Suhora and Krak\'ow Observatories in Poland, and at the Athens Observatory in Greece. Additionally, SMARTS observations were used to fill in some of the gaps in the data. The Lomb-Scargle Periodogram and the Weighted Wavelet Z-transform methods were employed to search for the possible QPOs in the resulting optical light curve of the source. Both the methods consistently yielded possible quasi-periodic signal around the periods of $\sim 400$ and $\sim 800$ days, the former one with a significance (over the underlying colored noise) of $\geq 99\%$. A number of likely explanations for such are discussed, with a preference given to a modulation of the jet production efficiency by highly magnetized accretion disks. This supports the previous findings and the interpretation reported recently in the literature for OJ 287 and other blazar sources.

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.

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.

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.