Axion-like particles (ALPs), predicted in theories beyond the Standard Model, can have observational effects on the transparency of the Universe to $\gamma$ rays in the presence of magnetic fields. In this work, we search for effects compatible with the existence of ALPs with 80 months of data from the ${\it Fermi}$ Large Area Telescope, by comparing the distributions of observed highest energy photons from sources beyond redshifts of z $\geq$ 0.1 with theoretical predictions in the presence of ALPs. We find no evidence for an increased $\gamma$-ray transparency due to ALPs and therefore we set limits on the ALPs parameters assuming a value of the intergalactic magnetic field strength of 1 nG. Photon-ALP couplings above $10^{-11}$ GeV$^{-1}$ are excluded for ALP masses $m_{a}$ $\lesssim3.0$ neV. These limits exclude a region of the parameter space not covered by other $\gamma$-ray telescopes and are compatible with constraints imposed by other experiments.

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.

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

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.

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.

On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

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.

We present a measurement of the cosmic-ray electron+positron spectrum between 7 GeV and 2 TeV performed with almost seven years of data collected with the Fermi Large Area Telescope. We find that the spectrum is well fit by a broken power law with a break energy at about 50 GeV. Above 50 GeV, the spectrum is well described by a single power law with a spectral index of $3.07 \pm 0.02 \; (\text{stat+syst}) \pm 0.04 \; (\text{energy measurement})$. An exponential cutoff lower than 1.8 TeV is excluded at 95\% CL.

In this contribution we review the recent progress in the modeling of Pulsar Wind Nebulae (PWN). We start with a brief overview of the relevant physical processes in the magnetosphere, the wind-zone and the inflated nebula bubble. Radiative signatures and particle transport processes obtained from 3D simulations of PWN are discussed in the context of optical and X-ray observations. We then proceed to consider particle acceleration in PWN and elaborate on what can be learned about the particle acceleration from the dynamical structures called "wisps" observed in the Crab nebula. We also discuss recent observational and theoretical results of gamma-ray flares and the inner knot of the Crab nebula, which had been proposed as the emission site of the flares. We extend the discussion to GeV flares from binary systems in which the pulsar wind interacts with the stellar wind from a companion star. The chapter concludes with a discussion of solved and unsolved problems posed by PWN.

We report on the Fermi-LAT detection of high-energy emission from the behind-the-limb (BTL) solar flares that occurred on 2013 October 11, and 2014 January 6 and September 1. The Fermi-LAT observations are associated with flares from active regions originating behind both the eastern and western limbs, as determined by STEREO. All three flares are associated with very fast coronal mass ejections (CMEs) and strong solar energetic particle events. We present updated localizations of the >100 MeV photon emission, hard X-ray (HXR)and EUV images, and broadband spectra from 10 keV to 10 GeV, as well as microwave spectra. We also provide a comparison of the BTL flares detected by Fermi-LAT with three on-disk flares and present a study of some of the significant quantities of these flares as an attempt to better understand the acceleration mechanisms at work during these occulted flares. We interpret the HXR emission to be due to electron bremsstrahlung from a coronal thin-target loop top with the accelerated electron spectra steepening at semirelativistic energies. The >100 MeV gamma-rays are best described by a pion-decay model resulting from the interaction of protons (and other ions) in a thick-target photospheric source. The protons are believed to have been accelerated (to energies >10 GeV) in the CME environment and precipitate down to the photosphere from the downstream side of the CME shock and landed on the front side of the Sun, away from the original flare site and the HXR emission.

