The uncharted territory of the high-redshift ($z\gtrsim3$) Universe holds the key to understand the evolution of quasars. In an attempt to identify the most extreme members of the quasar population, i.e., blazars, we have carried out a multi-wavelength study of a large sample of radio-loud quasars beyond $z=3$. Our sample consists of 9 $\gamma$-ray detected blazars and 133 candidate blazars selected based on the flatness of their soft X-ray spectra (0.3$-$10 keV photon index $\leq1.75$), including 15 with NuSTAR observations. The application of the likelihood profile stacking technique reveals that the high-redshift blazars are faint $\gamma$-ray emitters with steep spectra. The high-redshift blazars host massive black holes ($\langle \log~M_{\rm BH,~M{\odot}} \rangle>9$) and luminous accretion disks ($\langle L_{\rm disk} \rangle>10^{46}$ erg s$^{-1}$). Their broadband spectral energy distributions are found to be dominated by high-energy radiation indicating their jets to be among the most luminous ones. Focusing on the sources exhibiting resolved X-ray jets (as observed with the Chandra satellite), we find the bulk Lorentz factor to be larger with respect to other $z>3$ blazars, indicating faster moving jets. We conclude that the presented list of the high-redshift blazars may act as a reservoir for follow-up observations, e.g., with NuSTAR, to understand the evolution of relativistic jets at the dawn of the Universe.

What triggers collimated relativistic outflows or jets, from the centers of galaxies remains a fundamental question in astrophysics. The merging of two galaxies has been proposed to realize the conditions to successfully launch and drive such jets into the intergalactic medium. However, evidences for the operation of this mechanism are scarce. Here we report the first unambiguous detection of an ongoing merger of a narrow-line Seyfert 1 galaxy, TXS 2116$-$077, hosting a closely aligned, $\gamma$-ray emitting relativistic jet with a Seyfert 2 galaxy at a separation of $\sim$12 kpc, using the observations taken with 8.2 m Subaru telescope. Our subsequent followup observations with 10.4 m Gran Telescopio Canarias, 4.2 m William Herschel Telescope, and Chandra X-ray observatory have provided what is likely to be the first glimpse of the merging environment hosting a closely aligned relativistic jet. Our finding that the jet is considerably younger than the merger demonstrates that jet activity can be triggered by galaxy mergers and that $\gamma$-ray detected narrow-line Seyfert 1 galaxies represent the beginning phase of that activity. These results also highlight the crucial role of mergers in shaping the fate of galaxies in their cosmological evolution and are consistent with recent studies focused on the host galaxy imaging of this enigmatic class of active galactic nuclei.

A majority of the $\gamma$-ray emission from star-forming galaxies is generated by the interaction of high-energy cosmic rays with the interstellar gas and radiation fields. Star-forming galaxies are expected to contribute to both the extragalactic $\gamma$-ray background and the IceCube astrophysical neutrino flux. Using roughly 10\,years of $\gamma$-ray data taken by the {\it Fermi} Large Area Telescope, in this study we constrain the $\gamma$-ray properties of star-forming galaxies. We report the detection of 11 bona-fide $\gamma$-ray emitting galaxies and 2 candidates. Moreover, we show that the cumulative $\gamma$-ray emission of below-threshold galaxies is also significantly detected at $\sim$5\,$\sigma$ confidence. The $\gamma$-ray luminosity of resolved and unresolved galaxies is found to correlate with the total (8-1000\,$\mu$m) infrared luminosity as previously determined. Above 1\,GeV, the spectral energy distribution of resolved and unresolved galaxies is found to be compatible with a power law with a photon index of $\approx2.2-2.3$. Finally, we find that star-forming galaxies account for roughly 5\,\% and 3\,\% of the extragalactic $\gamma$-ray background and the IceCube neutrino flux, respectively.

Motivated by the identification of the blazar TXS 0506+056 as the first promising high-energy neutrino counterpart candidate, we search for additional neutrino blazars candidates among the Fermi-LAT detected blazars. We investigate the multi-wavelength behavior from radio to GeV gamma rays of blazars found to be in spatial coincidence with single high-energy neutrinos and lower-energy neutrino flare candidates. In addition, we compare the average gamma-ray emission of the potential neutrino-emitting sources to the entire sample of gamma-ray blazars. We find that neutrino-emitting blazar candidates are statistically compatible with both hypothesis of a linear correlation and of no correlation between neutrino and gamma-ray energy flux.

The recent spatial and temporal coincidence of the blazar TXS 0506+056 with the IceCube detected neutrino event IC-170922A, has opened up a realm of multi-messenger astronomy with blazar jets as a plausible site of cosmic-ray acceleration. After TXS 0505+056, for the first time, a second blazar, BZB J0955+3551, has recently been found to be spatially coincident with the IceCube detected neutrino event IC-200107A and undergoing its brightest X-ray flare measured so far. Here, we present the results of our multi-frequency campaign to study this peculiar event that includes observations from NuSTAR, Swift, NICER, and 10.4 m Gran Telescopio Canaris (GTC). The optical spectroscopic observation from GTC secured its redshift as $z=0.5573$. Both NuSTAR and NICER data reveal a rapid flux variability albeit at low-significance ($\lesssim3.5\sigma$). We explore the origin of the target photon field needed for the photo-pion production using analytical calculations and considering the observed optical-to-X-ray flux level and conclude that seed photons must come from outside the jet, similar to that reported for TXS 0506+056. We find that there could be a sub-dominant hadronic population responsible for the neutrino emission.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

