In this work, we have analyzed the $\gamma$-ray data in the energy range 100 MeV to 500 GeV from the \emph{Fermi}-Large Area Telescope (LAT) observations for the period August 4, 2008 to March 31, 2019. The $\gamma$-ray emission from AR Scorpii over the last decade is not statistically significant and therefore 2$\sigma$ upper limit on the integral flux above 100 MeV has been estimated. We reproduce the non-thermal broadband spectral energy distribution of AR Scorpii using an emission model having two synchrotron components due to the relativistic electrons in very high magnetic fields. The first component (Synchrotron-1) broadly describes the emissions at radio to high energy X-rays through the synchrotron radiation originating from a spherical region of radius $\sim$ 1.8$\times$10$^{10}$ cm and a magnetic field strength of $\sim$ 10$^3$ Gauss. The second component (Synchrotron-2) which reproduces the X-ray emission at lower energies and predicts the $\gamma$-ray emission, originates from another spherical region with radius $\sim$ 1.4$\times$10$^{10}$ cm and a magnetic field strength of $\sim$ 10$^6$ Gauss. The relativistic electron populations in both the emission regions are described by a smooth broken power law energy distribution. The $\gamma$-ray emission predicted by the Synchrotron-2 model is below the broadband sensitivity of the \emph{Fermi}-LAT and is also consistent with the 95$\%$ confidence level upper limit on the integral flux above 100 MeV derived from more than 10 years of observations. According to our model, the binary system AR Scorpii could be a $\gamma$-ray source, although its emission level must be below the current detection limit of the \emph{Fermi}-LAT.

We report the analysis of simultaneous multi-wavelength data of the high energy peaked blazar RGB J0710+591 from the LAXPC, SXT and UVIT instruments on-board AstroSat. The wide band X-ray spectrum (0.35 -- 30 keV) is modelled as synchrotron emission from a non-thermal distribution of high energy electrons. The spectrum is unusually curved, with a curvature parameter $\beta_p \sim 6.4$ for a log parabola particle distribution, or a high energy spectral index $p_2 > 4.5$ for a broken power-law distribution. The spectrum shows more curvature than an earlier quasi-simultaneous analysis of Swift-XRT/NuSTAR data where the parameters were $\beta_p \sim 2.2$ or $p_2 \sim 4$. It has long been known that a power-law electron distribution can be produced from a region where particles are accelerated under Fermi process and the radiative losses in acceleration site decide the maximum attainable Lorentz factor, $\gamma_{max}$. Consequently, this quantity decides the energy at which the spectrum curves steeply. We show that such a distribution provides a more natural explanation for the AstroSat data as well as the earlier XRT/NuSTAR observation, making this as the first well constrained determination of the photon energy corresponding to $\gamma_{max}$. This in turn provides an estimate of the acceleration time-scale as a function of magnetic field and Doppler factor. The UVIT observations are consistent with earlier optical/UV measurements and reconfirm that they plausibly correspond to a different radiative component than the one responsible for the X-ray emission.

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 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.

The X-ray spectral curvature of blazars is traditionally explained by an empirical log-parabola function characterized by three parameters, namely the flux, curvature and spectral index at a given energy. Since their exact relationship with the underlying physical quantities is unclear, interpreting the physical scenario of the source through these parameters is difficult. To attain an insight on the X-ray spectral shape, we perform a detailed study of the X-ray spectra of the blazar MKN 421, using an analytical model where the electron diffusion from the particle acceleration site is energy-dependent. The resultant synchrotron spectrum is again determined by three parameters, namely, the energy index of the escape time scale, the quantity connecting the electron energy to the observed photon energy and the normalization. The X-ray observations of MKN 421, during July 2012 - April 2013 by NuSTAR and Swift-XRT are investigated using this model and we find a significant correlation between model parameters and the observational quantities. Additionally, a strong anti-correlation is found between the fit parameters defining the spectral shape, which was not evident from earlier studies using empirical models. This indicates the flux variations in MKN 421 and possibly other blazars, may arise from a definite physical process that needs to be further investigated.

The observed log-normal flux distributions in the high energy emission from blazars have been interpreted as being due to variability stemming from non-linear multiplicative processes generated dynamically from the accretion disc. On the other hand, rapid minute scale variations in the flux point to a compact emitting region inside the jet, probably disconnected from the disc. In this work, we show that linear Gaussian variations of the intrinsic particle acceleration or escape time-scales can produce distinct non-Gaussian flux distributions, including log-normal ones. Moreover, the spectral index distributions can provide confirming evidence for the origin of the variability. Thus, modelling of the flux and index distributions can lead to quantitative identification of the micro-physical origin of the variability in these sources. As an example, we model the X-ray flux and index distribution of Mkn 421 obtained from over 9 years of MAXI observations and show that the variability in the X-ray emission is driven by Gaussian fluctuations of the particle acceleration process rather than that of the escape rate.

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.

