We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

With an aim to explore optical variability at diverse timescales in BL Lac source S5 0716+714, it was observed for 46 nights during 2013 January 14 to 2015 June 01 when it underwent two major outbursts. The observations were made using the 1.2-m Mount Abu InfraRed Observatory telescope mounted with a CCD camera. On 29 nights, the source was monitored for more than two hours, resulting in 6256 data points in R-band, to check for the intra-night variability. Observations in B, V and I bands with 159, 214, and 177 data points, respectively, along with daily averaged R-band data are used to address inter-night and long-term variability and the color behavior of S5 0716+71. The study suggests that the source shows significant intra-night variability with a duty cycle of more than 31% and night-to-night variations. The average brightness magnitudes in B, V, R & I bands were found to be 14.42(0.02), 14.02(0.01), 13.22(0.01) & 13.02(0.03), respectively, while S5 0716+714 was historically brightest with R = 11.68 mag on 2015 January 18, indicating that source was in relatively high state during this period. A mild bluer when brighter behavior, typical of BL Lacs, supports the shock-in-jet model. We notice larger amplitudes of variation when the source was relatively brighter. Based on the shortest time scale of variability and causality argument, upper bound on the size of the emission region is estimated to be 9:32x 10^{14} cm and the mass of the black hole to be 5.6 x10^8Ms.

Blazars outburst provide an opportunity to study their spectral behaviour, correlated variations, and structure of the jet. Such an unprecedented flaring activity in FSRQ CTA102, during 2016 - 2017, is used for a detailed study to understand variability mechanisms at short and long time scales, spectral behaviour in different energy regimes and to estimate sizes and location of the high-energy emitting region in the jet. Multiwavelength (MW) data (gamma-ray, X-ray, UV, optical and radio) for CTA102 during its outburst period, are obtained from Fermi-LAT, Swift-XRT/UVOT, Steward, Mt Abu Infrared Observatory and OVRO. These are analyzed to construct MW light curves, extract the spectral information, and to perform the correlated variability studies. Our study shows that CTA102 attained the highest ever flux levels across the EMS while flaring and otherwise, with rapid and prolonged activities at all the frequencies. A number of short term (3 to 8 days) and long term (>month) variability events are noticed across the EMS. We infer a redder when brighter trend in faint state and a bluer when brighter trend during a few optical flares. Based on the flux doubling timescale, the size of the gamma-ray emitting region is estimated as 8.76x10^{15}cm, located at a distance of about 5.58x10^{16} cm from the central engine. CTA102 was in extremely bright phase during 2016-17, due possibly to, successive high energy particle injections into the jet, creating shocks traveling down the jet which lead to the overall flux enhancement across the EMS. Alternatively, a decreasing viewing angle could also lead to the enhanced flux. The study reveals correlated variations in all the energy bands, with lags within time bins, indicating to co-spatial origin of the emissions. During the flaring event, a BWB color in optical and harder when brighter trend in X-ray and gamma-ray spectra are noticed.

PKS1510-089 is a flat spectrum radio quasar located at a redshift of 0.36. It is one of only a few such sources detected in very-high-energy (VHE, >100 GeV) gamma rays. Though PKS1510-089 is highly variable at GeV energies, until recently no variability has been observed in the VHE band. In 2015 May PKS1510-089 showed a high state in optical and in the GeV range. A VHE gamma-ray flare was detected with MAGIC at that time, showing the first instance of VHE gamma-ray flux variability on the time scale of days in this source. We will present the MAGIC results from this observation, discuss their temporal and spectral properties in the multi-wavelength context and present modelling of such emission in the external Compton scenario.