High synchrotron peak (HSP; $\nu_{sy}^{pk} > 10^{15}$ Hz) BL Lacs are some of the most extreme accelerators in the Universe. Those found at high redshifts ($z>1$) challenge our understanding of blazar evolution models and are crucial for cosmological measurements of the Extragalactic Background Light. In this paper, we study a high-$z$ BL Lac, 4FGL J2146.5-1344, detected to be at $z$=1.34 using the photometric dropout technique. We collected multi-wavelength data for this source from optical up to $\gamma$-rays, in order to study its spectral energy distribution (SED). In particular, this source was observed for the first time with NuSTAR, which accurately measures the synchrotron emission of this blazar up to 50 keV. Despite being classified as an HSP BL Lac object, the modeling of the SED reveals that this source likely belongs to the "masquerading BL Lac" class, which comprises of FSRQs appearing as disguised BL Lac objects.

With bolometric luminosities exceeding $10^{48}$ erg s$^{-1}$, powerful jets and supermassive black holes at their center, MeV blazars are some of the most extreme sources in the Universe. Recently, the Fermi-Large Area Telescope detected five new $\gamma$-ray emitting MeV blazars beyond redshift $z=3.1$. With the goal of precisely characterizing the jet properties of these extreme sources, we started a multiwavelength campaign to follow them up with joint NuSTAR, Swift and SARA observations. We observe six high-redshift quasars, four of them belonging to the new $\gamma$-ray emitting MeV blazars. Thorough X-ray analysis reveals spectral flattening at soft X-ray for three of these objects. The source NVSS J151002$+$570243 also shows a peculiar re-hardening of the X-ray spectrum at energies $E>6\,\rm keV$. Adopting a one-zone leptonic emission model, this combination of hard X-rays and $\gamma$-rays enables us to determine the location of the Inverse Compton peak and to accurately constrain the jet characteristics. In the context of the jet-accretion disk connection, we find that all six sources have jet powers exceeding accretion disk luminosity, seemingly validating this positive correlation even beyond $z>3$. Our six sources are found to have $10^9 \rm M_{\odot}$ black holes, further raising the space density of supermassive black holes in the redshift bin $z=[3,4]$.

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

The Hubble constant $H_{0}$ and matter density $\Omega_{m}$ of the Universe are measured using the latest $\gamma$-ray attenuation results from Fermi-LAT and Cherenkov telescopes. This methodology is based upon the fact that the extragalactic background light supplies opacity for very high energy photons via photon-photon interaction. The amount of $\gamma$-ray attenuation along the line of sight depends on the expansion rate and matter content of the Universe. This novel strategy results in a value of $H_{0}=67.4_{-6.2}^{+6.0}$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.14_{-0.07}^{+0.06}$. These estimates are independent and complementary to those based on the distance ladder, cosmic microwave background (CMB), clustering with weak lensing, and strong lensing data. We also produce a joint likelihood analysis of our results from $\gamma$ rays and these from more mature methodologies, excluding the CMB, yielding a combined value of $H_{0}=66.6\pm 1.6$~km~s$^{-1}$~Mpc$^{-1}$ and $\Omega_{m}=0.29\pm 0.02$.