The broad band optical/UV and X-ray spectra of blazars have been often modeled as synchrotron component arising from a broken power-law distribution of electrons. A broken power-law distribution is expected, since the high energy electrons undergo radiative losses effectively. The change in the energy index should then be approximately 1 and corresponds to a spectral index difference of 0.5. However, one of the long outstanding problems has been that the observed index change is significantly different. On the other hand, recent high quality observations of blazars suggest that their local spectra may not be a power-law, instead have a slight curvature and often represented by a log parabola model. Using XMM-Newton observations spanning over 12 years for the BL Lac PKS 2155-304, we show that the optical/UV and X-ray spectra can be well represented by a broken log parabola model. Further, we show that such a spectrum can indicate the energy dependence of the electron escape timescale from the main acceleration zone. This novel approach, besides addressing the observed difference in the photon spectral indices, also tries to explain the spectral turn over in far UV/soft X-rays as a result of the radiative losses.

The broadband spectral energy distribution(SED) of blazars is generally interpreted as radiation arising from synchrotron and inverse Compton mechanisms.Traditionally, the underlying source parameters responsible for these emission processes,like particle energy density, magnetic field, etc., are obtained through simple visual reproduction of the observed fluxes. However, this procedure is incapable of providing the confidence range on the estimated parameters. In this work, we propose an efficient algorithm to perform a statistical fit of the observed broadband spectrum of blazars using different emission models. Moreover, in this work we use the the observable quantities as the fit parameters, rather than the direct source parameters which govern the resultant SED. This significantly improves the convergence time and eliminates the uncertainty regarding the initial guess parameters. This approach also has an added advantage of identifying the degenerate parameters, which can be removed by including more observable information and/or additional constraints. A computer code developed based on this algorithm is implemented as an user-defined routine in the standard X-ray spectral fitting package, XSPEC. Further, we demonstrate the efficacy of the algorithm by fitting the well sampled SED of the blazar, 3C 279, during its gamma ray flare in 2014.

We present the results of a multi-frequency, time-averaged analysis of blazars included in the Candidate Gamma-ray Blazar Survey catalog. Our sample consists of 324 $\gamma$-ray detected ($\gamma$-ray loud) and 191 non $\gamma$-ray detected ($\gamma$-ray quiet) blazars, and we consider all the data up to 2016 April 1. We find that both the $\gamma$-ray loud and the $\gamma$-ray quiet blazar populations occupy similar regions in the WISE color-color diagram, and in the radio and X-ray bands $\gamma$-ray loud sources are brighter. A simple one-zone synchrotron inverse-Compton emission model is applied to derive the physical properties of both populations. We find that the central black hole mass and the accretion disk luminosity ($L_{\rm disk}$) computed from the modeling of the optical-UV emission with a Shakura-Sunyaev disk reasonably matches with that estimated from the optical spectroscopic emission-line information. A significantly larger Doppler boosting in the $\gamma$-ray loud blazars is noted, and their jets are more radiatively efficient. On the other hand, the $\gamma$-ray quiet objects are more MeV-peaked, thus could be potential targets for next-generation MeV missions. Our results confirm the earlier findings about the accretion-jet connection in blazars; however, many of the $\gamma$-ray quiet blazars tend to deviate from the recent claim that the jet power exceeds $L_{\rm disk}$ in blazars. A broadband study, considering a larger set of $\gamma$-ray quiet objects and also including BL Lacs, will be needed to confirm/reject this hypothesis and also to verify the evolution of the powerful high-redshift blazars into their low-power nearby counterparts.

The synchrotron peak of the X-ray bright High Energy Peaked Blazar (HBL) PKS 2155$-$304 occurs in the UV-EUV region and hence its X-ray emission (0.6--10 keV) lies mostly in the falling part of the synchrotron hump. We aim to study the X-ray emission of PKS 2155$-$304 during different intensity states in 2009$-$2014 using XMM$-$Newton satellite. We studied the spectral curvature of all of the observations to provide crucial information on the energy distribution of the non-thermal particles. Most of the observations show curvature or deviation from a single power-law and can be well modeled by a log parabola model. In some of the observations, we find spectral flattening after 6 keV. In order to find the possible origin of the X-ray excess, we built the Multi-band Spectral Energy distribution (SED). We find that the X-ray excess in PKS 2155--304 is difficult to fit in the one zone model but, could be easily reconciled in the spine/layer jet structure. The hard X-ray excess can be explained by the inverse Comptonization of the synchrotron photons (from the layer) by the spine electrons.

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.

