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.

We present an incremental version (4FGL-DR2, for Data Release 2) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first ten years of science data in the energy range from 50 MeV to 1 TeV, it uses the same analysis methods as the 4FGL catalog did for eight years of data. The spectral parameters, spectral energy distributions and associations are updated for all sources. Light curves are rebuilt for all sources with 1-year intervals (not 2-month intervals). Among the 5064 4FGL sources, 120 are formally below the detection threshold over 10 years (but are kept in the list), while 40 are newly associated. We report 723 new sources, mostly just above the detection threshold, among which two are considered identified and 342 have a plausible counterpart at other wavelengths.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

We present multi-frequency simultaneous VLBI radio observations of the flat spectrum radio quasar 1633+382 (4C~38.41) as part of the interferometric monitoring of gamma-ray bright active galactic nuclei (iMOGABA) program combined with additional observations in radio, optical, X-rays and $\gamma-$rays carried out between the period 2012 March - 2015 August. The monitoring of this source reveals a significant long-lived increase in its activity since approximately two years in the radio bands, which correlates with a similar increase on all other bands from sub-millimeter to $\gamma-$rays. A significant correlation is also found between radio fluxes and simultaneous spectral indices during this period. The study of the discrete correlation function (DCF) indicates time lags smaller than the $\sim40$ days uncertainties among both radio bands and also high-energy bands, and a time lag of $\sim$70 days, with $\gamma-$rays leading radio. We interpret that the high-energy and radio fluxes are arising from different emitting regions, located at $1\pm12$ and $40\pm13$ pc from the central engine respectively.

