3C 454.3 is frequently observed in flaring state. The long term light curve of this source has been analysed with 9 year (August 2008 - July 2017) data from Fermi LAT detector. We have identified five flares and one quiescent state. The flares have sub-structures with many peaks during flaring phase. We have estimated the rise and decay time of the flares and compared with flares of other similar sources. The modeling of gamma ray spectral energy distributions shows in most cases Log parabola function gives the best fit to the data. We have done time dependent leptonic modeling of two of the flares, for which simultaneous multi-wavelength data are available. These two long lasting flares Flare-2A and Flare-2D continued for 95 days and 133 days respectively. We have used the average values of Doppler factor, injected luminosity in electrons, size of the emission region and the magnetic field in the emission region in modeling these flares. The emission region is assumed to be in the broad line region in our single zone model. The energy losses (synchrotron, synchrotron self-Compton, external Compton) and escape of electrons from the emission region have been included while doing the modelling. Although, the total jet powers required to model these flares with leptonic model are higher compared to other sources, they are always found to be lower than the Eddington's luminosity of 3C 454.3. We also select some flaring peaks and show that time variation of the Doppler factor or the injected luminosity in electrons over short time scales can explain their light curves.

The first ever identification of a cosmic ray accelerator as the consequence of spacial and temporal correlation of IceCube event 170922A with flaring of a blazar TXS 0506+056 motivated us to look for other flaring blazars in Fermi-LAT 3FGL catalog, which could be correlated with IceCube high energy track events. We have studied the Fermi-LAT light curves of blazars correlated with neutrino track events. Among the eight sources identified within 2$\sigma$ angular uncertainty of the IceCube track events selected in our study, we find only one source 3FGL J2255+2409 was in flaring state during the neutrino detection. We have carried out a time dependent modelling of the multi-wavelength data from this blazar, and the neutrino event including leptonic energy losses and proton-proton interactions in its jet to determine whether it could be the origin of the neutrino event. Our lepto-hadronic model estimates a jet luminosity of $L_j = 3.6 \times10^{47}$ erg/sec during the neutrino phase of 3FGL/4FGL J2255+2409.

A multiwavelength temporal and spectral analysis of flares of 3C 279 during November 2017--July 2018 are presented in this work. Three bright gamma-ray flares were observed simultaneously in X-ray and Optical/UV along with a prolonged quiescent state. A "harder-when-brighter" trend is observed in both gamma-rays and X-rays during the flaring period. The gamma-ray light curve for all the flares are binned in one-day time bins and a day scale variability is observed. Variability time constrains the size and location of the emission region to 2.1$\times$10$^{16}$ cm and 4.4$\times$10$^{17}$ cm, respectively. The fractional variability reveals that the source is more than 100\% variable in gamma-rays and it decreases towards the lower energy. A cross-correlation study of the emission from different wavebands is done using the \textit{DCF} method, which shows a strong correlation between them without any time lags. The zero time lag between different wavebands suggest their co-spatial origin. This is the first time 3C 279 has shown a strong correlation between gamma-rays and X-rays emission with zero time lag. A single zone emission model was adopted to model the multiwavelength SEDs by using the publicly available code GAMERA. The study reveals that a higher jet power in electrons is required to explain the gamma-ray flux during the flaring state, as much as, ten times of that required for the quiescent state. However, more jet power in magnetic field has been observed during the quiescent state compared to the flaring state.

PKS 1510-089 is one of the most variable blazars in the third Fermi-LAT source catalog. During 2015, this source has shown four flares identified as flare A, B, C, and D in between three quiescent states Q1, Q2, and Q3. The multi-wavelength data from Fermi-LAT, Swift-XRT/UVOT, OVRO, and SMA observatory are used in our work to model these states. Different flux doubling times have been observed in different energy bands which indicate there could be multiple emission zones. The flux doubling time from the gamma-ray and X-ray light curves are found to be 10.6 hr, 2.5 days, and the average flux doubling time in the optical/UV band is 1 day. It is possible that the gamma-ray and optical/UV emission are produced in the same region whereas X-ray emission is coming from a different region along the jet axis. We have also estimated the discrete correlations functions (DCFs) among the light curves of different energy bands to infer about their emission regions. However, our DCF analysis does not show significant correlation in different energy bands though it shows peaks in some cases at small time lags. We perform a two-zone multi-wavelength time-dependent SED modeling with one emission zone located near the outer edge of the broad line region (BLR) and another further away in the dusty/molecular torus (DT/MT) region to study this high state.

In this work, I have presented a multi-frequency variability and correlation study of the blazar Ton 599, which was observed first time in flaring state at the end of 2017. Data from \textit{Fermi}-LAT, Swift-XRT/UVOT, Steward Observatory, and OVRO (15 GHz) is used, and it is found that the source is more variable in $\gamma$-ray and optical/UV than X-ray and radio. Large variations in the degree of polarization (DoP) and position angle (PA) is noticed during the flaring period. Maximum flux during $\gamma$-ray flare is found to be 12.63$\times$10$^{-7}$ at MJD 58057.5 from the 1-day bin light curve (LC), which is the maximum flux ever achieved by this source. It is further found that all the peaks of flare are very symmetric, which suggests the cooling time of electrons is much smaller than light crossing time. Using 1-day as a fast variability time, the size of the $\gamma$-ray emission region is estimated as 1.88$\times$10$^{16}$ cm. Two 42 GeV of photons are detected during the flare which puts a constraint on the location of the emission region, and it is found that the $\gamma$-ray emitting blob is located at the outer edge or outside the broad line region (BLR). A trend of increasing fractional variability towards higher energies is also seen. Strong correlations were seen between $\gamma$-ray, optical/UV, X-ray, and radio (15 GHz) emission. A small time lag between $\gamma$-ray and optical/UV suggest their emission to be co-spatial while lag of 27 days between $\gamma$-ray and OVRO (15 GHz) suggest two different emission zone separated by a distance of $\sim$ 5 pc.

We present a multi-wavelength study of the FSRQ CTA 102 using Fermi-LAT and simultaneous Swift-XRT/UVOT observations. The Fermi-LAT telescope detected one of the brightest flares from this object during Sep, 2016 to Mar, 2017. In the 190 days of observation period the source underwent four major flares. A detailed analysis of the temporal and spectral properties of these flares indicates the flare at MJD 57751.594 has a $\gamma$-ray flux of (30.12$\pm$4.48)$\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ (from 90 minutes binning) in the energy range of 0.1--300 GeV. This has been found to be the highest flux ever detected from CTA 102. Time dependent leptonic modelling of the pre-flare, rising state, flares and decaying state has been done. A single emission region of size $6.5\times 10^{16}$ cm has been used in our work to explain the multi-wavelength spectral energy distributions. During flares the luminosity in electrons increases nearly seventy times compared to the pre-flare state.

We have analyzed data from the Flat Spectrum Radio Quasar PKS 1510-089 collected over a period of 8 years from August 2008 to December 2016 with the Fermi-LAT. We have identified several flares of this highly variable source, studied their temporal and spectral properties in detail and compared with previous works on flares of PKS 1510-089. Five major flares and few sub-flares/sub-structures have been identified in our study. The fastest variability time is found to be 1.30$\pm$0.18 hr between MJD 55852.063 and 55852.188 where we estimate the minimum size of the emission region to be $4.85 \times 10^{15}$ cm. In most of the flares the spectral energy distributions are better fitted with Logparabolic distribution compared to simple Power law or Power law with exponential cut-offs. This has strong physics implications regarding the nature of the high energy gamma-ray emission region.