The observations of the prompt emission of gamma ray bursts (GRB) by GLAST Burst Monitor (GBM), on board Fermi Gamma-ray Space Telescope, suggest the presence of a significant thermal spectral component, whose origin is not well understood. Recently, it has been shown that for long duration GRBs, the spectral width as defined as the logarithm of the ratio of the energies at which the spectrum falls to half its peak value, lie in the range of 0.84-1.3 with a median value of 1.07. Thus, while most of the GRB spectra are found to be too narrow to be explained by synchrotron emission from an electron distribution, they are also significantly broader than a blackbody spectrum whose width should be 0.54. Here, we consider the possibility that an intrinsic thermal spectrum from a fire-ball like model, may be observed to be broadened if the system undergoes a rapid temperature evolution. We construct a toy-model to show that for bursts with durations in the range 5-70 seconds, the widths of their 1 second time-averaged spectra can be at the most $\lesssim 0.557$. Thus, while rapid temperature variation can broaden the detected spectral shape, the observed median value of $\sim 1.07$ requires that there must be significant sub-photospheric emission and/or an anisotropic explosion to explain the broadening for most GRB spectra.

We have observed Markarian 501 in Very High Energy (VHE) gamma-ray wavelength band for 70.6 hours from 15 April to 30 May, 2012 using TACTIC telescope. Detailed analysis of $\sim$66.3 hours of clean data revealed the presence of a TeV $\gamma$-ray signal (686$\pm$77 $\gamma$-ray events) from the source direction with a statistical significance of 8.89$\sigma$ above 850 GeV. Further, a total of 375 $\pm$ 47 $\gamma$-ray like events were detected in 25.2 hours of observation from 22 - 27 May, 2012 with a statistical significance of 8.05$\sigma$ indicating that the source has possibly switched over to a relatively high gamma-ray emission state. We have derived time-averaged differential energy spectrum of the state in the energy range 850 GeV - 17.24 TeV which fits well with a power law function of the form $dF/dE=f_{0}E^{-\Gamma}$ with $f_{0}= (2.27 \pm 0.38) \times 10 ^{-11} $ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$ and $\Gamma=2.57 \pm 0.15$. In order to investigate the source state, we have also used almost simultaneous multiwavelength observations viz: high energy data collected by $\it{Fermi}$-LAT, X-ray data collected by $\it{Swift}$-XRT and MAXI, optical and UV data collected by $\it{Swift}$-UVOT, and radio data collected by OVRO, and reconstructed broad-band Spectral Energy Distribution (SED). The obtained SED supports leptonic model (homogeneous single zone) for VHE gamma-ray emission involving synchrotron and synchrotron self Compton (SSC) processes.

Strong X-ray and $\gamma$--ray flares have been detected in February 2010 from the high synchrotron peaked blazar Mrk 421 (z=0.031). With the motivation of understanding the physics involved in this flaring activity, we study the variability of the source in X-ray and $\gamma$--ray energy bands during the period February 10-23, 2010 (MJD 55237-55250). We use near simultaneous X-ray data collected by \emph{MAXI}, \emph{Swift}-XRT and $\gamma$--ray data collected by \emph{Fermi}-LAT and \emph{TACTIC} along with the optical V-band observations by \emph{SPOL} at Steward Observatory. We observe that the variation in the one day averaged flux from the source during the flare is characterized by fast rise and slow decay. Besides, the TeV $\gamma$--ray flux shows a strong correlation with the X-ray flux, suggesting the former to be an outcome of synchrotron self Compton emission process. To model the observed X-ray and $\gamma$--ray light curves, we numerically solve the kinetic equation describing the evolution of particle distribution in the emission region. The injection of particle distribution into the emission region, from the putative acceleration region, is assumed to be a time dependent power law. The synchrotron and synchrotron self Compton emission from the evolving particle distribution in the emission region are used to reproduce the X-ray and $\gamma$--ray flares successfully. Our study suggests that the flaring activity of Mrk 421 can be an outcome of an efficient acceleration process associated with the increase in underlying non-thermal particle distribution.

A detailed multi-epoch study of the broadband spectral behaviour of the very high energy (VHE) source, 1ES\,1011+496, provides us with valuable information regarding the underlying particle distribution. Simultaneous observations of the source at optical/ UV/ X-ray/ $\gamma$-ray during three different epochs, as obtained from Swift-UVOT/ Swift-XRT/ Fermi-LAT, are supplemented with the information available from the VHE telescope array, HAGAR. The longterm flux variability at the Fermi-LAT energies is clearly found to be lognormal. It is seen that the broadband spectral energy distribution (SED) of 1ES\,1011+496 can be successfully reproduced by synchrotron and synchrotron self Compton emission models. Notably, the observed curvature in the photon spectrum at X-ray energies demands a smooth transition of the underlying particle distribution from a simple power law to a power law with an exponential cutoff or a smooth broken power law distribution, which may possibly arise when the escape of the particles from the main emission region is energy dependent. Specifically, if the particle escape rate is related to its energy as $E^{0.5}$ then the observed photon spectrum is consistent with the ones observed during the various epochs.

