We used the data from the Fermi Gamma-ray Space Telescope to characterise the {\gamma}-ray flux variability of blazars on month-like time scales. Our sample consists of 1120 blazars of which 481 are flat spectrum radio quasars (FSRQs) and 639 are BL Lac objects (BL Lacs). We generated monthly binned light curves of our sample for a period of approximately nine years from 2008 August to 2017 December and quantified variability by using excess variance (F_var). On month-like time scales, 371/481 FSRQs are variable ( 80%), while only about 50% (304/639) of BL Lacs are variable. This suggests that FSRQs are more variable than BL Lac objects. We find a mean F_var of 0.55 +/- 0.33 and 0.47 +/- 0.29 for FSRQs and BL Lacs respectively. Large F_var in FSRQs is also confirmed from the analysis of the ensemble structure function. By Dividing our sample of blazars based on the position of the synchrotron peak in their broad-band spectral energy distribution, we find that the low synchrotron peaked (LSP) sources have the largest mean F_var value of 0.54 +/- 0.32 while the intermediate synchrotron peaked (ISP) and high synchrotron peaked (HSP) sources have mean F_var values of 0.45 +/- 0.25 and 0.47 +/- 0.33 respectively. On month-like time scales, we find FSRQs to show a high duty cycle (DC) of variability of 66% relative to BL Lacs that show a DC of 36%. We find that both the F_var and time scale of variability ({\tau}) do not correlate with the black hole mass. We note that F_var is found to be weakly correlated with Doppler factor ({\delta}) and {\tau} is also weakly correlated with {\delta}. Most of the sources in our sample have {\tau} of the order of days, which might be related to processes in the jet. We find marginal difference in the distribution of {\tau} between FSRQs and BL Lacs.

Blazars show optical and $\gamma$-ray flux variations that are generally correlated, although there are exceptions. Here we present anomalous behaviour seen in the blazar 3C 454.3 based on an analysis of quasi-simultaneous data at optical, UV, X-ray and $\gamma$-ray energies, spanning about 9 years from August 2008 to February 2017. We have identified four time intervals (epochs), A, B, D and E, when the source showed large-amplitude optical flares. In epochs A and B the optical and $\gamma$-ray flares are correlated, while in D and E corresponding flares in $\gamma$-rays are weak or absent. In epoch B the degree of optical polarization strongly correlates with changes in optical flux during a short-duration optical flare superimposed on one of long duration. In epoch E the optical flux and degree of polarization are anti-correlated during both the rising and declining phases of the optical flare. We carried out broad-band spectral energy distribution (SED) modeling of the source for the flaring epochs A,B, D and E, and a quiescent epoch, C. Our SED modeling indicates that optical flares with absent or weak corresponding $\gamma$-ray flares in epochs D and E could arise from changes in a combination of parameters, such as the bulk Lorentz factor, magnetic field and electron energy density, or be due to changes in the location of the $\gamma$-ray emitting regions.

We present the results of our broadband study of the $\gamma$-ray emitting narrow line Seyfert 1 (NLSy1) galaxy SBS 0846+513 ($z=0.585$). This includes multi-band flux variations, $\gamma$-ray spectral analysis, broad band spectral energy distribution (SED) modeling, and intranight optical variability (INOV) observations carried over 6 nights between 2012 November and 2013 March using the 2 m Himalayan Chandra Telescope and the 1.3 m telescope at Devasthal. Multiple episodes of flaring activity are seen in the $\gamma$-ray light curve of the source which are also reflected in the observations at lower frequencies. A statistically significant curvature is noticed in the seven years averaged $\gamma$-ray spectrum, thus indicating its similarity with powerful flat spectrum radio quasars (FSRQs). Modeling the SEDs with a one zone leptonic emission model hints the optical-UV spectrum to be dominated by synchrotron radiation, whereas, inverse Compton scattering of broad line region photons reproduces the $\gamma$-ray part of the SEDs. The source was found to be variable on all the six nights of optical observations with a variation of $\sim$0.3 magnitude within a single night, coinciding with a high $\gamma$-ray activity state. The observed large amplitude INOV clearly indicates the presence of a closely aligned beamed relativistic jet in SBS 0846+513. Our broadband study supports the recent claims in literature that $\gamma$-ray emitting NLSy1 galaxies are similar to blazars and constitute the low black hole mass counterparts to FSRQs.