S5 0716+714 is a well known BL Lac object, one of the brightest and most active blazars. The discovery in the Very High Energy band (VHE, E > 100 GeV) by MAGIC happened in 2008. In January 2015 the source went through the brightest optical state ever observed, triggering MAGIC follow-up and a VHE detection with 13{\sigma} significance (ATel 6999). Rich multiwavelength coverage of the flare allowed us to construct the broad-band spectral energy distribution of S5 0716+714 during its brightest outburst. In this work we will present the preliminary analysis of MAGIC and Fermi-LAT data of the flaring activity in January and February 2015 for the HE (0.1 < HE < 300 GeV) and VHE band, together with radio (Mets\"ahovi, OVRO, VLBA, Effelsberg), sub-millimeter (SMA), optical (Tuorla, Perkins, Steward, AZT-8+ST7, LX-200, Kanata), X-ray and UV (Swift-XRT and UVOT), in the same time-window and discuss the time variability of the multiwavelength light curves during this impressive outburst.

A prominent outburst of the flat spectrum radio quasar 3C~454.3 was observed in 2014 June with the \emph{Fermi} Large Area Telescope. This outburst was characterized by a three-stage light-curve pattern---plateau, flare and post-flare---that occurred from 2014 May to July, in a similar pattern as observed during the exceptional outburst in 2010 November. The highest flux of the outburst reported in this paper occurred during 2014 June 7--29, showing a multiple-peak structure in the light-curves. The average flux in these 22 days was found to be $F[E > 100~\mathrm{MeV}] = (7.2 \pm 0.2) \times 10^{-6}$~ph~cm$^{-2}$~s$^{-1}$, with a spectral index, for a simple power law, of $\Gamma = 2.04 \pm 0.01$. That made this outburst the first $\gamma$-ray high state of 3C~454.3 ever to be detected by \emph{Fermi} with such a hard spectrum over several days. The highest flux was recorded on 2014 June 15, in a 3 hr bin, at MJD 56823.5625, at a level of $F[E > 100~\mathrm{MeV}] = (17.6 \pm 1.9) \times 10^{-6}$~ph~cm$^{-2}$~s$^{-1}$. The rise time of one of the short subflares was found to be $T_r= 1200 \pm 700$~s at MJD = 56827, when the flux increased from 4 to 12 $\times 10^{-6}$~ph~cm$^{-2}$~s$^{-1}$. Several photons above 20 GeV were collected during this outburst, including one at 45 GeV on MJD 56827, constraining the $\gamma$-ray emission region to be located close to the outer boundary of the broad-line region, leading to fast flux variability.

Using high-resolution radio imaging with VLBI techniques, the TANAMI program has been observing the parsec-scale radio jets of southern (declination south of -30{\deg}) gamma-ray bright AGN simultaneously with Fermi/LAT monitoring of their gamma-ray emission. We present the radio and gamma-ray properties of the TANAMI sources based on one year of contemporaneous TANAMI and Fermi/LAT data. A large fraction (72%) of the TANAMI sample can be associated with bright gamma-ray sources for this time range. Association rates differ for different optical classes with all BL Lacs, 76% of quasars and just 17% of galaxies detected by the LAT. Upper limits were established on the gamma-ray flux from TANAMI sources not detected by LAT. This analysis led to the identification of three new Fermi sources whose detection was later confirmed. The gamma-ray and radio luminosities are related by $L_\gamma \propto L_r^{0.89+-0.04}$. The brightness temperatures of the radio cores increase with the average gamma-ray luminosity, and the presence of brightness temperatures above the inverse Compton limit implies strong Doppler boosting in those sources. The undetected sources have lower gamma/radio luminosity ratios and lower contemporaneous brightness temperatures. Unless the Fermi/LAT-undetected blazars are strongly gamma-ray-fainter than the Fermi/LAT-detected ones, their gamma-ray luminosity should not be significantly lower than the upper limits calculated here.

We report for the first time a gamma-ray and multi-wavelength nearly-periodic oscillation in an active galactic nucleus. Using the Fermi Large Area Telescope (LAT) we have discovered an apparent quasi-periodicity in the gamma-ray flux (E >100 MeV) from the GeV/TeV BL Lac object PG 1553+113. The marginal significance of the 2.18 +/-0.08 year-period gamma-ray cycle is strengthened by correlated oscillations observed in radio and optical fluxes, through data collected in the OVRO, Tuorla, KAIT, and CSS monitoring programs and Swift UVOT. The optical cycle appearing in ~10 years of data has a similar period, while the 15 GHz oscillation is less regular than seen in the other bands. Further long-term multi-wavelength monitoring of this blazar may discriminate among the possible explanations for this quasi-periodicity.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi-LAT (3LAC) is presented. It is based on the third Fermi-LAT catalog (3FGL) of sources detected between 100 MeV and 300 GeV with a Test Statistic (TS) greater than 25, between 2008 August 4 and 2012 July 31. The 3LAC includes 1591 AGNs located at high Galactic latitudes (|b|>10{\deg}), a 71% increase over the second catalog based on 2 years of data. There are 28 duplicate associations, thus 1563 of the 2192 high-latitude gamma-ray sources of the 3FGL catalog are AGNs. Most of them (98%) are blazars. About half of the newly detected blazars are of unknown type, i.e., they lack spectroscopic information of sufficient quality to determine the strength of their emission lines. Based on their gamma-ray spectral properties, these sources are evenly split between flat-spectrum radio quasars (FSRQs) and BL~Lacs. The most abundant detected BL~Lacs are of the high-synchrotron-peaked (HSP) type. About 50% of the BL~Lacs have no measured redshifts. A few new rare outliers (HSP-FSRQs and high-luminosity HSP BL~Lacs) are reported. The general properties of the 3LAC sample confirm previous findings from earlier catalogs. The fraction of 3LAC blazars in the total population of blazars listed in BZCAT remains non-negligible even at the faint ends of the BZCAT-blazar radio, optical and X-ray flux distributions, which is a clue that even the faintest known blazars could eventually shine in gamma rays at LAT-detection levels. The energy-flux distributions of the different blazar populations are in good agreement with extrapolation from earlier catalogs.

The third catalog of active galactic nuclei (AGNs) detected by the Fermi Large Area Telescope (3LAC) is presented. It is based on the third catalog (3FGL,\cite{3FGL}) of sources detected with a test statistic greater than 25, using the first 4 years of data. The 3LAC includes 1591 AGNs located at high ($|b|>10^\circ$) Galactic latitudes (with 28 duplicate associations, thus corresponding to 1563 gamma-ray sources among 2192 sources in the 3FGL catalog), providing $71\%$ more sources with respect to the 2FGL. Various properties, such as gamma-ray fluxes and photon power law spectral indices, redshifts, gamma-ray luminosities, variability, and their correlations are presented and discussed for the different blazar and non-blazar classes.

Fermi-LAT analyses show that the gamma-ray photon spectral indices Gamma_gamma of a large sample of blazars correlate with the vFv peak synchrotron frequency v_s according to the relation Gamma_gamma = d - k log v_s. The same function, with different constants d and k, also describes the relationship between Gamma_gamma and peak Compton frequency v_C. This behavior is derived analytically using an equipartition blazar model with a log-parabola description of the electron energy distribution (EED). In the Thomson regime, k = k_EC = 3b/4 for external Compton processes and k = k_SSC = 9b/16 for synchrotron self-Compton (SSC) processes, where b is the log-parabola width parameter of the EED. The BL Lac object Mrk 501 is fit with a synchrotron/SSC model given by the log-parabola EED, and is best fit away from equipartition. Corrections are made to the spectral-index diagrams for a low-energy power-law EED and departures from equipartition, as constrained by absolute jet power. Analytic expressions are compared with numerical values derived from self-Compton and external Compton scattered gamma-ray spectra from Ly alpha broad-line region and IR target photons. The Gamma_gamma vs. v_s behavior in the model depends strongly on b, with progressively and predictably weaker dependences on gamma-ray detection range, variability time, and isotropic gamma-ray luminosity. Implications for blazar unification and blazars as ultra-high energy cosmic-ray sources are presented. Arguments by Ghisellini et al. (2014) that the jet power exceeds the accretion luminosity depend on the doubtful assumption that we are viewing at the Doppler angle.

Flat spectrum radio quasars (FSRQs) are bright active galactic nuclei surrounded by gas clouds within a UV-visible intense radiation field that form the so-called broad line region (BLR). These objects emit relativistic jets from a region close to the central supermassive black hole and through the BLR. The Fermi-Large Area Telescope (Fermi-LAT) is sensitive to gamma-ray photons from $\sim$30 MeV to more than 300 GeV. We have performed spectral analysis of bright FSRQs in a 5.5 year (2008-2014) data sample collected by Fermi-LAT, using the new Pass 8 event selection and instrument response function. Also, our study of flaring episodes in a limited time range brings interesting results while compared to the full 5.5 year data samples.

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope routinely detects the highly dust-absorbed, reddened, and MeV-peaked flat spectrum radio quasar PKS 1830-211 (z=2.507). Its apparent isotropic gamma-ray luminosity (E>100 MeV) averaged over $\sim$ 3 years of observations and peaking on 2010 October 14/15 at 2.9 X 10^{50} erg s^{-1}, makes it among the brightest high-redshift Fermi blazars. No published model with a single lens can account for all of the observed characteristics of this complex system. Based on radio observations, one expects time delayed variability to follow about 25 days after a primary flare, with flux about a factor 1.5 less. Two large gamma-ray flares of PKS 1830-211 have been detected by the LAT in the considered period and no substantial evidence for such a delayed activity was found. This allows us to place a lower limit of about 6 on the gamma rays flux ratio between the two lensed images. Swift XRT observations from a dedicated Target of Opportunity program indicate a hard spectrum and with no significant correlation of X-ray flux with the gamma-ray variability. The spectral energy distribution can be modeled with inverse Compton scattering of thermal photons from the dusty torus. The implications of the LAT data in terms of variability, the lack of evident delayed flare events, and different radio and gamma-ray flux ratios are discussed. Microlensing effects, absorption, size and location of the emitting regions, the complex mass distribution of the system, an energy-dependent inner structure of the source, and flux suppression by the lens galaxy for one image path may be considered as hypotheses for understanding our results.

We have been using the 0.76-m Katzman Automatic Imaging Telescope (KAIT) at Lick Observatory to optically monitor a sample of 157 blazars that are bright in gamma rays, being detected with high significance ($\ge 10\sigma$) in one year by the Large Area Telescope (LAT) on the {\it Fermi Gamma-ray Space Telescope}. We attempt to observe each source on a 3-day cadence with KAIT, subject to weather and seasonal visibility. The gamma-ray coverage is essentially continuous. KAIT observations extend over much of the 5-year {\it Fermi} mission for several objects, and most have $>100$ optical measurements spanning the last three years. These blazars (flat-spectrum radio quasars and BL~Lac objects) exhibit a wide range of flaring behavior. Using the discrete correlation function (DCF), here we search for temporal relationships between optical and gamma-ray light curves in the 40 brightest sources in hopes of placing constraints on blazar acceleration and emission zones. We find strong optical--gamma-ray correlation in many of these sources at time delays of $\sim 1$ to $\sim 10$ days, ranging between $-40$ and +30 days. A stacked average DCF of the 40 sources verifies this correlation trend, with a peak above 99% significance indicating a characteristic time delay consistent with 0 days. These findings strongly support the widely accepted leptonic models of blazar emission. However, we also find examples of apparently uncorrelated flares (optical flares with no gamma-ray counterpart and gamma-ray flares with no optical counterpart) that challenge simple, one-zone models of blazar emission. Moreover, we find that flat-spectrum radio quasars tend to have gamma rays leading the optical, while intermediate and high synchrotron peak blazars with the most significant peaks have smaller lags/leads.

The {\gamma}-ray sky can be decomposed into individually detected sources, diffuse emission attributed to the interactions of Galactic cosmic rays with gas and radiation fields, and a residual all-sky emission component commonly called the isotropic diffuse {\gamma}-ray background (IGRB). The IGRB comprises all extragalactic emissions too faint or too diffuse to be resolved in a given survey, as well as any residual Galactic foregrounds that are approximately isotropic. The first IGRB measurement with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) used 10 months of sky-survey data and considered an energy range between 200 MeV and 100 GeV. Improvements in event selection and characterization of cosmic-ray backgrounds, better understanding of the diffuse Galactic emission, and a longer data accumulation of 50 months, allow for a refinement and extension of the IGRB measurement with the LAT, now covering the energy range from 100 MeV to 820 GeV. The IGRB spectrum shows a significant high-energy cutoff feature, and can be well described over nearly four decades in energy by a power law with exponential cutoff having a spectral index of $2.32\pm0.02$ and a break energy of $(279\pm52)$ GeV using our baseline diffuse Galactic emission model. The total intensity attributed to the IGRB is $(7.2\pm0.6) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$ sr$^{-1}$ above 100 MeV, with an additional $+15$%/$-30$% systematic uncertainty due to the Galactic diffuse foregrounds.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