We present a systematic characterization of multi-wavelength emission from blazar PKS 1510-089 using well-sampled data at infrared(IR)-optical, X-ray and $\gamma$-ray energies. The resulting flux distributions, except at X-rays, show two distinct lognormal profiles corresponding to a high and a low flux level. The dispersions exhibit energy dependent behavior except for the LAT $\gamma$-ray and optical B-band. During the low level flux states, it is higher towards the peak of the spectral energy distribution, with $\gamma$-ray being intrinsically more variable followed by IR and then optical, consistent with mainly being a result of varying bulk Lorentz factor. On the other hand, the dispersions during the high state are similar in all bands expect optical B-band, where thermal emission still dominates. The centers of distributions are a factor of $\sim 4$ apart, consistent with anticipation from studies of extragalactic $\gamma$-ray background with the high state showing a relatively harder mean spectral index compared to the low state.

The well studied blazar Markarian 421 (Mrk 421, $z$=0.031) was the subject of an intensive multi-wavelength campaign when it flared in 2013 April. The recorded X-ray and very high energy (VHE, E$>$100 GeV) $\gamma$-ray fluxes are the highest ever measured from this object. At the peak of the activity, it was monitored by the hard X-ray focusing telescope {\it Nuclear Spectroscopic Telescope Array} ({\it NuSTAR}) and {\it Swift} X-Ray Telescope (XRT). In this work, we present a detailed variability analysis of {\it NuSTAR} and {\it Swift}-XRT observations of Mrk 421 during this flaring episode. We obtained the shortest flux doubling time of 14.01$\pm$5.03 minutes, which is the shortest hard X-ray (3$-$79 keV) variability ever recorded from Mrk 421 and is on the order of the light crossing time of the black hole's event horizon. A pattern of extremely fast variability events superposed on slowly varying flares is found in most of the {\it NuSTAR} observations. We suggest that these peculiar variability patterns may be explained by magnetic energy dissipation and reconnection in a fast moving compact emission region within the jet. Based on the fast variability, we derive a lower limit on the magnetic field strength of $B \ge 0.73 \delta_1^{-2/3} \, \nu_{19}^{1/3}$~G, where $\delta_1$ is the Doppler factor in units of 10, and $\nu_{19}$ is the characteristic X-ray synchrotron frequency in units of $10^{19}$~Hz.

CGRaBS J0809+5341, a high redshift blazar at z = 2.144, underwent a giant optical outburst on 2014 April 19 when it brightened by $\sim$5 mag and reached an unfiltered apparent magnitude of 15.7 mag. This implies an absolute magnitude of -30.5 mag, making it one of the brightest quasars in the Universe. This optical flaring triggered us to carry out observations during the decaying part of the flare covering a wide energy range using the {\it Nuclear Spectroscopic Telescope Array}, {\it Swift}, and ground based optical facilities. For the first time, the source is detected in $\gamma$-rays by the Large Area Telescope onboard the {\it Fermi Gamma-Ray Space Telescope}. A high optical polarization of $\sim$10\% is also observed. Using the Sloan Digital Sky Survey spectrum, accretion disk luminosity and black hole mass are estimated as $1.5 \times 10^{45}$ erg s$^{-1}$ and $10^{8.4}~M_{\odot}$ respectively. Using a single zone leptonic emission model, we reproduce the spectral energy distribution of the source during the flaring activity. This analysis suggests that the emission region is probably located outside the broad line region, and the jet becomes radiatively efficient. We also show that the overall properties of CGRaBS J0809+5341 seems not to be in agreement with the general properties observed in high redshift blazars up to now.

The well studied blazar 3C 279 underwent a giant $\gamma$-ray outburst in 2014 March-April. The measured $\gamma$-ray flux (1.21 $\pm$ 0.10 $\times$ 10$^{-5}$ ph cm$^{-2}$ s$^{-1}$ in 0.1-300 GeV energy range) is the highest detected from 3C 279 by Fermi Large Area Telescope. Hour scale $\gamma$-ray flux variability are observed, with a flux doubling time as short as 1.19 $\pm$ 0.36 hours detected during one flare. The $\gamma$-ray spectrum is found to be curved at peak of the flare suggesting low probability of detecting very high energy (VHE; E $>$ 100 GeV) emission, which is further confirmed by the Very Energetic Radiation Imaging Telescope Array System observations. The $\gamma$-ray flux increased by more than an order in comparison to low activity state and the flare consists of multiple sub-structures having fast rise and slow decay profile. The flux enhancement is seen in all the wavebands though at a lesser extent compared to $\gamma$-rays. During the flare, a considerable amount of the kinetic jet power gets converted to $\gamma$-rays and the jet becomes radiatively efficient. A one zone leptonic emission model is used to reproduce the flare and we find increase in the bulk Lorentz factor as a major cause of the outburst. From the observed fast variability, lack of VHE detection, and the curved $\gamma$-ray spectrum, we conclude that the location of the emission region cannot be far out from the broad line region (BLR) and contributions from both BLR and torus photons are required to explain the observed $\gamma$-ray spectrum.

