Determining whether the flux distribution of an Astrophysical source is a Gaussian or a log-normal, provides key insight into the nature of its variability. For lightcurves of moderate length ($< 10^3$), a useful first analysis is to test the Gaussianity of the flux and logarithm of the flux, by estimating the skewness and applying the Anderson-Darling (AD) method. We perform extensive simulations of lightcurves with different lengths, variability, Gaussian measurement errors and power spectrum index $\beta$ (i.e. $P(f) \propto f^{-\beta}$), to provide a prescriptionand guidelines for reliable use of these two tests. We present empirical fits for the expected standard deviation of skewness and tabulated AD test critical values for $\beta = 0.5$ and $1.0$, which differ from the values given in the literature which are for white noise ($\beta = 0$). Moreover, we show that for white noise, for most practical situations, these tests are meaningless, since binning in time alters the flux distribution. For $\beta \gtrsim 1.5$, the skewness variance does not decrease with length and hence the tests are not reliable. Thus, such tests can be applied only to systems with $\beta\gtrsim 0.5$ and $\beta \lesssim 1.0$. As an example of the prescription given in this work, we reconfirm that the Fermi data of the blazar, 3FGL\,J0730.2-1141, shows that its $\gamma$ ray flux is consistent with a log-normal distribution and not with a Gaussian one.

We present a detailed study of flux and index distributions of three blazars (one FSRQ and two BL Lacs) by using 16 years of Rossi X-ray Timing Explorer archival data. The three blazars were chosen such that their flux and index distributions have sufficient number of data points ($\geq$90) with relatively less uncertainty ($\bar{\sigma_{err}^{2}}/\sigma^{2}$<0.2) in light curves. Anderson-Darling (AD) test and histogram fitting shows that flux distribution of FSRQ 3C273 is log-normal, while its photon index distribution is Gaussian. This result is consistent with linear Gaussian perturbation in the particle acceleration time-scale, which produces log-normal distribution in flux. However, for two BL Lacs viz. Mkn501 and Mkn421, AD test shows that their flux distributions are neither Gaussian nor log-normal, and their index distributions are non-normal. The histogram fitting of Mkn501 and Mkn421, suggests that their flux distributions are more likely to be a bi-modal, and their index distributions are double Gaussian. Since, Sinha et al. (2018) had shown that Gaussian distribution of index produces a log-normal distribution in flux, double Gaussian distribution of index in Mkn501 and Mkn421 indicates that their flux distributions are probably double log-normal. Observation of double log-normal flux distribution with double Gaussian distribution in index reaffirms two flux states hypothesis. Further, the difference observed in the flux distribution of FSRQ (3C273) and BL Lacs (Mkn501 and Mkn421) at X-rays, suggest that the low energy emitting electrons have a single log-normal flux distribution while the high energy ones have a double log-normal flux distribution.

We present a detailed temporal and spectral study of the blazar 3C\,279 using multi-wavelength observations from Swift-XRT, Swift-UVOT and Fermi-LAT during a flare in 2018 January. The temporal analysis of $\gamma$-ray light curve indicates a lag of $\sim 1$ d between the 0.1--3 GeV and 3--500 GeV emission. Additionally, the $\gamma$-ray light curve shows asymmetry with slow rise--fast decay in energy band 0.1--3 GeV and fast rise--slow decay in the 3--500 GeV band. We interpret this asymmetry as a result of shift in the Compton spectral peak. This inference is further supported by the correlation studies between the flux and the parameters of the log-parabola fit to the source spectra in the energy range 0.1--500 GeV. We found that the flux correlates well with the peak spectral energy and the log-parabola fit parameters show a hard index with large curvature at high flux states. Interestingly, the hardest index with large curvature was synchronous with a very high energy flare detected by H.E.S.S. Our study of the spectral behavior of the source suggests that $\gamma$-ray emission is most likely to be associated with the Compton up-scattering of IR photons from the dusty environment. Moreover, the fit parameters indicate the increase in bulk Lorentz factor of emission region to be a dominant cause for the flux enhancement.

We present the results of the $\gamma$-ray flux distribution study on the brightest blazars which are observed by the \emph{Fermi}-LAT. We selected 50 brightest blazars based on the maximum number of detection reported in the LAT third AGN catalog. We performed standard unbinned maximum likelihood analysis on the LAT data during the period between August 2008 and December 2016, in order to obtain the average monthly flux. After quality cuts, blazars for which at least 90\% of the total flux was survived were selected for the further study, and this includes 19 FSRQs and 19 BL Lacs. The Anderson-Darling and $\chi^2$ tests suggest that the integrated monthly flux follow a log-normal distribution for all sources, except for three FSRQs for which neither a normal nor a log-normal distribution was preferred. A double log-normal flux distribution tendency were observed in these sources, though it has to be confirmed with improved statistics. We also found that, the standard deviation of the log-normal flux distribution increases with the mean spectral index of the blazar, and can be fitted with a line of slope 0.24$\pm$0.04. We repeat our study on three additional brightest unclassified blazars to identify their flux distribution properties. Based on the features of their log-normal flux distribution, we infer these unclassified blazars may be closely associated with FSRQs. We also highlight that considering the log-normal behavior of the flux distribution of blazars, averaging their long term flux in linear scale can largely under estimate the nominal flux and this discrepancy can propagate down to the estimation of source parameters through spectral modeling.

We perform a detailed spectral study of a recent flaring activity from the Flat Spectrum Radio Quasar (FSRQ), 3C\,454.3, observed simultaneously in optical, UV, X-ray and $\gamma$-ray energies during 16 to 28 August, 2015. The source reached its peak $\gamma$-ray flux of $(1.9\pm0.2)\times\,10^{-05} \; {\rm ph\,cm^{-2}\,s^{-1}}$ on 22 August. The time averaged broadband spectral energy distribution (SED) is obtained for three time periods, namely "flaring state"; covering the peak $\gamma$-ray flux, "post flaring state"; immediately following the peak flare and "quiescent state"; separated from the flaring event and following the post flaring state. The SED corresponding to the flaring state is investigated using different emission models involving synchrotron, synchrotron self Compton (SSC) and external Compton (EC) mechanisms. Our study suggests that the X-ray and $\gamma$-ray emission from 3C\,454.3 cannot be attributed to a single emission mechanism and instead, one needs to consider both SSC and EC mechanisms. Moreover, the target photon energy responsible for the EC process corresponds to an equivalent temperature of 564 K, suggesting that the flare location lies beyond the broad line emitting region of the FSRQ. SED fitting of the other two flux states further supports these inferences.