Blazar spectral models generally have numerous unconstrained parameters, leading to ambiguous values for physical properties like Doppler factor delta or fluid magnetic field B'. To help remedy this problem, a few modifications of the standard leptonic blazar jet scenario are considered. First, a log-parabola function for the electron distribution is used. Second, analytic expressions relating energy loss and kinematics to blazar luminosity and variability, written in terms of equipartition parameters, imply delta, B', and the principal electron Lorentz factor gamma'_pk. The external radiation field in a blazar is approximated by Ly alpha radiation from the broad line region (BLR) and ~0.1 eV infrared radiation from a dusty torus. When used to model 3C 279 SEDs from 2008 and 2009 reported by Hayashida et al. (2012), we derive delta ~ 20-30, B' ~ few G, and total (IR + BLR) external radiation field energy densities u ~ 0.01 - 0.001 erg/cm^3, implying an origin of the gamma-ray emission site in 3C 279 at the outer edges of the BLR. This is consistent with the gamma-ray emission site being located at a distance R <~ Gamma^2 c t_{var} ~ 0.1 (\Gamma/30)^2 (t_{var}/10^4 s) pc from the black hole powering 3C 279's jets, where t_{var} is the variability time scale of the radiation in the source frame, and at farther distances for narrow-jet and magnetic_reconnection models. Excess >~ 5 GeV gamma-ray emission observed with Fermi LAT from 3C 279 challenge the model, opening the possibility of hadronic origins of the emission. For low hadronic content, absolute jet powers of ~10% of the Eddington luminosity are calculated.