We present a high time resolution study of the two brightest $\gamma$-ray outbursts from a blazar PKS 1222+216 observed by the \textit{Fermi} Large Area Telescope (LAT) in 2010. The $\gamma$-ray light-curves obtained in four different energy bands: 0.1--3, 0.1--0.3, 0.3--1 and 1--3 GeV, with time bin of 6 hr, show asymmetric profiles with a similar rise time in all the bands but a rapid decline during the April flare and a gradual one during the June. The light-curves during the April flare show $\sim 2$ days long plateau in 0.1--0.3 GeV emission, erratic variations in 0.3--1 GeV emission, and a daily recurring feature in 1--3 GeV emission until the rapid rise and decline within a day. The June flare shows a monotonic rise until the peak, followed by a gradual decline powered mainly by the multi-peak 0.1--0.3 GeV emission. The peak fluxes during both the flares are similar except in the 1--3 GeV band in April which is twice the corresponding flux during the June flare. Hardness ratios during the April flare indicate spectral hardening in the rising phase followed by softening during the decay. We attribute this behavior to the development of a shock associated with an increase in acceleration efficiency followed by its decay leading to spectral softening. The June flare suggests hardening during the rise followed by a complicated energy dependent behavior during the decay. Observed features during the June flare favor multiple emission regions while the overall flaring episode can be related to jet dynamics.

The very high energy (VHE) gamma ray spectral index of high energy peaked blazars correlates strongly with its corresponding redshift whereas no such correlation is observed in the X-ray or the GeV bands. We attribute this correlation to a result of photon-photon absorption of TeV photons with the extragalactic background light (EBL) and utilizing this, we compute the allowed flux range for the EBL, which is independent of previous estimates. The observed VHE spectrum of the sources in our sample can be well approximated by a power-law, and if the de-absorbed spectrum is also assumed to be a power law, then we show that the spectral shape of EBL will be $\epsilon n(\epsilon) \sim k log(\frac{\epsilon}{\epsilon_p}) $. We estimate the range of values for the parameters defining the EBL spectrum, $k$ and $\epsilon_p$, such that the correlation of the intrinsic VHE spectrum with redshift is nullified. The estimated EBL depends only on the observed correlation and the assumption of a power law source spectrum. Specifically, it does not depend on the spectral modeling or radiative mechanism of the sources, nor does it depend on any theoretical shape of the EBL spectrum obtained through cosmological calculations. The estimated EBL spectrum is consistent with the upper and lower limits imposed by different observations. Moreover, it also agrees closely with the theoretical estimates obtained through cosmological evolution models.

We present results from multi-wavelength study of intense flaring activity from a high frequency peaked BL Lac object Mrk 421. The source was observed in its flaring state on February 16, 2010 with the $TACTIC$ at energies above 1.5 TeV. Near simultaneous multi-wavelength data were obtained from high energy (MeV-GeV) $\gamma$--ray observations with \emph{Fermi}--LAT, X--ray observations by the \emph{Swift} and \emph{MAXI} satellites, optical V-band observation by SPOL at \emph{Steward Observatory} and radio 15 GHz observation at OVRO 40 meter-telescope. We have performed a detailed spectral and temporal analysis of $TACTIC$, \emph{Fermi}--LAT and \emph{Swift}--XRT observations of Mrk 421 during February 10--23, 2010 (MJD 55237-55250). The flaring activity of the source is studied by investigating the properties of daily light curves from radio to $TeV$ energy range and we present the correlation and variability analysis in each energy band. The $TeV$ flare detected by $TACTIC$ on February 16, 2010 is well correlated with the activity in lower energy bands. The differential energy spectrum of the source, in the energy range 1.5-11 TeV, as measured by $TACTIC$ on this night is described by a power law ($dN/dE \propto E^{-\Gamma}$) with spectral index $\Gamma = 2.6\pm0.3$. After accounting for absorption of $TeV$ photons by low energy extragalactic background light photons via pair production, the intrinsic $TeV$ spectrum reveals a power law index of $2.3\pm0.3$. Finally the broad band spectral energy distribution of the source in flaring state is reproduced using a simple emission model involving synchrotron and synchrotron self Compton processes. The obtained parameters are then used to understand the energetics of the source during the flaring episode.