This article reviews our current understanding about $\gamma$-ray detected narrow-line Seyfert 1 ($\gamma$-NLSy1) galaxies. The detection with the Large Area Telescope onboard {\it Fermi}~Gamma-ray Space Telescope has provided the strongest evidence for the presence of closely aligned relativistic jet in these intriguing active galactic nuclei (AGN) and opened up a realm to explore the physical conditions needed to launch the jet in a different central engine and host galaxy environment than that is known for blazars. Promising results acquired from various multi-wavelength campaigns are converging to a scenario in which the $\gamma$-NLSy1 galaxies can be considered as `young' blazars. These enigmatic sources hold the key to unravel the jet triggering mechanism and evolution of the AGN phase of a galaxy, in general. As such, $\gamma$-NLSy1s should be considered as one of the top priority targets for next generation observational facilities.

We present a likelihood profile stacking technique based on the {\it Fermi}-Large Area Telescope (LAT) data to explore the $\gamma$-ray characteristics of {\it Fermi}-LAT undetected astrophysical populations. The pipeline is applied to a sample of $\gamma$-ray unresolved extreme blazars, i.e., sources with the highest synchrotron peak frequencies ($\nu_{\rm Syn}^{\rm peak}\geqslant 10^{17}$ Hz), and we report a cumulative $\gamma$-ray detection with more than 32$\sigma$ confidence for 2 degrees of freedom. Comparing the generated stacked $\gamma$-ray spectrum with the sensitivity limits of the upcoming Cherenkov Telescope Array (CTA), we find that the {\it Fermi}-LAT undetected population of such extreme blazars, on average, may remain well below the CTA detection threshold due to their faintness and extragalactic background light (EBL) absorption. However, $\gamma$-ray detected blazars belonging to the same class are promising candidates for CTA observations. The EBL corrected stacked spectra of these sources do not show any softening up to 1 TeV. This finding suggests the inverse Compton peak of extreme blazars to lie above 1 TeV, thus indicating a hard intrinsic TeV spectrum. Our analysis also predicts that at 100 GeV, at least $\sim$10\% of the diffuse extragalactic $\gamma$-ray background originates from the $\gamma$-ray undetected extreme blazars. These results highlight the effectiveness of the developed stacking technique to explore the uncharted territory of $\gamma$-ray undetected astrophysical objects.

The recently released 105-month {\it Swift}-Burst Alert Telescope (BAT) all-sky hard X-ray survey catalog presents an opportunity to study astrophysical objects detected in the deepest look at the entire hard X-ray (14$-$195 keV) sky. Here we report the results of a multifrequency study of 146 blazars from this catalog, quadrupling the number compared to past studies, by utilizing recent data from the {\it Fermi}-Large Area Telescope (LAT), Swift-BAT, and archival measurements. In our $\gamma$-ray analysis of $\sim$10 years of the LAT data, 101 are found as $\gamma$-ray emitters, whereas, 45 remains LAT undetected. We model the broadband spectral energy distributions with a synchrotron-inverse Compton radiative model. On average, BAT detected sources host massive black holes ($M_{\rm bh}\sim10^9$ M$_{\odot}$) and luminous accretion disks ($L_{\rm d}\sim10^{46}$ erg s$^{-1}$). At high-redshifts ($z>2$), BAT blazars host more powerful jets with luminous accretion disks compared to those detected only with the {\it Fermi}-LAT. We find good agreement in the black hole masses derived from the single-epoch optical spectroscopic measurements and standard accretion disk modeling approaches. Other physical properties of BAT blazars are similar to those known for {\it Fermi}-LAT detected objects.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission (unresolved gamma-ray background, UGRB) below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This work presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with $\sim$3.7$\sigma$ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55 $\pm$ 0.23 and 1.86 $\pm$ 0.15.

MeV blazars are the most luminous persistent sources in the Universe and emit most of their energy in the MeV band. These objects display very large jet powers and accretion luminosities and are known to host black holes with a mass often exceeding $10^9 M_{\odot}$. An MeV survey, performed by a new generation MeV telescope which will bridge the entire energy and sensitivity gap between the current generation of hard X-ray and gamma-ray instruments, will detect $>$1000 MeV blazars up to a redshift of $z=5-6$. Here we show that this would allow us: 1) to probe the formation and growth mechanisms of supermassive black holes at high redshifts, 2) to pinpoint the location of the emission region in powerful blazars, 3) to determine how accretion and black hole spin interplay to power the jet.

The Extragalactic Background Light (EBL) can be probed via the absorption imprint it leaves in the spectra of gamma-ray sources ($\gamma\gamma \rightarrow e^-e^+$). We recently developed a dedicated technique to reconstruct the EBL, and its evolution with redshift, from $\gamma$ ray optical depth data using a large sample of blazars detected by the $Fermi$ Large Area Telescope. Here, we extend this dataset to the TeV regime using ground-based Cherenkov observations of 38 blazars and report the first homogeneous measurement of the EBL spectral intensity covering the ultraviolet to infrared wavelengths ($\sim$0.1-100$\mathrm{\mu m}$). A minimal EBL throughout the wavelength range with respect to integrated galaxy light is found, allowing little additional unresolved emission from faint or truly diffuse populations setting an upper limit of $\lesssim 4~{\rm nW\cdot m^{-2}sr^{-1}}$ at 1.4 ${\rm \mu m}$. In particular, the cosmic optical background (COB) at $z=0$ is found to be $27.8_{-2.0}^{+2.1}~{\rm nW\cdot m^{-2}sr^{-1}}$. This work lays the foundation for accurate gamma-ray measurements of the EBL across its whole spectral range using a combination of GeV and TeV data.

We present the results of a detailed multi-wavelength study of a sample of sixteen narrow-line Seyfert 1 (NLSy1) galaxies known in $\gamma$-rays so far. This includes a focused X-ray spectral analysis, to compare them with a more general NLSy1 population, and a broadband spectral energy distribution (SED) modeling considering the leptonic radiative processes to understand their similarity/dissimilarity with blazars. The X-ray spectra of $\gamma$-NLSy1 galaxies exhibit similarities both with typical NLSy1 sources below $\lesssim$2 keV and blazars at higher energies. We also find weak evidences of excess absorption in the soft X-ray spectra of 3 sources and a tentative detection of Fe line in 1H~0323+342. A comparison of the broadband parameters derived from the SED modeling reveals the similarity of $\gamma$-NLSy1 objects with blazars, in particular more with flat spectrum radio quasars. The $\gamma$-NLS1 galaxies host relatively low power jets with small bulk Lorentz factors with respect to blazars and this explains the paucity of such sources in the $\gamma$-ray catalogs and identification of an even fewer parent population members. Many of the observed features of these intriguing objects suggest they host low-mass black holes. The accretion rate in Eddington units can be considered as the driving factor to explain the physical properties of $\gamma$-NLSy1 sources, similar to that which is known for general blazar population.

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

We report the first results obtained from our campaign to characterize the intranight-optical variability (INOV) properties of {\it Fermi}~detected blazars, using the observations from the recently commissioned 1.3 m J C Bhattacharya telescope (JCBT). During the first run, we were able to observe 17 blazars in the Bessel $R$ filter for $\sim$137 hrs. Using $C$ and scaled $F$-statistics, we quantify the extent of INOV and derive the duty cycle (DC) which is the fraction of time during which a source exhibits a substantial flux variability. We find a high DC of 40\% for BL Lac objects and the flat spectrum radio quasars are relatively less variable (DC $\sim15$\%). However, when estimated for blazars sub-classes, a high DC of $\sim$59\% is found in low synchrotron peaked (LSP) blazars, whereas, intermediate and high synchrotron peaked objects have a low DC of $\sim$11\% and 13\%, respectively. We find evidences about the association of the high amplitude INOV with the $\gamma$-ray flaring state. We also notice a high polarization during the elevated INOV states (for the sources that have polarimetric data available), thus supporting the jet based origin of the observed variability. We plan to enlarge the sample and utilize the time availability from the small telescopes, such as 1.3 m JCBT, to strengthen/verify the results obtained in this work and those existed in the literature.

High-redshift ($z>2$) blazars are the most powerful members of the blazar family. Yet, only a handful of them have both X-ray and $\gamma$-ray detection, thereby making it difficult to characterize the energetics of the most luminous jets. Here, we report, for the first time, the Fermi-Large Area Telescope detection of the significant $\gamma$-ray emission from the high-redshift blazar DA 193 ($z=2.363$). Its time-averaged $\gamma$-ray spectrum is soft ($\gamma$-ray photon index = $2.9\pm0.1$) and together with a relatively flat hard X-ray spectrum (14$-$195 keV photon index = $1.5\pm0.4$), DA 193 presents a case to study a typical high-redshift blazar with inverse Compton peak being located at MeV energies. An intense GeV flare was observed from this object in the first week of 2018 January, a phenomenon rarely observed from high-redshift sources. What makes this event a rare one is the observation of an extremely hard $\gamma$-ray spectrum (photon index = $1.7\pm0.2$), which is somewhat unexpected since high-redshift blazars typically exhibit a steep falling spectrum at GeV energies. The results of our multi-frequency campaign, including both space- (Fermi, NuSTAR, and Swift) and ground-based (Steward and Nordic Optical Telescope) observatories, are presented and this peculiar $\gamma$-ray flare is studied within the framework of a single-zone leptonic emission scenario.

The light emitted by all galaxies over the history of the Universe produces the extragalactic background light (EBL) at ultraviolet, optical, and infrared wavelengths. The EBL is a source of opacity for $\gamma$ rays via photon-photon interactions, leaving an imprint in the spectra of distant $\gamma$-ray sources. We measure this attenuation using {739} active galaxies and one gamma-ray burst detected by the {\it Fermi} Large Area Telescope. This allows us to reconstruct the evolution of the EBL and determine the star-formation history of the Universe over 90\% of cosmic time. Our star-formation history is consistent with independent measurements from galaxy surveys, peaking at redshift $z\sim2$. Upper limits of the EBL at the epoch of re-ionization suggest a turnover in the abundance of faint galaxies at $z\sim 6$.

Here we explore the disk-jet connection in the broad-line radio quasar 4C+74.26, utilizing the results of the multiwavelength monitoring of the source. The target is unique in that its radiative output at radio wavelengths is dominated by a moderately-beamed nuclear jet, at optical frequencies by the accretion disk, and in the hard X-ray range by the disk corona. Our analysis reveals a correlation (local and global significance of 96\% and 98\%, respectively) between the optical and radio bands, with the disk lagging behind the jet by $250 \pm 42$ days. We discuss the possible explanation for this, speculating that the observed disk and the jet flux changes are generated by magnetic fluctuations originating within the innermost parts of a truncated disk, and that the lag is related to a delayed radiative response of the disk when compared with the propagation timescale of magnetic perturbations along relativistic outflow. This scenario is supported by the re-analysis of the NuSTAR data, modelled in terms of a relativistic reflection from the disk illuminated by the coronal emission, which returns the inner disk radius $R_{\rm in}/R_{\rm ISCO} =35^{+40}_{-16}$. We discuss the global energetics in the system, arguing that while the accretion proceeds at the Eddington rate, with the accretion-related bolometric luminosity $L_{\rm bol} \sim 9 \times 10^{46}$ erg s$^{-1}$ $\sim 0.2 L_{\rm Edd}$, the jet total kinetic energy $L_\textrm{j} \sim 4 \times 10^{44}$ erg s$^{-1}$, inferred from the dynamical modelling of the giant radio lobes in the source, constitutes only a small fraction of the available accretion power.

The HAWC (High Altitude Water Cherenkov) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100~GeV) gamma-ray sources based on 507 days of observation. Among these, there are nineteen sources that are not associated with previously known TeV sources. We have studied fourteen of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1~TeV-30~TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected fourteen new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected GeV gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.

We present the results of a study of the time-averaged spectral energy distributions (SEDs) of eight flat spectrum radio quasars (FSRQs) present in the second catalog of high energy sources detected beyond 50 GeV by the {\it Fermi}-Large Area Telescope (2FHL). Both leptonic and hadronic scenarios are adopted to explain the multi-wavelength SEDs and we find them to be marginally consistent with the 2FHL spectra above 50 GeV. We derive the expected degree of X-ray and $\gamma$-ray polarizations both for the average and elevated activity states and note that: (i) a hadronic radiative model consistently predicts a higher degree of high energy polarization compared to leptonic ones, and (ii) the X-ray polarization degree is higher than the $\gamma$-ray polarization in the leptonic scenario, but similar to the $\gamma$-ray polarization if the observed radiation is powered by hadronic processes. From the leptonic modeling, the location of the $\gamma$-ray emitting region is found to be at the outer edge of the broad line region (BLR) and is consistent with the $\gamma\gamma$ opacity estimates for the $\gamma$-ray absorption by the BLR. We conclude that a majority of the FSRQs could be detected by the upcoming Cherenkov Telescope Array, though future high energy polarimeters will be able to detect them only during elevated activity states, which could provide supportive evidence for the hadronic origin of the X-ray and $\gamma$-ray emission.

BL Lacertae (BL Lac) objects are the prominent members of the third {\it Fermi} Large Area Telescope catalog of $\gamma$-ray sources. Half of the BL Lac population ($\sim$ 300) lack redshift measurements, which is due to the absence of lines in their optical spectrum, thereby making it difficult to utilize spectroscopic methods. Our photometric drop-out technique can be used to establish the redshift for a fraction of these sources. This work employed 6 filters mounted on the $Swift$-UVOT and 4 optical filters on two telescopes, the 0.65 m SARA-CTIO in Chile and 1.0 m SARA-ORM in the Canary Islands, Spain. A sample of 15 sources was extracted from the $Swift$ archival data for which 6 filter UVOT observations were conducted. By complementing the {\it Swift} observations with the SARA ones, we were able to discover two high redshift sources: 3FGL J1155.4-3417 and 3FGL J1156.7-2250 at $z=1.83^{+0.10}_{-0.13}$ and $z=1.73^{+0.11}_{-0.19}$ , respectively, resulting from the dropouts in the powerlaw template fits to these data. The discoveries add to the important (26 total) sample of high-redshift BL Lacs. While the sample of high-z BL Lacs is still rather small, these objects do not seem to fit well within known schemes of the blazar population and represent the best probes of the extragalactic background light.

The spatial extension of a gamma-ray source is an essential ingredient to determine its spectral properties as well as its potential multi-wavelength counterpart. The capability to spatially resolve gamma-ray sources is greatly improved by the newly delivered Fermi-Large Area Telescope (LAT) Pass 8 event-level analysis which provides a greater acceptance and an improved point spread function, two crucial factors for the detection of extended sources. Here, we present a complete search for extended sources located within 7 degrees from the Galactic plane, using 6 years of LAT data above 10 GeV. We find 46 extended sources and provide their morphological and spectral characteristics. This constitutes the first catalog of hard LAT extended sources, named the Fermi Galactic Extended Source Catalog, which allows a thorough study of the properties of the Galactic plane in the sub-TeV domain.

The detection of significant $\gamma$-ray emission from radio-loud narrow line Seyfert 1 (NLSy1s) galaxies enables us to study jets in environments different than those in blazars. However, due to the small number of known $\gamma$-ray emitting NLSy1 ($\gamma$-NLSy1) galaxies, a comprehensive study could not be performed. Here we report the first detection of significant $\gamma$-ray emission from four active galactic nuclei (AGN), recently classified as NLSy1 from their Sloan Digital Sky Survey (SDSS) optical spectrum. Three flat spectrum radio quasars (FSRQs) present in the third Large Area Telescope AGN catalog (3LAC) are also found as $\gamma$-NLSy1 galaxies. Comparing the $\gamma$-ray properties of these objects with 3LAC blazars reveals their spectral shapes to be similar to FSRQs, however, with low $\gamma$-ray luminosity ($\lesssim10^{46-47}$ erg s$^{-1}$). In the Wide-field Infrared Survey Explorer color-color diagram, these objects occupy a region mainly populated by FSRQs. Using the H$_{\beta}$ emission line parameters, we find that on average $\gamma$-NLSy1 have smaller black hole masses than FSRQs at similar redshifts. In the low-resolution SDSS image of one of the $\gamma$-NLSy1 source, we find the evidence of an extended structure. We conclude by noting that overall many observational properties of $\gamma$-NLSy1 sources are similar to FSRQs and therefore, these objects could be their low black hole mass counterparts, as predicted in the literature.

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 extragalactic background light (EBL), from ultraviolet to infrared wavelengths, is predominantly due to emission from stars, accreting black holes and reprocessed light due to Galactic dust. The EBL can be studied through the imprint it leaves, via $\gamma$-$\gamma$ absorption of high-energy photons, in the spectra of distant $\gamma$-ray sources. The EBL has been probed through the search for the attenuation it produces in the spectra of BL Lacertae (BL Lac) objects and individual $\gamma$-ray bursts (GRBs). GRBs have significant advantages over blazars for the study of the EBL especially at high redshifts. Here we analyze a combined sample of twenty-two GRBs, detected by the Fermi Large Area Telescope between 65 MeV and 500 GeV. We report a marginal detection (at the ~2.8 $\sigma$ level) of the EBL attenuation in the stacked spectra of the source sample. This measurement represents a first constraint of the EBL at an effective redshift of ~1.8. We combine our results with prior EBL constraints and conclude that Fermi-LAT is instrumental to constrain the UV component of the EBL. We discuss the implications on existing empirical models of EBL evolution.

We present the second catalog of flaring gamma-ray sources (2FAV) detected with the Fermi All-sky Variability Analysis (FAVA), a tool that blindly searches for transients over the entire sky observed by the Large Area Telescope (LAT) on board the \textit{Fermi} Gamma-ray Space Telescope. With respect to the first FAVA catalog, this catalog benefits from a larger data set, the latest LAT data release (Pass 8), as well as from an improved analysis that includes likelihood techniques for a more precise localization of the transients. Applying this analysis on the first 7.4 years of \textit{Fermi} observations, and in two separate energy bands 0.1$-$0.8 GeV and 0.8$-$300 GeV, a total of 4547 flares has been detected with a significance greater than $6\sigma$ (before trials), on the time scale of one week. Through spatial clustering of these flares, 518 variable gamma-ray sources are identified. Likely counterparts, based on positional coincidence, have been found for 441 sources, mostly among the blazar class of active galactic nuclei. For 77 2FAV sources, no likely gamma-ray counterpart has been found. For each source in the catalog, we provide the time, location, and spectrum of each flaring episode. Studying the spectra of the flares, we observe a harder-when-brighter behavior for flares associated with blazars, with the exception of BL Lac flares detected in the low-energy band. The photon indexes of the flares are never significantly smaller than 1.5. For a leptonic model, and under the assumption of isotropy, this limit suggests that the spectrum of the freshly accelerated electrons is never harder than $p\sim$2.

The most powerful sources among the blazar family are MeV blazars. Often detected at $z>2$, they usually display high X- and \gm-ray luminosities, larger-than-average jet powers and black hole masses $\gtrsim 10^9 M_{\odot}$. In the present work we perform a multiwavelength study of three high redshift blazars: 3FGL J0325.5+2223 ($z=2.06$), 3FGL J0449.0+1121 ($z= 2.15$), and 3FGL J0453.2$-$2808 ($z=2.56$), analysing quasi simultaneous data from GROND, \swift-UVOT and XRT, \nustar, and \fermi-LAT. Our main focus is on the hard X-ray band recently unveiled by \nustar~(3$-$79 keV) where these objects show a hard spectrum which enables us to constrain the inverse Compton peak and the jet power. We found that all three targets resemble the most powerful blazars, with the synchrotron peak located in the sub-millimeter range and the inverse Compton peak in the MeV range, and therefore belong to the MeV blazar class. Using a simple one zone leptonic emission model to reproduce the spectral energy distributions, we conclude that a simple combination of synchrotron and accretion disk emission reproduces the infrared-optical spectra while the X-ray to \gm-ray part is well reproduced by the inverse Compton scattering of low energy photons supplied by the broad line region. The black hole masses for each of the three sources are calculated to be $\gtrsim 4 \times 10^{8} M_{\odot}$. The three studied sources have jet power at the level of, or beyond, the accretion luminosity.

The gamma-ray blazar OJ 287 was in a high activity state during December 2015 - February 2016. Coinciding with this high brightness state, we observed this source for photometry on 40 nights in R-band and for polarimetry on 9 epochs in UBVRI bands. During the period of our observations, the source brightness varied between $13.20 \pm 0.04$ to $14.98 \pm 0.04$ mag and the degree of polarization (P ) fluctuated between $6.0 \pm 0.3$% and $28.3 \pm 0.8$% in R-band. Focusing on intra-night optical variability (INOV), we find a duty cycle of about 71% using $\chi^2$-statistics, similar to that known for blazars. From INOV data, the shortest variability time scale is estimated to be $142 \pm 38$ min yielding a lower limit of the observed Doppler factor $\delta_0 = 1.17$, the magnetic field strength $B \le 3.8$ G and the size of the emitting region Rs < $2.28 \times 10^{14}$ cm. On inter-night timescales, a significant anti-correlation between R-band flux and P is found. The observed P at U-band is generally larger than that observed at longer wavelength bands suggesting a wavelength dependent polarization. Using V -band photometric and polarimetric data from Steward Observatory obtained during our monitoring period we find a varied correlation between P and V-band brightness. While an anticorrelation is seen between P and V -band mag at sometimes, no correlation is seen at other times, thereby, suggesting the presence of more than one short-lived shock components in the jet of OJ 287.

In a first systematic effort to characterize the intra-night optical variability (INOV) of different classes of narrow line Seyfert 1 (NLSy1) galaxies, we have carried out observations on a sample of radio-loud (RL) and radio-quiet (RQ) NLSy1 galaxies. The RL-NLSy1 galaxies are further divided into {\gamma}-ray loud (GL) and {\gamma}-ray quiet (GQ) NLSy1 galaxies. Our sample consists of four sets, each set consisting of a RQ-NLSy1, a GQ-NLSy1 and a GL-NLSy1 galaxy, closely matched in redshift and optical luminosity. Our observations on both RQ and GQ-NLSy1 galaxies consist of a total of 19 nights, whereas the data for GL-NLSy1 galaxies (18 nights) were taken from literature published earlier by us. This enabled us to do a comparison of the duty cycle (DC) of different classes of NLSy1 galaxies. Using power-enhanced F-test, with a variability threshold of 1%, we find DCs of about 55%, 39% and 0% for GL, GQ and RQ-NLSy1 galaxies respectively. The high DC and large amplitude of INOV (24.0 +/- 13.7%) shown by GL-NLSy1 galaxies relative to the other two classes might be due to their inner aligned relativistic jets having large bulk Lorentz factors. The null DC of RQ-NLSy1 galaxies could mean the presence of low power and/or largely misaligned jets in them. However, dividing RL-NLSy1 galaxies into low and high optical polarization sources, we find that sources with large polarization show somewhat higher DCs (69%) and amplitudes (29%) compared to those with low polarization. This points to a possible link between INOV and optical polarization.

Determining redshifts for BL Lacertae (BL Lac) objects using the traditional spectroscopic method is challenging due to the absence of strong emission lines in their optical spectra. We employ the photometric dropout technique to determine redshifts for this class of blazars using the combined 13 broad-band filters from Swift-UVOT and the multi-channel imager GROND at the MPG 2.2 m telescope at ESO's La Silla Observatory. The wavelength range covered by these 13 filters extends from far ultraviolet to the near-Infrared. We report results on 40 new Fermi detected BL Lacs with the photometric redshifts determinations for 5 sources, with 3FGL J1918.2-4110 being the most distant in our sample at z=2.16. Reliable upper limits are provided for 20 sources in this sample. Using the highest energy photons for these Fermi-LAT sources, we evaluate the consistency with the Gamma-ray horizon due to the extragalactic background light.

Detection of periodicity in the broad-band non-thermal emission of blazars has so far been proven to be elusive. However, there are a number of scenarios which could lead to quasi-periodic variations in blazar light curves. For example, orbital or thermal/viscous period of accreting matter around central supermassive black holes could, in principle, be imprinted in the multi-wavelength emission of small-scale blazar jets, carrying as such crucial information about plasma conditions within the jet launching regions. In this paper, we present the results of our time series analysis of $\sim 9.2$ year-long, and exceptionally well-sampled optical light curve of the BL Lac OJ 287. The study primarily uses the data from our own observations performed at the Mt. Suhora and Krak\'ow Observatories in Poland, and at the Athens Observatory in Greece. Additionally, SMARTS observations were used to fill in some of the gaps in the data. The Lomb-Scargle Periodogram and the Weighted Wavelet Z-transform methods were employed to search for the possible QPOs in the resulting optical light curve of the source. Both the methods consistently yielded possible quasi-periodic signal around the periods of $\sim 400$ and $\sim 800$ days, the former one with a significance (over the underlying colored noise) of $\geq 99\%$. A number of likely explanations for such are discussed, with a preference given to a modulation of the jet production efficiency by highly magnetized accretion disks. This supports the previous findings and the interpretation reported recently in the literature for OJ 287 and other blazar sources.

We present a multi-wavelength study of four high redshift blazars, S5 0014+81 ($z=3.37$), CGRaBS J0225+1846 ($z=2.69$), BZQ J1430+4205 ($z=4.72$), and 3FGL J1656.2$-$3303 ($z=2.40$), using the quasi-simultaneous data from {\it Swift}, {\it NuSTAR}, and {\it Fermi}-Large Area Telescope (LAT) and also the archival {\it XMM-Newton} observations. Other than 3FGL J1656.2$-$3303, none of the sources were known as $\gamma$-ray emitters and our analysis of $\sim$7.5 years of LAT data reveals the first time detection of the statistically significant $\gamma$-ray emission from CGRaBS J0225+1846. We generate the broadband spectral energy distributions (SED) of all the objects, centering at the epoch of {\it NuSTAR} observations and reproduce them using a one zone leptonic emission model. The optical$-$UV emission in all the objects can be explained by the radiation from the accretion disk, whereas, X-ray to $\gamma$-ray window of the SEDs are found to be dominated by the inverse Compton scattering off the broad line region photons. All of them host billion solar mass black holes. Comparing the accretion disk luminosity and the jet power of these sources with a large sample of blazars, we find them to occupy high disk luminosity-jet power regime. We also investigate the X-ray spectral properties of the sources in detail with a major focus on studying the causes of soft X-ray deficit, a feature generally seen in high redshift radio-loud quasars. We summarize that this feature could be explained based on the intrinsic curvature in the jet emission rather than due to external effects predicted in the earlier studies, such as host galaxy and/or warm absorption.

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.

The $\gamma$-ray loud narrow line Seyfert 1 ($\gamma$-NLSy1) galaxy PKS 1502+036 ($z=0.409$) exhibited its first $\gamma$-ray outburst on 2015 December 20. In the energy range of 0.1-300 GeV, the highest flux measured by {\it Fermi}-Large Area Telescope is (3.90 $\pm$ 1.52) $\times$ 10$^{-6}$ ph cm$^{-2}$ s$^{-1}$, which is the highest $\gamma$-ray flux ever detected from this object. The associated spectral shape is soft ($\Gamma_{0.1-300~{\rm GeV}}=2.57\pm0.17$) and this corresponds to an isotropic $\gamma$-ray luminosity of (1.2 $\pm$ 0.6) $\times$ 10$^{48}$ erg s$^{-1}$. We generate the broadband spectral energy distribution (SED) during the GeV flare and reproduce it using a one zone leptonic emission model. The optical-UV spectrum can be explained by a combination of synchrotron and the accretion disk emission, whereas, the X-ray to $\gamma$-ray SED can be satisfactorily reproduced by inverse-Compton scattering of thermal photons originated from the torus. The derived SED parameters hint for the increase in the bulk Lorentz factor as a major cause of the flare and the location of the emission region is estimated as outside the broad line region but still inside torus. A comparison of the GeV flaring SED of PKS 1502+036 with that of two other $\gamma$-NLSy1 galaxies, namely, 1H 0323+342 ($z=0.061$) and PMN J0948+0022 ($z=0.585$), and also with FSRQ 3C 279 ($z=0.536$) has led to the conclusion that the GeV flaring SEDs of $\gamma$-NLSy1 galaxies resemble with FSRQs and a major fraction of their bolometric luminosity is emitted at $\gamma$-ray energies.

The blazar 3C 279 exhibited a giant $\gamma$-ray outburst in 2013 December and 2014 April. Apart from the very fast $\gamma$-ray flux variability, the spectral nature of the flares were also found to vary significantly between these two flaring events. A prominent curvature in the $\gamma$-ray spectrum was noticed in 2014 April flare, on the other hand, the 2013 December displayed an extreme hardening of the spectrum. These observations, thus, put strong constraints on our understanding of the underlying particle acceleration mechanisms.

MeV blazars are a sub--population of the blazar family, exhibiting larger--than--average jet powers, accretion luminosities and black hole masses. Because of their extremely hard X--ray continua, these objects are best studied in the X-ray domain. Here, we report on the discovery by the $Fermi$ Large Area Telescope and subsequent follow-up observations with $NuSTAR$, $Swift$ and GROND of a new member of the MeV blazar family: PMN J0641$-$0320. Our optical spectroscopy provides confirmation that this is a flat--spectrum radio quasar located at a redshift of $z=1.196$. Its very hard $NuSTAR$ spectrum (power--law photon index of $\sim$1 up to $\sim$80 keV) indicates that the emission is produced via inverse Compton scattering off photons coming from outside the jet.The overall spectral energy distribution of PMN J0641$-$0320 is typical of powerful blazars and by reproducing it with a simple one-zone leptonic emission model we find the emission region to be located either inside the broad line region or within the dusty torus.

The blazar 3C 279 exhibited twin $\gamma$-ray flares of similar intensity in 2013 December and 2014 April. In this work, we present a detailed multi-wavelength analysis of the 2013 December flaring event. Multi-frequency observations reveal the uncorrelated variability patterns with X-ray and optical-UV fluxes peaking after the $\gamma$-ray maximum. The broadband spectral energy distribution (SED) at the peak of the $\gamma$-ray activity shows a rising $\gamma$-ray spectrum but a declining optical-UV flux. This observation along with the detection of uncorrelated variability behavior rules out the one-zone leptonic emission scenario. We, therefore, adopt two independent methodologies to explain the SED: a time dependent lepto-hadronic modeling and a two-zone leptonic radiative modeling approach. In the lepto-hadronic modeling, a distribution of electrons and protons subjected to a randomly orientated magnetic field produces synchrotron radiation. Electron synchrotron is used to explain the IR to UV emission while proton synchrotron emission is used to explain the high energy $\gamma$-ray emission. A combination of both electron synchrotron self Compton emission and proton synchrotron emission is used to explain the X-ray spectral break seen during the later stage of the flare. In the two-zone modeling, we assume a large emission region emitting primarily in IR to X-rays and $\gamma$-rays to come primarily from a fast moving compact emission region. We conclude by noting that within a span of 4 months, 3C 279 has shown the dominance of a variety of radiative processes over each other and this reflects the complexity involved in understanding the physical properties of blazar jets in general.

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.

An exceptional $\gamma$-ray outburst from 3C 279 is detected by {\it Fermi}-Large Area Telescope (LAT) in 2015 June. In the energy range of 0.1$-$300 GeV, the highest flux measured is (39.1$\pm$2.5) $\times$ 10$^{-6}$ \phflux, which is the highest $\gamma$-ray flux ever detected from 3C 279, exceeding the previous historically brightest flare observed by {\it EGRET} in 1996. The high activity period consists of three major flares with the last one being the brightest. All but one flares show a faster rise and slower decay pattern and at the peak of the activity, the $\gamma$-ray spectrum is found to show a clear signature of break/curvature. The obtained spectral parameters hint for the peak of the inverse Compton emission to lie in the LAT energy range (around $\sim$1 GeV) which is in contrast to that seen during the 2013 December and 2014 April $\gamma$-ray flares of 3C 279. From the $\gamma\gamma$ pair opacity arguments, the minimum Doppler factor is estimated to be 14 and the location of the $\gamma$-ray emitting region is found to be either at the outer edge of the broad line region or farther out from it.

A broadband study of the high redshift blazar S5 0836+71 (z = 2.172) is presented. Multi-frequency light curves show multiple episodes of X-ray and $\gamma$-ray flares, while optical-UV fluxes show little variations. During the GeV outburst, the highest $\gamma$-ray flux measured is (5.22 $\pm$ 1.10) $\times$ 10$^{-6}$ ph cm$^{-2}$ s$^{-1}$ in the range of 0.1-300 GeV, which corresponds to an isotropic $\gamma$-ray luminosity of (1.62 $\pm$ 0.44) $\times$ 10$^{50}$ erg s$^{-1}$, thereby making this as one of the most luminous $\gamma$-ray flare ever observed from any blazar. A fast $\gamma$-ray flux rising time of $\sim$3 hours is also noticed which is probably the first measurement of hour scale variability detected from a high redshift (z > 2) blazar. The various activity states of S5 0836+71 are reproduced under the assumption of single zone leptonic emission model. In all the states, the emission region is located inside the broad line region, and the optical-UV radiation is dominated by the accretion disk emission. The modeling parameters suggests the enhancement in bulk Lorentz factor as a primary cause of the $\gamma$-ray flare. The high X-ray activity with less variable $\gamma$-ray counterpart can be due to emission region to be located relatively closer to the black hole where the dominating energy density of the disk emission results in higher X-ray flux due to inverse-Compton scattering of disk photons.

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 detailed analysis of the gamma-ray flux variability and spectral properties of the five radio-loud narrow line Seyfert 1 (RL-NLSy1) galaxies, detected by the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope, namely 1H 0323+342, SBS 0846+513, PMN J0948+0022, PKS 1502+036, and PKS 2004$-$447. The first three sources show significant flux variations including the rapid variability of a few hours by 1H 0323+342. The average gamma-ray spectrum of 1H 0323+342 and PMN J0948+0022 shows deviation from a simple power law (PL) behavior, whereas for other three sources, PL model gives better fit. The spectra of 1H 0323+342, SBS 0846+513, and PMN J0948+0022, respectively in low, flaring, and moderately active states, show significant curvature. Such curvature in the gamma-ray spectrum of 1H 0323+342 and PMN J0948+0022 can be due to emission region located inside the broad line region (BLR) where the primary mechanism of the gamma-ray emission is inverse-Compton (IC) scattering of BLR photons occurring in the Klein-Nishina regime. The gamma-ray emission of SBS 0846+513 is explained by IC scattering of dusty torus photons which puts the emission region to be outside the BLR and thus under the Thomson regime. Therefore, the observed curvature of SBS 0846+513 could be intrinsic to the particle energy distribution. The presence of curvature in the gamma-ray spectrum and flux variability amplitudes of some of the RL-NLSy1 galaxies suggest that these sources could be akin to low/moderate jet power flat spectrum radio quasars.

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.

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.

The Large Area Telescope (LAT) onboard the {\it Fermi} Gamma Ray Space Telescope has detected \gamma-ray emission from some Narrow Line Seyfert 1(NLSy1) galaxies. This indicates the presence of relativistic jets in these sources similar to blazars and radio galaxies. In an attempt to have an idea of the intranight optical variability (INOV) characteristics of these \gamma-ray loud NLSy1 galaxies, we have carried out optical flux monitoring observations of three NLSy1 galaxies detected by {\it Fermi}/LAT: 1H 0323+342, PMN J0948+0022 and PKS 1502+036. These optical monitoring observations showed the presence of rapid optical flux variations in these sources. The intranight differential light curves of these sources have revealed flux variations on time scales of hours with amplitudes of variability > 3 percent for most of the time. On using the {\it F}-statistics to classify the variability nature of these sources, we obtained a duty cycle (DC) of INOV of ~85 percent. Alternatively, the more commonly used {\it C}-statistics gave a DC of INOV of ~57 percent. Such high DC of INOV are characteristics of the BL Lac class of AGN. The results of our monitoring observations thus indicate that there is similarity in the INOV nature of \gamma-ray loud NLSy1 galaxies and BL Lac objects, arguing strongly for the presence of relativistic jets aligned closely to the observers line of sight. Moreover, our dense monitoring observations on some of the nights have led to the clear detection of some mini-flares superimposed on the flux variations during the night over timescales as short as 12 minutes. The detection of short timescale flux variability in the sources studied here is clearly due to stronger time compression leading to the jets in these sources having large Doppler factors, similar to that of the inner jets of TeV blazars.