Using data from the Fermi Large Area Telescope (LAT), we report the first clear gamma-ray measurement of a delay between flares from the gravitationally lensed images of a blazar. The delay was detected in B0218+357, a known double-image lensed system, during a period of enhanced gamma-ray activity with peak fluxes consistently observed to reach >20-50 times its previous average flux. An auto-correlation function analysis identified a delay in the gamma-ray data of 11.46 +/- 0.16 days (1 sigma) that is ~1 day greater than previous radio measurements. Considering that it is beyond the capabilities of the LAT to spatially resolve the two images, we nevertheless decomposed individual sequences of superposing gamma-ray flares/delayed emissions. In three such ~8-10 day-long sequences within a ~4-month span, considering confusion due to overlapping flaring emission and flux measurement uncertainties, we found flux ratios consistent with ~1, thus systematically smaller than those from radio observations. During the first, best-defined flare, the delayed emission was detailed with a Fermi pointing, and we observed flux doubling timescales of ~3-6 hrs implying as well extremely compact gamma-ray emitting regions.

We present a \gamma-ray photon flux and spectral variability study of the flat-spectrum radio quasar 3C 273 over a rapid flaring activity period between September 2009 to April 2010. Five major flares are observed in the source during this period. The most rapid flare observed in the source has a flux doubling time of 1.1 hr. The rapid \gamma-ray flares allow us to constrain the location and size of the \gamma-ray emission region in the source. The \gamma \gamma-opacity constrains the Doppler factor, $\delta_{\gamma} \geq$ 10 for the highest energy (15 GeV) photon observed by the {\it Fermi}-Large Area Telescope (LAT). Causality arguments constrain the size of the emission region to 1.6$\times 10^{15}$ cm. The \gamma-ray spectra measured over this period show clear deviations from a simple power law with a break in 1-2 GeV energy range. We discuss possible explanations for the origin of the \gamma-ray spectral breaks. Our study suggests that the \gamma-ray emission region in 3C 273 is located within the broad line region ($<$1.6 pc). The spectral behavior and temporal characteristics of the individual flares indicate the presence of multiple shock scenarios at the base of the jet.

Observations performed with the Fermi-LAT telescope have revealed the presence of a spectral break in the GeV spectrum of flat-spectrum radio quasars (FSRQs) and other low- and intermediate-synchrotron peaked blazars. We propose that this feature can be explained by Compton scattering of broad-line region (BLR) photons by a non-thermal population of electrons described by a log-parabolic function. We consider in particular a scenario in which the energy densities of particles, magnetic field, and soft photons in the emitting region are close to equipartition. We show that this model can satisfactorily account for the overall spectral energy distribution of the FSRQ 3C 454.3, reproducing the GeV spectal cutoff due to Klein-Nishina effects and a curving electron distribution.