We present a multi-wavelength study of the radio-loud narrow line Seyfert 1 galaxy (NLSy1), 1H 0323+342, detected by Fermi Gamma Ray Space Telescope. Multi-band light curves show many orphan X-ray and optical flares having no corresponding {\gamma}-ray counterparts. Such anomalous variability behavior can be due to different locations of the emission region from the central source. During a large flare, {\gamma}-ray flux doubling time scale as small as $\sim$ 3 hours is noticed. We built spectral energy distribution (SED) during different activity states and modeled them using an one-zone leptonic model. The shape of the optical/UV component of the SEDs is dominated by accretion disk emission in all the activity states. In the X-ray band, significant thermal emission from the hot corona is inferred during quiescent and first flaring states, however, during subsequent flares, non-thermal jet component dominates. The {\gamma}-ray emission in all the states can be well explained by inverse-Compton scattering of accretion disk photons reprocessed by the broad line region. The source showed violent intra-night optical variability, coinciding with one of the high {\gamma}-ray activity states. An analysis of the overall X-ray spectrum fitted with an absorbed power-law plus relativistic reflection component hints for the presence of Fe K-{\alpha} line and returns a high black hole spin value of a=0.96 $\pm$ 0.14. We argue that 1H 0323+342 possesses dual characteristics, akin to flat spectrum radio quasars (FSRQs) as well as radio-quiet NLSy1s, though at a low jet power regime compared to powerful FSRQs.

The $\gamma$-ray flare of PKS 1222+216, observed in June 2010, is interpreted as an outcome of jet dynamics at recollimation zone. We obtained the $\gamma$-ray light-curves in three different energy bands, namely, 100--300 MeV, 300 MeV--1 GeV and 1--3 GeV from observations by the \emph{Fermi} Large Area Telescope (LAT). We also use the \emph{Swift}--XRT flux from 0.3--10 keV obtained from archival data. We supplement these with the 0.07--0.4 TeV observations with MAGIC telescope, available in the literature. The detection of source at very high energy (VHE, $E>100$ GeV) with a differential photon spectral index of $2.7\pm0.3$ and the rapid variability associated with it suggests that the emission arises from a compact region located beyond the broad line emitting region. The plausible $\gamma$-ray emission mechanism can then be inverse Compton scattering of IR photons from obscuring torus. Further, the decay time of LAT flare cannot be explained by considering simple radiative loss mechanisms. Hence, to interpret the LAT light curves, we develop a model where the broadband emission originates from a compact region, arising plausibly from the compression of jet matter at the recollimation zone. The flare is then expressed as an outcome of jet deceleration probably associated with this focusing effect. The parameters of the model are further constrained by reproducing the broadband spectral energy distribution of the source obtained during the flare episode. Our study suggests that the particle energy density exceeds magnetic energy density by a large factor which in turn may cause rapid expansion of the emission region. However, near equipartition can be achieved towards the end of LAT flare during which the compact emission region would have expanded to the size of jet cross-section.

Very high energy (VHE) $\gamma$-rays above 10$'$s of GeV energy, emitted from distant blazars, are attenuated by photons from the extragalactic background light (EBL). Unfortunately, neither the EBL nor the intrinsic blazar spectrum is accurately known to derive one quantity from the other. In this work we use a homogeneous one zone model involving synchrotron, synchrotron self Compton (SSC) and external Compton (EC) emission mechanisms to estimate the intrinsic VHE spectra of blazars. The model is applied on three VHE blazars, namely PKS2155-304, RGB J0710+591 and 3C279, for which simultaneous multi-wavelength data are available from various observations. The predicted values of the intrinsic VHE fluxes are then compared with the observations by imaging atmospheric Cherenkov telescopes to determine the optical depth of VHE $\gamma$-rays. On comparing these optical depth values with those predicted by four different EBL models, we observe a somewhat pronounced systematic deviation for PKS2155-304 and 3C279 at higher energies, especially for the EBL model proposed by Finke et al.(2010). We attribute this deviation to be an outcome of either the failure of the extrapolation of blazar SED to VHE energies and/or due to various assumptions buried in the EBL models.

The broadband spectrum of a BL Lac object, OJ 287, from radio to $\gamma$-rays obtained during a major $\gamma$-ray flare detected by \emph{Fermi} in 2009 are studied to understand the high energy emission mechanism during this episode. Using a simple one-zone leptonic model, incorporating synchrotron and inverse Compton emission processes, we show that the explanation of high energy emission from X-rays to $\gamma$-rays, by considering a single emission mechanism, namely, synchrotron self-Compton (SSC) or external Compton (EC) requires unlikely physical conditions. However, a combination of both SSC and EC mechanisms can reproduce the observed high energy spectrum satisfactorily. Using these emission mechanisms we extract the physical parameters governing the source and its environment. Our study suggests that the emission region of OJ 287 is surrounded by a warm infrared (IR) emitting region of $\sim 250 \, K$. Assuming this region as a spherical cloud illuminated by an accretion disk, we obtain the location of the emission region to be $\sim 9 pc$. This supports the claim that the $\gamma$-ray emission from OJ 287 during the 2009 flare arises from a location far away from the central engine as deduced from millimeter-gamma ray correlation study and very long baseline array images.