The addition of a 28 m Cherenkov telescope (CT5) to the H.E.S.S. array extended the experiment's sensitivity to lower energies. The lowest energy threshold is obtained using monoscopic analysis of data taken with CT5, providing access to gamma-ray energies below 100 GeV. Such an extension of the instrument's energy range is particularly beneficial for studies of Active Galactic Nuclei with soft spectra, as expected for those at a redshift > 0.5. The high-frequency peaked BL Lac objects PKS 2155-304 (z = 0.116) and PG 1553+113 (0.43 < z < 0.58) are among the brightest objects in the gamma-ray sky, both showing clear signatures of gamma-ray absorption at E > 100 GeV interpreted as being due to interactions with the extragalactic background light (EBL). Multiple observational campaigns of PKS 2155-304 and PG 1553+113 were conducted during 2013 and 2014 using the full H.E.S.S. II instrument. A monoscopic analysis of the data taken with the new CT5 telescope was developed along with an investigation into the systematic uncertainties on the spectral parameters. The energy spectra of PKS 2155-304 and PG 1553+113 were reconstructed down to energies of 80 GeV for PKS 2155-304, which transits near zenith, and 110 GeV for the more northern PG 1553+113. The measured spectra, well fitted in both cases by a log-parabola spectral model (with a 5.0 sigma statistical preference for non-zero curvature for PKS 2155-304 and 4.5 sigma for PG 1553+113), were found consistent with spectra derived from contemporaneous Fermi-LAT data, indicating a sharp break in the observed spectra of both sources at E ~ 100 GeV. When corrected for EBL absorption, the intrinsic H.E.S.S. II mono and Fermi-LAT spectrum of PKS 2155-304 was found to show significant curvature. For PG 1553+113, however, no significant detection of curvature in the intrinsic spectrum could be found within statistical and systematic uncertainties.

We report on Chandra X-ray Observatory (CXO) observations of the pulsar wind nebula (PWN) associated with PSR B0355+54 (eight observations with a 395 ks total exposure, performed over an 8 month period). We investigated the spatial and spectral properties of the emission coincident with the pulsar, compact nebula (CN), and extended tail. We find that the CN morphology can be interpreted in a way that suggests a small angle between the pulsar spin axis and our line-of-sight, as inferred from the radio data. On larger scales, emission from the 7' (2 pc) tail is clearly seen. We also found hints of two faint extensions nearly orthogonal to the direction of the pulsar's proper motion. The spectrum extracted at the pulsar position can be described with an absorbed power-law + blackbody model. The nonthermal component can be attributed to magnetospheric emission, while the thermal component can be attributed to emission from either a hot spot (e.g., a polar cap) or the entire neutron star surface. Surprisingly, the spectrum of the tail shows only a slight hint of cooling with increasing distance from the pulsar. This implies either a low magnetic field with fast flow speed, or particle re-acceleration within the tail. We estimate physical properties of the PWN and compare the morphologies of the CN and the extended tail with those of other bow shock PWNe observed with long CXO exposures.

The morphology of young Pulsar Wind Nebulae (PWN) is largely determined by the properties of the wind injected by the pulsar. We have used a recent parametrization of the wind obtained from Force Free Electrodynamics simulations of pulsar magnetospheres to simulate nebulae for different sets of pulsar parameters. We performed axisymmetric Relativistic Magnetohydrodynamics simulations to test the morphology dependence of the nebula on the obliquity of the pulsar and on the magnetization of the pulsar wind. We compare these simulations to the morphology of the Vela and Crab PWN. We find that the morphology of Vela can be reproduced qualitatively if the pulsar obliquity angle is alpha ~45deg and the magnetization of the wind is high (sigma_0 ~ 3.0). A morphology similar to the one of the Crab Nebula is only obtained for low magnetization simulations with alpha >~ 45deg. Interestingly, we find that Kelvin-Helmholtz instabilities produce small scale turbulences downstream of the reverse shock of the pulsar wind.

The field of gamma-ray astronomy has seen important progress during the last decade, yet there exists so far no common software framework for the scientific analysis of gamma-ray telescope data. We propose to fill this gap by means of the GammaLib software, a generic library that we have developed to support the analysis of gamma-ray event data. GammaLib has been written in C++ and all functionality is available in Python through an extension module. On top of this framework we have developed the ctools software package, a suite of software tools that enables building of flexible workflows for the analysis of Imaging Air Cherenkov Telescope event data. The ctools are inspired by science analysis software available for existing high-energy astronomy instruments, and they follow the modular ftools model developed by the High Energy Astrophysics Science Archive Research Center. The ctools have been written in Python and C++, and can be either used from the command line, via shell scripts, or directly from Python. In this paper we present the GammaLib and ctools software versions 1.0 that have been released end of 2015. GammaLib and ctools are ready for the science analysis of Imaging Air Cherenkov Telescope event data, and also support the analysis of Fermi-LAT data and the exploitation of the COMPTEL legacy data archive. We propose to use ctools as the Science Tools software for the Cherenkov Telescope Array Observatory.

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first 7 years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