The GeV observations by Fermi-LAT give us the opportunity to characterize the high-energy emission (100 MeV - 300 GeV) variability properties of the BL Lac object S5 0716+714. In this study, we performed flux and spectral analysis of more than 3 year long (August 2008 to April 2012) Fermi-LAT data of the source. During this period, the source exhibits two different modes of flux variability with characteristic timescales of ~75 and ~140 days, respectively. We also notice that the flux variations are characterized by a weak spectral hardening. The GeV spectrum of the source shows a clear deviation from a simple power law, and is better explained by a broken power law. Similar to other bright Fermi blazars, the break energy does not vary with the source flux during the different activity states. We discuss several possible scenarios to explain the observed spectral break.

Densely time sampled multi-frequency flux measurements of the extreme BL Lac object S5 0716+714 over the past three years allow us to study its broad-band variability, and the detailed underlying physics, with emphasis on the location and size of the emitting regions and the evolution with time. We study the characteristics of some prominent mm-/gamma-ray flares in the context of the shock-in-jet model and investigate the location of the high energy emission region. The rapid rise and decay of the radio flares is in agreement with the formation of a shock and its evolution, if a geometrical variation is included in addition to intrinsic variations of the source. We find evidence for a correlation between flux variations at gamma-ray and radio frequencies. A two month time-delay between gamma-ray and radio flares indicates a non-cospatial origin of gamma-rays and radio flux variations in S5 0716+714.

We present the results of a series of radio, optical, X-ray and gamma-ray observations of the BL Lac object S50716+714 carried out between April 2007 and January 2011. The multi-frequency observations were obtained using several ground and space based facilities. The intense optical monitoring of the source reveals faster repetitive variations superimposed on a long-term variability trend at a time scale of ~350 days. Episodes of fast variability recur on time scales of ~ 60-70 days. The intense and simultaneous activity at optical and gamma-ray frequencies favors the SSC mechanism for the production of the high-energy emission. Two major low-peaking radio flares were observed during this high optical/gamma-ray activity period. The radio flares are characterized by a rising and a decaying stage and are in agreement with the formation of a shock and its evolution. We found that the evolution of the radio flares requires a geometrical variation in addition to intrinsic variations of the source. Different estimates yield a robust and self-consistent lower limits of \delta > 20 and equipartition magnetic field B_eq > 0.36 G. Causality arguments constrain the size of emission region \theta < 0.004 mas. We found a significant correlation between flux variations at radio frequencies with those at optical and gamma-rays. The optical/GeV flux variations lead the radio variability by ~65 days. The longer time delays between low-peaking radio outbursts and optical flares imply that optical flares are the precursors of radio ones. An orphan X-ray flare challenges the simple, one-zone emission models, rendering them too simple. Here we also describe the spectral energy distribution modeling of the source from simultaneous data taken through different activity periods.

We have performed an analysis of the diffuse gamma-ray emission with the Fermi Large Area Telescope in the Milky Way Halo region searching for a signal from dark matter annihilation or decay. In the absence of a robust dark matter signal, constraints are presented. We consider both gamma rays produced directly in the dark matter annihilation/decay and produced by inverse Compton scattering of the e+e- produced in the annihilation/decay. Conservative limits are derived requiring that the dark matter signal does not exceed the observed diffuse gamma-ray emission. A second set of more stringent limits is derived based on modeling the foreground astrophysical diffuse emission using the GALPROP code. Uncertainties in the height of the diffusive cosmic-ray halo, the distribution of the cosmic-ray sources in the Galaxy, the index of the injection cosmic-ray electron spectrum and the column density of the interstellar gas are taken into account using a profile likelihood formalism, while the parameters governing the cosmic-ray propagation have been derived from fits to local cosmic-ray data. The resulting limits impact the range of particle masses over which dark matter thermal production in the early Universe is possible, and challenge the interpretation of the PAMELA/Fermi-LAT cosmic ray anomalies as annihilation of dark matter.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