Variable gamma-ray emission has been discovered in five Radio-loud Narrow Line Seyfert 1 (NLSy1) galaxies by the Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope. This has clearly demonstrated that these NLSy1 galaxies do have relativistic jets similar to two other cases of gamma-ray emitting Active Galactic Nuclei (AGN), namely blazars and radio galaxies. We present here our results on the multi-band analysis of two gamma-ray emitting NLSy1 galaxies namely PKS 1502+036 (z = 0.409) and PKS 2004-447 (z = 0.240) using archival data. We generate multi-band long term light curves of these sources, build their spectral energy distribution (SED) and model them using an one zone leptonic model. They resemble more to the SEDs of the flat spectrum radio quasar (FSRQ) class of AGN. We then compare the SEDs of these two sources with two other Fermi detected AGN along the traditional blazar sequence, namely the BL Lac Mrk 421 (z = 0.03) and the FSRQ 3C 454.3 (z = 0.86). The SEDs of both PKS 1502+036 and PKS 2004-447 are found to be intermediate to the SEDs of Mrk 421 and 3C 454.3. In the gamma-ray spectral index v/s gamma-ray luminosity plane, both these NLSy1 galaxies occupy a distinct position, wherein, they have luminosity between Mrk 421 and 3C 454.3, however steep gamma-ray spectra similar to 3C 454.3. Their Compton dominance as well as their X-ray spectral slope also lie between Mrk 421 and 3C 454.3. We argue that the physical properties of both PKS 1502+036 and PKS 2004$-$447 are in general similar to blazars and intermediate between FSRQs and BL Lac objects and these sources thus could fit into the traditional blazar sequence.

We have observed the blazar Markarian 421 with the TACTIC $\gamma$-ray telescope at Mt. Abu, India, from 22 November 2009 to 16 May 2010 for 265 hours. Detailed analysis of the data so recorded revealed presence of a TeV $\gamma$-ray signal with a statistical significance of 12.12$\sigma$ at $E_{\gamma}\geq$ 1 TeV. We have estimated the time averaged differential energy spectrum of the source in the energy range 1.0 - 16.44 TeV. The spectrum fits well with the power law function of the form ($dF/dE=f_0 E^{-\Gamma}$) with $f_0=(1.39\pm0.239)\times 10^{-11}cm^{-2}s^{-1}TeV^{-1}$ and $\Gamma=2.31\pm0.14$.

We model the simultaneous observations of the flat spectrum radio quasar 3C 279 at radio, optical, X-ray and very high energy (VHE) gamma ray energies during 2006 flare using a simple one zone leptonic model. We consider synchrotron emission due to cooling of a non-thermal electron distribution in an equipartition magnetic field and inverse Compton emission due to the scattering off synchrotron photons (SSC) and external soft photons (EC) by the same distribution of electrons. We show that the VHE gamma ray flux cannot be explained by SSC process thereby suggesting the EC mechanism as a plausible emission process at this energy. The EC scattering of BLR photons to VHE energies will be in Klein-Nishina regime predicting a steep spectrum which is contrary to the observations. However the infrared photons from the dusty torus can be boosted to VHE energies with the scattering process remaining in the Thomson regime. Though EC process can successfully explain the observed VHE flux, it require a magnetic field much lower than the equipartition value to reproduce the observed X-ray flux. Hence we attribute the X-ray emission due to SSC process. We derive the physical parameters of 3C 279 considering the above mentioned emission processes. In addition we assume the size of emission region constrained by a variability timescale of one day. This model can successfully reproduce the broadband spectrum of 3C 279 but predicts substantially large flux at MeV-GeV energies.

The high resolution image of the jet of the BL Lac object MKN501 in radio, show a limb-brightened feature. An explanation of this feature as an outcome of differential Doppler boosting of jet spine and jet boundary due to transverse velocity structure of the jet requires large viewing angle. However this inference contradicts with the constraints derived from the high energy $\gamma$-ray studies unless the jets bends over a large angle immediately after the $\gamma$-ray zone (close to the central engine). In this letter we propose an alternate explanation to the limb-brightened feature of MKN501 by considering the diffusion of electrons accelerated at the boundary shear layer into the jet medium and this consideration does not require large viewing angle. Also the observed difference in the spectral index at the jet boundary and jet spine can be understood within the frame work of shear acceleration.

The flux and the spectral index in X-ray energy band from the knots of M87 jet as observed by {\it{Chandra}} indicate a possible synchrotron origin but cannot be explained by simple one zone models with continuous injection of non-thermal electrons. In this letter we propose a two-zone model to explain the observed spectra of the knots of M87 jet. We consider the synchrotron emission from a region with tangled magnetic field where relativistic non-thermal electrons are continuously injected in from an associated acceleration region. The acceleration region is assumed to be compact zone possibly around a shock front. A power-law distribution of electron is injected into the acceleration region and are accelerated to a maximum energy determined by the acceleration time scale and the loss processes. With the present model we are able to explain the overall broadband features of the knots of M87 jet. Also the present model predicts a change in spectral index at ultraviolet energies and future observations at these energies can be used to constrain the parameters involved in the model.