Most of the celestial gamma rays detected by the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope originate from the interstellar medium when energetic cosmic rays interact with interstellar nucleons and photons. Conventional point and extended source studies rely on the modeling of this diffuse emission for accurate characterization. We describe here the development of the Galactic Interstellar Emission Model (GIEM) that is the standard adopted by the LAT Collaboration and is publicly available. The model is based on a linear combination of maps for interstellar gas column density in Galactocentric annuli and for the inverse Compton emission produced in the Galaxy. We also include in the GIEM large-scale structures like Loop I and the Fermi bubbles. The measured gas emissivity spectra confirm that the cosmic-ray proton density decreases with Galactocentric distance beyond 5 kpc from the Galactic Center. The measurements also suggest a softening of the proton spectrum with Galactocentric distance. We observe that the Fermi bubbles have boundaries with a shape similar to a catenary at latitudes below 20 degrees and we observe an enhanced emission toward their base extending in the North and South Galactic direction and located within 4 degrees of the Galactic Center.

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.

To uniformly determine the properties of supernova remnants (SNRs) at high energies, we have developed the first systematic survey at energies from 1 to 100 GeV using data from the Fermi Large Area Telescope. Based on the spatial overlap of sources detected at GeV energies with SNRs known from radio surveys, we classify 30 sources as likely GeV SNRs. We also report 14 marginal associations and 245 flux upper limits. A mock catalog in which the positions of known remnants are scrambled in Galactic longitude, allows us to determine an upper limit of 22% on the number of GeV candidates falsely identified as SNRs. We have also developed a method to estimate spectral and spatial systematic errors arising from the diffuse interstellar emission model, a key component of all Galactic Fermi LAT analyses. By studying remnants uniformly in aggregate, we measure the GeV properties common to these objects and provide a crucial context for the detailed modeling of individual SNRs. Combining our GeV results with multiwavelength (MW) data, including radio, X-ray, and TeV, demonstrates the need for improvements to previously sufficient, simple models describing the GeV and radio emission from these objects. We model the GeV and MW emission from SNRs in aggregate to constrain their maximal contribution to observed Galactic cosmic rays.

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.

Galaxy clusters are one of the prime sites to search for dark matter (DM) annihilation signals. Depending on the substructure of the DM halo of a galaxy cluster and the cross sections for DM annihilation channels, these signals might be detectable by the latest generation of $\gamma$-ray telescopes. Here we use three years of Fermi Large Area Telescope (LAT) data, which are the most suitable for searching for very extended emission in the vicinity of nearby Virgo galaxy cluster. Our analysis reveals statistically significant extended emission which can be well characterized by a uniformly emitting disk profile with a radius of 3\deg that moreover is offset from the cluster center. We demonstrate that the significance of this extended emission strongly depends on the adopted interstellar emission model (IEM) and is most likely an artifact of our incomplete description of the IEM in this region. We also search for and find new point source candidates in the region. We then derive conservative upper limits on the velocity-averaged DM pair annihilation cross section from Virgo. We take into account the potential $\gamma$-ray flux enhancement due to DM sub-halos and its complex morphology as a merging cluster. For DM annihilating into $b\overline{b}$, assuming a conservative sub-halo model setup, we find limits that are between 1 and 1.5 orders of magnitude above the expectation from the thermal cross section for $m_{\mathrm{DM}}\lesssim100\,\mathrm{GeV}$. In a more optimistic scenario, we exclude $\langle \sigma v \rangle\sim3\times10^{-26}\,\mathrm{cm^{3}\,s^{-1}}$ for $m_{\mathrm{DM}}\lesssim40\,\mathrm{GeV}$ for the same channel. Finally, we derive upper limits on the $\gamma$-ray-flux produced by hadronic cosmic-ray interactions in the inter cluster medium. We find that the volume-averaged cosmic-ray-to-thermal pressure ratio is less than $\sim6\%$.