Recent detections of the starburst galaxies M82 and NGC 253 by gamma-ray telescopes suggest that galaxies rapidly forming massive stars are more luminous at gamma-ray energies compared to their quiescent relatives. Building upon those results, we examine a sample of 69 dwarf, spiral, and luminous and ultraluminous infrared galaxies at photon energies 0.1-100 GeV using 3 years of data collected by the Large Area Telescope (LAT) on the \textit{Fermi Gamma-ray Space Telescope} (\textit{Fermi}). Measured fluxes from significantly detected sources and flux upper limits for the remaining galaxies are used to explore the physics of cosmic rays in galaxies. We find further evidence for quasi-linear scaling relations between gamma-ray luminosity and both radio continuum luminosity and total infrared luminosity which apply both to quiescent galaxies of the Local Group and low-redshift starburst galaxies (conservative $P$-values $\lesssim0.05$ accounting for statistical and systematic uncertainties). The normalizations of these scaling relations correspond to luminosity ratios of $\log(L_{0.1-100 \rm{GeV}}/L_{1.4 \rm{GHz}}) = 1.7 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ and $\log(L_{0.1-100 \rm{GeV}}/L_{8-1000 \mu\rm{m}}) = -4.3 \pm 0.1_{\rm (statistical)} \pm 0.2_{\rm (dispersion)}$ for a galaxy with a star formation rate of 1 $M_{\odot}$ yr$^{-1}$, assuming a Chabrier initial mass function. Using the relationship between infrared luminosity and gamma-ray luminosity, the collective intensity of unresolved star-forming galaxies at redshifts $0<z<2.5$ above 0.1 GeV is estimated to be 0.4-2.4 $\times 10^{-6}$ ph cm$^{-2}$ s$^{-1}$ sr$^{-1}$ (4-23% of the intensity of the isotropic diffuse component measured with the LAT). We anticipate that $\sim10$ galaxies could be detected by their cosmic-ray induced gamma-ray emission during a 10-year \textit{Fermi} mission.

The Second Catalog of Blazars and other Active Galactic Nuclei detected by the Fermi/LAT (2LAC) includes about 1100 sources, 886 of which comprise the Clean Sample. The general properties of the different populations of sources classified according to the strength of their emission lines (FSRQs, BL Lacs) or the estimated position of the synchrotron peak are reviewed.

We present a method enabling the creation of constant-uncertainty/constant-significance light curves with the data of the Fermi-Large Area Telescope (LAT). The adaptive-binning method enables more information to be encapsulated within the light curve than with the fixed-binning method. Although primarily developed for blazar studies, it can be applied to any sources. This method allows the starting and ending times of each interval to be calculated in a simple and quick way during a first step. The reported mean flux and spectral index (assuming the spectrum is a power-law distribution) in the interval are calculated via the standard LAT analysis during a second step. The absence of major caveats associated with this method has been established by means of Monte-Carlo simulations. We present the performance of this method in determining duty cycles as well as power-density spectra relative to the traditional fixed-binning method.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

The high-frequency peaked BL Lac object PKS 2005-489 was the target of a multi-wavelength campaign with simultaneous observations in the TeV gamma-ray (H.E.S.S.), GeV gamma-ray (Fermi/LAT), X-ray (RXTE, Swift), UV (Swift) and optical (ATOM, Swift) bands. This campaign was carried out during a high flux state in the synchrotron regime. The flux in the optical and X-ray bands reached the level of the historical maxima. The hard GeV spectrum observed with Fermi/LAT connects well to the very high energy (VHE, E>100GeV) spectrum measured with H.E.S.S. with a peak energy between ~5 and 500 GeV. Compared to observations with contemporaneous coverage in the VHE and X-ray bands in 2004, the X-ray flux was ~50 times higher during the 2009 campaign while the TeV gamma-ray flux shows marginal variation over the years. The spectral energy distribution during this multi-wavelength campaign was fit by a one zone synchrotron self-Compton model with a well determined cutoff in X-rays. The parameters of a one zone SSC model are inconsistent with variability time scales. The variability behaviour over years with the large changes in synchrotron emission and small changes in the inverse Compton emission does not warrant an interpretation within a one-zone SSC model despite an apparently satisfying fit to the broadband data in 2009.

Some questions raised by Fermi-LAT data about blazars are summarized, along with attempts at solutions within the context of leptonic models. These include both spectral and statistical questions, including the origin of the GeV breaks in low-synchrotron peaked blazars, the location of the gamma-ray emission sites, the correlations in the spectral energy distributions with luminosity, and the difficulty of synchrotron/SSC models to fit the spectra of some TeV blazars.

We investigate the Fermi LAT gamma-ray and 15 GHz VLBA radio properties of a joint gamma-ray- and radio-selected sample of AGNs obtained during the first 11 months of the Fermi mission (2008 Aug 4 - 2009 Jul 5). Our sample contains the brightest 173 AGNs in these bands above declination -30 deg. during this period, and thus probes the full range of gamma-ray loudness (gamma-ray to radio band luminosity ratio) in the bright blazar population. The latter quantity spans at least four orders of magnitude, reflecting a wide range of spectral energy distribution (SED) parameters in the bright blazar population. The BL Lac objects, however, display a linear correlation of increasing gamma-ray loudness with synchrotron SED peak frequency, suggesting a universal SED shape for objects of this class. The synchrotron self-Compton model is favored for the gamma-ray emission in these BL Lacs over external seed photon models, since the latter predict a dependence of Compton dominance on Doppler factor that would destroy any observed synchrotron SED peak - gamma-ray loudness correlation. The high-synchrotron peaked (HSP) BL Lac objects are distinguished by lower than average radio core brightness temperatures, and none display large radio modulation indices or high linear core polarization levels. No equivalent trends are seen for the flat-spectrum radio quasars (FSRQ) in our sample. Given the association of such properties with relativistic beaming, we suggest that the HSP BL Lacs have generally lower Doppler factors than the lower-synchrotron peaked BL Lacs or FSRQs in our sample.

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the Active Galactic Nuclei (AGN) detected by Fermi during its first year of operation, with the largest datasets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the VLA and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate-data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the cm radio and the broad band (E>100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability <1e-7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6e-6 to 9.0e-8. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability <1e-7) for both FSRQs and BL Lacs separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