We present the results from a multiwavelength study of the flaring activity in HBL, 1ES 1959+650, during January 2015-June 2016. The source underwent significant flux enhancements showing two major outbursts (March 2015 and October 2015) in optical, UV, X-rays and gamma-rays. Normally, HBLs are not very active but 1ES 1959+650 has shown exceptional outburst activity across the whole electromagnetic spectrum (EMS). We used the data from Fermi-LAT, Swift-XRT & UVOT and optical data from Mt. Abu InfraRed Observatory (MIRO) along with archival data from Steward Observatory to look for possible connections between emissions at different energies and the nature of variability during flaring state. During October 2015 outburst, thirteen nights of optical follow-up observations showed brightest and the faintest nightly averaged V-band magnitudes as 14.45(0.03) and 14.85(0.02), respectively. In optical, the source showed a hint of optical intra-night variability during the outburst. A significant short-term variability in optical during MJD 57344 to MJD 57365 and in gamma-rays during MJD 57360 and MJD 57365 was also noticed. Multiwavelength study suggests the flaring activity at all frequencies to be correlated in general, albeit with diverse flare durations. We estimated the strength of the magnetic field as 4.21 G using the time-lag between optical and UV bands as synchrotron cooling time scale (2.34 hrs). The upper limits on the sizes of both the emission regions, gamma-ray and optical, are estimated to be of the order of 10^16cm using shortest variability time scales. The quasi-simultaneous flux enhancements in 15 GHz and VHE gamma-ray emissions indicates to a fresh injection of plasma into the jet, which interacts with a standing sub-mm core resulting in co-spatial emissions across the EMS. The complex and prolonged behavior of the second outburst in October 2015 is discussed in detail.

We present results based on the observations of the blazar 3C 66A from 2005 November 06 to 2016 February 14 in the BVR and I broadbands using 1.2m telescope of the Mt. Abu InfraRed Observatory (MIRO). The source was observed on 160 nights out of which on 89 nights it was monitored for more than 1 hr to check for the presence of any intra-day variability (IDV). The blazar 3C 66A exhibited significant variations in the optical flux on short and long term time scales. However, unlike highly variable S5 0716+71, it showed IDV duty cycle of about 8% only. Our statistical studies suggest the IDV time scales ranging from $\sim$ 37 min to about 3.12 hours and, in one case, a possibility of the quasi-periodic variations with characteristic timescale of $\sim$ 1.4 hrs. The IDV amplitudes in R$-$band were found to vary from 0.03 mag to as much as 0.6 mag, with larger amplitude of variation when the source was relatively fainter. The typical rate of the flux variation was estimated to be $\sim$0.07 mag hr$^{-1}$ in both, the rising and the falling phases. However, the rates of the brightness variation as high as 1.38 mag/hr were also detected. The shortest timescale of the variation, 37 min, sets an upper limit of $6.92 \times 10^{14}$ cm on the size of the emission region and about $3.7 \times 10^8\ \mathrm{M}_{\odot}$ as an estimate of the mass of the black hole, assuming the origin of the rapid optical variability is in close vicinity of the central SMBH. The long-term study suggests a mild bluer-when-brighter behavior, typical for BL Lacs.

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.

We present a detailed investigation of the flaring activity observed from a BL Lac object, S5 0716+714 , during its brightest ever optical state in the second half of January 2015. Observed almost simultaneously in the optical, X-rays and {\gamma}-rays, a significant change in the degree of optical polarization (PD) and a swing in the position angle (PA) of polarization were recorded. A detection in the TeV (VHE) was also reported by the MAGIC consortium during this flaring episode. Two prominent sub-flares, peaking about 5-days apart, were seen in almost all the energy bands. The multi-wavelength light-curves, spectral energy distribution (SED) and polarization are modeled using the time-dependent code developed by Zhang et al. (2014). This model assumes a straight jet threaded by large scale helical magnetic fields taking into account the light travel time effects, incorporating synchrotron flux and polarization in 3D geometry. The rapid variation in PD and rotation in PA are most likely due to re-connections happening in the emission region in the jet, as suggested by the change in the ratio of toroidal to poloidal components of magnetic field during quiescent and flaring states.