We search for evidence of dark matter (DM) annihilation in the isotropic gamma-ray background (IGRB) measured with 50 months of Fermi Large Area Telescope (LAT) observations. An improved theoretical description of the cosmological DM annihilation signal, based on two complementary techniques and assuming generic weakly interacting massive particle (WIMP) properties, renders more precise predictions compared to previous work. More specifically, we estimate the cosmologically-induced gamma-ray intensity to have an uncertainty of a factor ~20 in canonical setups. We consistently include both the Galactic and extragalactic signals under the same theoretical framework, and study the impact of the former on the IGRB spectrum derivation. We find no evidence for a DM signal and we set limits on the DM-induced isotropic gamma-ray signal. Our limits are competitive for DM particle masses up to tens of TeV and, indeed, are the strongest limits derived from Fermi LAT data at TeV energies. This is possible thanks to the new Fermi LAT IGRB measurement, which now extends up to an energy of 820 GeV. We quantify uncertainties in detail and show the potential this type of search offers for testing the WIMP paradigm with a complementary and truly cosmological probe of DM particle signals.

Gamma-ray observations give us a direct view into the most extreme environments of the universe. They help us to study astronomical particle accelerators as supernovae remnants, pulsars, active galaxies or gamma-ray bursts and help us to understand the propagation of cosmic rays through our Milky Way. This article summarizes the status of gamma-ray observations from space; it is the write-up of a rapporteur talk given at the 34th ICRC in The Hague, The Netherlands. The primary instrument used in the presented studies is the Large Area Telescope on-board the Fermi Spacecraft, which images the whole gamma-ray sky at photon energies between 20 MeV and 2 TeV. The Fermi mission is currently in its 8th year of observations. This article will review many of the exciting discoveries made in this time, focusing on the most recent ones.

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.

Supernovae (SNe) exploding in a dense circumstellar medium (CSM) are hypothesized to accelerate cosmic rays in collisionless shocks and emit GeV gamma rays and TeV neutrinos on a time scale of several months. We perform the first systematic search for gamma-ray emission in Fermi LAT data in the energy range from 100 MeV to 300 GeV from the ensemble of 147 SNe Type IIn exploding in dense CSM. We search for a gamma-ray excess at each SNe location in a one year time window. In order to enhance a possible weak signal, we simultaneously study the closest and optically brightest sources of our sample in a joint-likelihood analysis in three different time windows (1 year, 6 months and 3 months). For the most promising source of the sample, SN 2010jl (PTF10aaxf), we repeat the analysis with an extended time window lasting 4.5 years. We do not find a significant excess in gamma rays for any individual source nor for the combined sources and provide model-independent flux upper limits for both cases. In addition, we derive limits on the gamma-ray luminosity and the ratio of gamma-ray-to-optical luminosity ratio as a function of the index of the proton injection spectrum assuming a generic gamma-ray production model. Furthermore, we present detailed flux predictions based on multi-wavelength observations and the corresponding flux upper limit at 95% confidence level (CL) for the source SN 2010jl (PTF10aaxf).

One of the most intriguing results from the gamma-ray instruments in orbit has been the detection of powerful flares from the Crab Nebula. These flares challenge our understanding of pulsar wind nebulae and models for particle acceleration. We report on the portion of a multiwavelength campaign using Keck, HST, and Chandra concentrating on a small emitting region, the Crab's inner knot, located a fraction of an arcsecond from the pulsar. We find that the knot's radial size, tangential size, peak flux, and the ratio of the flux to that of the pulsar are correlated with the projected distance of the knot from the pulsar. A new approach, using singular value decomposition for analyzing time series of images, was introduced yielding results consistent with the more traditional methods while some uncertainties were substantially reduced. We exploit the characterization of the knot to discuss constraints on standard shock-model parameters that may be inferred from our observations assuming the inner knot lies near to the shocked surface. These include inferences as to wind magnetization, shock shape parameters such as incident angle and poloidal radius of curvature, as well as the IR/optical emitting particle enthalpy fraction. We find that while the standard shock model gives good agreement with observation in many respects, there remain two puzzles: (a) The observed angular size of the knot relative to the pulsar--knot separation is much smaller than expected; (b) The variable, yet high degree of polarization reported is difficult to reconcile with a highly relativistic downstream flow.

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.

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.