Recent claims that the strength B_IGMF of the intergalactic magnetic field (IGMF) is >~ 1e-15 G are based on upper limits to the expected cascade flux in the GeV band produced by blazar TeV photons absorbed by the extragalactic background light. This limit depends on an assumption that the mean blazar TeV flux remains constant on timescales >~2 (B_ IGMF/1e-18 G)^2 / (E/{10 GeV})^2 yr for an IGMF coherence length ~ 1 Mpc, where E is the measured photon energy. Restricting TeV activity of 1ES 0229+200 to ~3 -- 4 years during which the source has been observed leads to a more robust lower limit of B_IGMF >~ 1e-18 G, which can be larger by an order of magnitude if the intrinsic source flux above ~5 -- 10 TeV from 1ES 0229+200 is strong.

In this paper we report on the two-year-long Fermi-LAT observation of the peculiar blazar 4C +21.35 (PKS 1222+216). This source was in a quiescent state from the start of science operations of the Fermi Gamma-ray Space Telescope in 2008 August until 2009 September, and then became more active, with gradually increasing flux and some moderately-bright flares. In 2010 April and June, 4C +21.35 underwent a very strong GeV outburst composed of several major flares characterized by rise and decay timescales of the order of a day. During the outburst, the GeV spectra of 4C +21.35 displayed a broken power-law form with spectral breaks observed near 1-3 GeV photon energies. We demonstrate that, at least during the major flares, the jet in 4C +21.35 carried a total kinetic luminosity comparable to the total accretion power available to feed the outflow. We also discuss the origin of the break observed in the flaring spectra of 4C +21.35. We show that, in principle, a model involving annihilation of the GeV photons on the He II Lyman recombination continuum and line emission of "broad line region" clouds may account for such. However, we also discuss the additional constraint provided by the detection of 4C +21.35 at 0.07-0.4 TeV energies by the MAGIC telescope, which coincided with one of the GeV flares of the source. We argue that there are reasons to believe that the $lesssim$,TeV emission of 4C +21.35 (as well as the GeV emission of the source, if co-spatial), is not likely to be produced inside the broad line region zone of highest ionization ($sim 10^{17}$,cm from the nucleus), but instead originates further away from the active center, namely around the characteristic scale of the hot dusty torus surrounding the 4C +21.35 nucleus ($sim 10^{19}$,cm).