The TACTIC $\gamma$-ray telescope, equipped with a light collector of area $\sim$9.5m$^2$ and a medium resolution imaging camera of 349-pixels, has been in operation at Mt.Abu, India since 2001. This paper describes the main features of its various subsystems and its overall performance with regard to (a) tracking accuracy of its 2-axes drive system, (b) spot size of the light collector, (c) back-end signal processing electronics and topological trigger generation scheme, (d) data acquisition and control system and (e) relative and absolute gain calibration methodology. Using a trigger field of view of 11$\times$11 pixels ($\sim$ 3.4$^\circ$$\times3.4^\circ$), the telescope records a cosmic ray event rate of $\sim$2.5 Hz at a typical zenith angle of 15$^\circ$. Monte Carlo simulation results are also presented in the paper for comparing the expected performance of the telescope with actual observational results. The consistent detection of a steady signal from the Crab Nebula above $\sim$1.2 TeV energy, at a sensitivity level of $\sim5.0\sigma$ in $\sim$25 h, alongwith excellent matching of its energy spectrum with that obtained by other groups, reassures that the performance of the TACTIC telescope is quite stable and reliable. Furthermore, encouraged by the detection of strong $\gamma$-ray signals from Mrk 501 (during 1997 and 2006 observations) and Mrk 421 (during 2001 and 2005-2006 observations), we believe that there is considerable scope for the TACTIC telescope to monitor similar TeV $\gamma$-ray emission activity from other active galactic nuclei on a long term basis.

The TACTIC $\gamma$-ray telescope has observed Mrk 421 on 66 clear nights from Dec. 07, 2005 to Apr. 30, 2006, totalling $\sim$ 202 hours of on-source observations. Here, we report the detection of flaring activity from the source at $\geq$ 1 TeV energy and the time-averaged differential $\gamma$-ray spectrum in the energy range 1-11 TeV for the data taken between Dec. 27, 2005 to Feb. 07, 2006 when the source was in a relatively higher state as compared to the rest of the observation period. Analysis of this data spell, comprising about $\sim$97h reveals the presence of a $\sim 12.0 \sigma$ $\gamma$-ray signal with daily flux of $>$ 1 Crab unit on several days. A pure power law spectrum with exponent $-3.11\pm0.11$ as well as a power law spectrum with an exponential cutoff $(\Gamma = -2.51\pm0.26$ and $E_0=(4.7\pm2.1) TeV)$ are found to provide reasonable fits to the inferred differential spectrum within statistical uncertainties. We believe that the TeV light curve presented here, for nearly 5 months of extensive coverage, as well as the spectral information at $\gamma$-ray energies of $>$ 5 TeV provide a useful input for other groups working in the field of $\gamma$-ray astronomy.

The dynamics of relativistically moving blobs ejected out of a central AGN, are considered. It is assumed that the collision between two blobs are completely inelastic, such that the bulk kinetic energy lost in the collision is used to energize electrons to relativistic energies via acceleration in internal shocks which are formed due to the collision. These high energy electrons which are produced on a time-scale corresponding to the collision time-scale, cool by radiative losses due to synchrotron and Inverse Compton processes. The model is applied to the radio/X-ray knots of several AGN. For three of these sources we have analyzed long (> 40 ksec) Chandra observations and report on constrains on the X-ray spectral indices. In the framework of this model the AGN are inferred to sporadically eject relativistic blobs on time-scales ranging from 10^{11} to 10^{12} secs for different sources. It is shown that the collision time-scales can be longer than the age of the knot and, hence, non-thermal electrons are continuously being injected into the system. This continuous injection, in contrast to an instantaneous one time injection, gives rise to a characteristic spectral break, rather than a high-energy cutoff in the spectrum.

We consider the evolution of a spherically expanding plasma cloud, where there is continuous injection of non-thermal electrons. We compute the time dependent electron distribution and resultant photon spectra taking into account synchrotron, adiabatic and inverse Compton cooling. This model is different from previous works where, instead of a continuous injection of particles, a short injection period was assumed. We apply this model to the radio/optical knots in the large scale jets of AGN, detected in X-rays by {\it Chandra} and find that the overall broadband spectral features can be reproduced. It is shown that for some sources, constraints on the X-ray spectral index (by a longer {\it Chandra} observation) will be able to differentiate between the different models. This in turn will put a strong constraint on the acceleration mechanism active in these sources.

The results of our observational campaigns on the two extragalactic sources Mkn 501 and Mkn 421, carried out with the Imaging Element of the TACTIC array, during March-May, 1998 and April-May, 1999, are presented. The results indicate that the two BL Lac objects (Mkn 501 and Mkn 421) were in a `low' gamma-ray emission state during both epochs of our observations.