We present new Jansky Very Large Array (VLA) radio images of the Crab Nebula at 5.5 GHz, taken at two epochs separated by 6 days about two months after a gamma-ray flare in 2012 July. We find no significant change in the Crab's radio emission localized to a region of <2 light-months in radius, either over the 6-day interval between our present observations or between the present observations and ones from 2001. Any radio counterpart to the flare has a radio luminosity of <~ $2 \times 10^{-4}$ times that of the nebula. Comparing our images to one from 2001, we do however find changes in radio brightness, up to 10% in amplitude, which occur on decade timescales throughout the nebula. The morphology of the changes is complex suggesting both filamentary and knotty structures. The variability is stronger, and the timescales likely somewhat shorter, nearer the centre of the nebula. We further find that even with the excellent uv~coverage and signal-to-noise of the VLA, deconvolution errors are much larger than the noise, being up to 1.2% of peak brightness of the nebula in this particular case.

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

The Crab nebula and its pulsar (referred to together as "Crab") have historically played a central role in astrophysics. True to their legacy, several unique discoveries have been made recently. The Crab was found to emit gamma-ray pulsations up to energies of 400 GeV, beyond what was previously expected from pulsars. Strong gamma-ray flares, of durations of a few days were discovered from within the nebula, while the source was previously expected to be stable in flux on these time scales. Here we review these intriguing and suggestive developments. In this context we give an overview of the observational properties of the Crab and our current understanding of pulsars and their nebulae.

Context. It has been suggested that the bow shocks of runaway stars are sources of high-energy gamma rays (E > 100 MeV). Theoretical models predicting high-energy gamma-ray emission from these sources were followed by the first detection of non-thermal radio emission from the bow shock of BD+43$^\deg$ 3654 and non-thermal X-ray emission from the bow shock of AE Aurigae. Aims. We perform the first systematic search for MeV and GeV emission from 27 bow shocks of runaway stars using data collected by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (Fermi). Methods. We analysed 57 months of Fermi-LAT data at the positions of 27 bow shocks of runaway stars extracted from the Extensive stellar BOw Shock Survey catalogue (E-BOSS). A likelihood analysis was performed to search for gamma-ray emission that is not compatible with diffuse background or emission from neighbouring sources and that could be associated with the bow shocks. Results. None of the bow shock candidates is detected significantly in the Fermi-LAT energy range. We therefore present upper limits on the high-energy emission in the energy range from 100 MeV to 300 GeV for 27 bow shocks of runaway stars in four energy bands. For the three cases where models of the high-energy emission are published we compare our upper limits to the modelled spectra. Our limits exclude the model predictions for Zeta Ophiuchi by a factor $\approx$ 5.