The high-frequency-peaked BL Lacertae object RGB J0710+591 was observed in the very high-energy (VHE; E > 100 GeV) wave band by the VERITAS array of atmospheric Cherenkov telescopes. The observations, taken between 2008 December and 2009 March and totaling 22.1 hr, yield the discovery of VHE gamma rays from the source. RGB J0710+591 is detected at a statistical significance of 5.5 standard deviations (5.5{\sigma}) above the background, corresponding to an integral flux of (3.9 +/- 0.8) x 10-12 cm-2 s-1 (3% of the Crab Nebula's flux) above 300 GeV. The observed spectrum can be fit by a power law from 0.31 to 4.6 TeV with a photon spectral index of 2.69 +/- 0.26stat +/- 0.20sys. These data are complemented by contemporaneous multiwavelength data from the Fermi Large Area Telescope, the Swift X-ray Telescope, the Swift Ultra-Violet and Optical Telescope, and the Michigan-Dartmouth-MIT observatory. Modeling the broadband spectral energy distribution (SED) with an equilibrium synchrotron self-Compton model yields a good statistical fit to the data. The addition of an external-Compton component to the model does not improve the fit nor brings the system closer to equipartition. The combined Fermi and VERITAS data constrain the properties of the high-energy emission component of the source over 4 orders of magnitude and give measurements of the rising and falling sections of the SED.

The first published Fermi large area telescope (Fermi-LAT) measurement of the isotropic diffuse gamma-ray emission is in good agreement with a single power law, and is not showing any signature of a dominant contribution from dark matter sources in the energy range from 20 to 100 GeV. We use the absolute size and spectral shape of this measured flux to derive cross section limits on three types of generic dark matter candidates: annihilating into quarks, charged leptons and monochromatic photons. Predicted gamma-ray fluxes from annihilating dark matter are strongly affected by the underlying distribution of dark matter, and by using different available results of matter structure formation we assess these uncertainties. We also quantify how the dark matter constraints depend on the assumed conventional backgrounds and on the Universe's transparency to high-energy gamma-rays. In reasonable background and dark matter structure scenarios (but not in all scenarios we consider) it is possible to exclude models proposed to explain the excess of electrons and positrons measured by the Fermi-LAT and PAMELA experiments. Derived limits also start to probe cross sections expected from thermally produced relics (e.g. in minimal supersymmetry models) annihilating predominantly into quarks. For the monochromatic gamma-ray signature, the current measurement constrains only dark matter scenarios with very strong signals.

Dark matter (DM) particle annihilation or decay can produce monochromatic $\gamma$-rays readily distinguishable from astrophysical sources. $\gamma$-ray line limits from 30 GeV to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a $\gamma$-ray line analysis, and integrated over most of the sky. We obtain $\gamma$-ray line flux upper limits in the range $0.6-4.5\times 10^{-9}\mathrm{cm}^{-2}\mathrm{s}^{-1}$, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

(Abridged) We have conducted a detailed investigation of the broad-band spectral properties of the \gamma-ray selected blazars of the Fermi LAT Bright AGN Sample (LBAS). By combining our accurately estimated Fermi gamma-ray spectra with Swift, radio, infra-red, optical and other hard X-ray/gamma-ray data, collected within three months of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous Spectral Energy Distributions (SED) for 48 LBAS blazars.The SED of these gamma-ray sources is similar to that of blazars discovered at other wavelengths, clearly showing, in the usual Log $\nu $ - Log $\nu$ F$_\nu$ representation, the typical broad-band spectral signatures normally attributed to a combination of low-energy synchrotron radiation followed by inverse Compton emission of one or more components. We have used these SEDs to characterize the peak intensity of both the low and the high-energy components. The results have been used to derive empirical relationships that estimate the position of the two peaks from the broad-band colors (i.e. the radio to optical and optical to X-ray spectral slopes) and from the gamma-ray spectral index. Our data show that the synchrotron peak frequency $\nu_p^S$ is positioned between 10$^{12.5}$ and 10$^{14.5}$ Hz in broad-lined FSRQs and between $10^{13}$ and $10^{17}$ Hz in featureless BL Lacertae objects.We find that the gamma-ray spectral slope is strongly correlated with the synchrotron peak energy and with the X-ray spectral index, as expected at first order in synchrotron - inverse Compton scenarios. However, simple homogeneous, one-zone, Synchrotron Self Compton (SSC) models cannot explain most of our SEDs, especially in the case of FSRQs and low energy peaked (LBL) BL Lacs. (...)

The recent detection by the Fermi gamma-ray space telescope of high-energy gamma-rays from the radio galaxy NGC 1275 makes the observation of the very high energy (VHE: E > 100 GeV) part of its broadband spectrum particularly interesting, especially for the understanding of active galactic nuclei (AGN) with misaligned multi-structured jets. The radio galaxy NGC 1275 was recently observed by VERITAS at energies above 100 GeV for about 8 hours. No VHE gamma-ray emission was detected by VERITAS from NGC 1275. A 99% confidence level upper limit of 2.1% of the Crab Nebula flux level is obtained at the decorrelation energy of approximately 340 GeV, corresponding to 19% of the power-law extrapolation of the Fermi Large Area Telescope (LAT) result.

We have measured the gamma-ray fluxes of the blazars Mrk 421 and Mrk 501 in the energy range between 50 and 350 GeV (1.2 to 8.3 x 10^25 Hz). The detector, called CELESTE, used first 40, then 53 heliostats of the former solar facility "Themis" in the French Pyrenees to collect Cherenkov light generated in atmospheric particle cascades. The signal from Mrk 421 is often strong. We compare its flux with previously published multi-wavelength studies and infer that we are straddling the high energy peak of the spectral energy distribution. The signal from Mrk 501 in 2000 was weak (3.4 sigma). We obtain an upper limit on the flux from 1ES 1426+428 of less than half that of the Crab flux near 100 GeV. The data analysis and understanding of systematic biases have improved compared to previous work, increasing the detector's sensitivity.

If dark matter is made of neutralinos, annihilation of such Majorana particles should produce high energy cosmic rays, especially in galaxy halo high density regions like galaxy centres. M31 (Andromeda) is our nearest neighbour spiral galaxy, and both its high mass and its low distance make it a source of interest for the indirect search for dark matter through gamma-ray detection. The ground based atmospheric Cherenkov telescope CELESTE observed M31 from 2001 to 2003, in the mostly unexplored energy range 50-500 GeV. These observations provide an upper limit on the flux above 50 GeV around $10^{-10}\rm{cm}^{-2}\rm{s}^{-1}$ in the frame of supersymmetric dark matter, and more generally on any gamma emission from M31.