Current theories predict relativistic hadronic particle populations in clusters of galaxies in addition to the already observed relativistic leptons. In these scenarios hadronic interactions give rise to neutral pions which decay into $\gamma$ rays, that are potentially observable with the Large Area Telescope (LAT) on board the Fermi space telescope. We present a joint likelihood analysis searching for spatially extended $\gamma$-ray emission at the locations of 50 galaxy clusters in 4 years of Fermi-LAT data under the assumption of the universal cosmic-ray model proposed by Pinzke & Pfrommer (2010). We find an excess at a significance of $2.7\sigma$ which upon closer inspection is however correlated to individual excess emission towards three galaxy clusters: Abell 400, Abell 1367 and Abell 3112. We discuss these cases in detail and conservatively attribute the emission to unmodeled background (for example, radio galaxies within the clusters). Through the combined analysis of 50 clusters we exclude hadronic injection efficiencies in simple hadronic models above 21% and establish limits on the cosmic-ray to thermal pressure ratio within the virial radius, $R_{200}$, to be below 1.2-1.4% depending on the morphological classification. In addition we derive new limits on the $\gamma$-ray flux from individual clusters in our sample.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly-discovered features may offer a clue to the origin of high-energy CRs. We use the ${\it Fermi}$ Large Area Telescope observations of the $\gamma$-ray emission from the Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range $\sim 90~$GeV-$6~$TeV (derived from a photon energy range $15~$GeV-$1~$TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index $2.68 \pm 0.04$ and $2.61 \pm 0.08$ above $\sim 200~$GeV, respectively.

We present deep VERITAS observations of the blazar PKS 1424+240, along with contemporaneous Fermi Large Area Telescope, Swift X-ray Telescope and Swift UV Optical Telescope data between 2009 February 19 and 2013 June 8. This blazar resides at a redshift of $z\ge0.6035$, displaying a significantly attenuated gamma-ray flux above 100 GeV due to photon absorption via pair-production with the extragalactic background light. We present more than 100 hours of VERITAS observations from three years, a multiwavelength light curve and the contemporaneous spectral energy distributions. The source shows a higher flux of (2.1$\pm0.3$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV in 2009 and 2011 as compared to the flux measured in 2013, corresponding to (1.02$\pm0.08$)$\times10^{-7}$ ph m$^{-2}$s$^{-1}$ above 120 GeV. The measured differential very high energy (VHE; $E\ge100$ GeV) spectral indices are $\Gamma=$3.8$\pm$0.3, 4.3$\pm$0.6 and 4.5$\pm$0.2 in 2009, 2011 and 2013, respectively. No significant spectral change across the observation epochs is detected. We find no evidence for variability at gamma-ray opacities of greater than $\tau=2$, where it is postulated that any variability would be small and occur on longer than year timescales if hadronic cosmic-ray interactions with extragalactic photon fields provide a secondary VHE photon flux. The data cannot rule out such variability due to low statistics.

The dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

We report on a bright flare in the Crab Nebula detected by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. The period of significantly increased luminosity occurred in 2013 March and lasted for approximately 2 weeks. During this period, we observed flux variability on timescales of approximately 5\,hours. The combined photon flux above 100 MeV from the pulsar and its nebula reached a peak value of $(12.5\pm 0.8)\cdot 10^{-6}$\,cm$^{-2}$\,s$^{-1}$ on 2013 March 6. This value exceeds the average flux by almost a factor of 6 and implies a $\sim20$ times higher flux for the synchrotron component of the nebula alone. This is the second brightest flare observed from this source. Spectral and temporal analysis of the LAT data collected during the outburst reveal a rapidly varying synchrotron component of the Crab Nebula while the pulsar emission remains constant in time.

Pulsar wind nebulae (PWNe) have been established as the most populous class of TeV gamma-ray emitters. Since launch, the Fermi Large Area Telescope (LAT)identified five high-energy (100MeV <E< 100 GeV) gamma-ray sources as PWNe, and detected a large number of PWNe candidates, all powered by young and energetic pulsars. The wealth of multi-wavelength data available and the new results provided by Fermi-LAT give us an opportunity to find new PWNe and to explore the radiative processes taking place in known ones. The TeV gamma-ray unidentifiedsources (UNIDs) are the best candidates for finding new PWNe. Using 45 months of Fermi-LAT data for energies above 10 GeV, an analysis was performed near the position of 58TeV PWNe and UNIDs within 5deg of the Galactic Plane to establish new constraints on PWNe properties and find new clues on the nature of UNIDs. Of the 58 sources, 30 were detected, and this work provides their gamma-rayfluxes for energies above 10 GeV. The spectral energy distributions (SED) andupper limits, in the multi-wavelength context, also provide new information on the source nature and can help distinguish between emission scenarios, i.e. between classification as a pulsar candidate or as a PWN candidate. Six new GeV PWNe candidates are described in detail and compared with existing models. A population study of GeV PWNe candidates as a function of the pulsar/PWN system characteristics is presented.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

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.

We present results from our analysis of Chandra X-ray Observatory, W. M. Keck Observatory, and Karl G. Jansky Very Large Array (VLA) images of the Crab Nebula that were contemporaneous with the gamma-ray flare of 2011 April. Despite hints in the X-ray data, we find no evidence for statistically significant variations that pinpoint the specific location of the flares within the Nebula. The Keck observations extend this conclusion to the "inner knot", i.e., the feature within an arcsecond of the pulsar. The VLA observations support this conclusion. We also discuss theoretical implications of the gamma-ray flares and suggest that the most dramatic gamma-ray flares are due to radiation-reaction-limited synchrotron emission associated with sudden, dissipative changes in the current system sustained by the central pulsar.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

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.

We report on the gamma-ray observations of giant molecular clouds Orion A and B with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. The gamma-ray emission in the energy band between \sim100 MeV and \sim100 GeV is predicted to trace the gas mass distribution in the clouds through nuclear interactions between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production cross-section for the nuclear interaction is known to \sim10% precision which makes the LAT a powerful tool to measure the gas mass column density distribution of molecular clouds for a known CR intensity. We present here such distributions for Orion A and B, and correlate them with those of the velocity integrated CO intensity (WCO) at a 1{\deg} \times1{\deg} pixel level. The correlation is found to be linear over a WCO range of ~10 fold when divided in 3 regions, suggesting penetration of nuclear CRs to most of the cloud volumes. The Wco-to-mass conversion factor, Xco, is found to be \sim2.3\times10^20 cm-2(K km s-1)-1 for the high-longitude part of Orion A (l > 212{\deg}), \sim1.7 times higher than \sim1.3 \times 10^20 found for the rest of Orion A and B. We interpret the apparent high Xco in the high-longitude region of Orion A in the light of recent works proposing a non-linear relation between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than the H2 column density in the region making the gas "darker" to Wco.

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.

We present time-resolved broad-band observations of the quasar 3C 279 obtained from multi-wavelength campaigns conducted during the first two years of the Fermi Gamma-ray Space Telescope mission. While investigating the previously reported gamma-ray/optical flare accompanied by a change in optical polarization, we found that the optical emission appears delayed with respect to the gamma-ray emission by about 10 days. X-ray observations reveal a pair of `isolated' flares separated by ~90 days, with only weak gamma-ray/optical counterparts. The spectral structure measured by Spitzer reveals a synchrotron component peaking in the mid-infrared band with a sharp break at the far-infrared band during the gamma-ray flare, while the peak appears in the mm/sub-mm band in the low state. Selected spectral energy distributions are fitted with leptonic models including Comptonization of external radiation produced in a dusty torus or the broad-line region. Adopting the interpretation of the polarization swing involving propagation of the emitting region along a curved trajectory, we can explain the evolution of the broad-band spectra during the gamma-ray flaring event by a shift of its location from ~ 1 pc to ~ 4 pc from the central black hole. On the other hand, if the gamma-ray flare is generated instead at sub-pc distance from the central black hole, the far-infrared break can be explained by synchrotron self-absorption. We also model the low spectral state, dominated by the mm/sub-mm peaking synchrotron component, and suggest that the corresponding inverse-Compton component explains the steady X-ray emission.

Dark matter particle annihilation or decay can produce monochromatic gamma-ray lines and contribute to the diffuse gamma-ray background. Flux upper limits are presented for gamma-ray spectral lines from 7 to 200 GeV and for the diffuse gamma-ray background from 4.8 GeV to 264 GeV obtained from two years of Fermi Large Area Telescope data integrated over most of the sky. We give cross section upper limits and decay lifetime lower limits for dark matter models that produce gamma-ray lines or contribute to the diffuse spectrum, including models proposed as explanations of the PAMELA and Fermi cosmic-ray data.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

The Large Area Telescope on board the Fermi satellite observed a gamma-ray flare in the Crab nebula lasting for approximately nine days in April of 2011. The source, which at optical wavelengths has a size of ~11 ly across, doubled its gamma-ray flux within eight hours. The peak photon flux was (186 +- 6) 10-7 cm-2 s-1 above 100 MeV, which corresponds to a 30-fold increase compared to the average value. During the flare, a new component emerged in the spectral energy distribution, which peaked at an energy of (375 +- 26) MeV at flare maximum. The observations imply that the emission region was likely relativistically beamed toward us and that variations in its motion are responsible for the observed spectral variability.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

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.

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.

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.