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.

We report on the timing observations of the millisecond pulsar PSR J2055+3829 originally discovered as part of the SPAN512 survey conducted with the Nan\c{c}ay Radio Telescope. The pulsar has a rotational period of 2.089 ms, and is in a tight 3.1 hr orbit around a very low mass ($0.023 \leq m_c \lesssim 0.053$ M$_\odot$, 90\% c.l.) companion. Our 1.4 GHz observations reveal the presence of eclipses of the pulsar's radio signal caused by the outflow of material from the companion, for a few minutes around superior conjunction of the pulsar. The very low companion mass, the observation of radio eclipses, and the detection of time variations of the orbital period establish PSR J2055+3829 as a `black widow' (BW) pulsar. Inspection of the radio signal from the pulsar during ingress and egress phases shows that the eclipses in PSR J2055+3829 are asymmetric and variable, as is commonly observed in other similar systems. More generally, the orbital properties of the new pulsar are found to be very similar to those of other known eclipsing BW pulsars. No gamma-ray source is detected at the location of the pulsar in recent \textit{Fermi}-LAT source catalogs. We used the timing ephemeris to search ten years of \textit{Fermi} Large Area Telescope (LAT) data for gamma-ray pulsations, but were unable to detect any, possibly because of the pulsar's large distance. We finally compared the mass functions of eclipsing and non-eclipsing BW pulsars and confirmed previous findings that eclipsing BWs have higher mass functions than their non-eclipsing counterparts. Larger inclinations could explain the higher mass functions of eclipsing BWs. On the other hand, the mass function distributions of Galactic disk and globular cluster BWs appear to be consistent, suggesting, despite the very different environments, the existence of common mechanisms taking place in the last stages of evolution of BWs.

The current version of the Fermi Large Area Telescope data (P8R2) has been publicly available since June 2015, with the caveat that the residual background of all event classes, except ULTRACLEANVETO, was not fully isotropic: it was enhanced by a factor ~2 at 1-3 GeV within ~20 deg of the Ecliptic compared to the poles. By investigating the residual background using data only, we were able to find two sources of residual background: one due to non-interacting heavy ions and one due to cosmic-ray electrons leaking through the ribbons of the Anti-Coincidence Detector, the latter source being responsible for the background anisotropy. A set of simple cuts allows us to reject these events while losing less than 1% of the SOURCE class acceptance. This new selection has been used to produce a new version of the LAT data (P8R3).

PMN J1603$-$4904 is a likely member of the rare class of $\gamma$-ray emitting young radio galaxies. Only one other source, PKS 1718$-$649, has been confirmed so far. These objects, which may transition into larger radio galaxies, are a stepping stone to understanding AGN evolution. It is not completely clear how these young galaxies, seen edge-on, can produce high-energy $\gamma$-rays. PMN J1603$-$4904 has been detected by TANAMI Very Long Baseline Interferometry (VLBI) observations and has been followed-up with multiwavelength observations. A Fermi/LAT $\gamma$-ray source has been associated with it in the LAT catalogs. We have obtained Chandra observations of the source in order to consider the possibility of source confusion, due to the relatively large positional uncertainty of Fermi/LAT. The goal was to investigate the possibility of other X-ray bright sources in the vicinity of PMN J1603$-$4904 that could be counterparts to the $\gamma$-ray emission. With Chandra/ACIS, we find no other sources in the uncertainty ellipse of Fermi/LAT data, which includes an improved localization analysis of 8 years of data. We further study the X-ray fluxes and spectra. We conclude that PMN J1603$-$4904 is indeed the second confirmed $\gamma$-ray bright young radio galaxy.

We describe an image-based method that uses two radio criteria, compactness and spectral index, to identify promising pulsar candidates among Fermi Large Area Telescope (LAT) unassociated sources. These criteria are applied to those radio sources from the Giant Metrewave Radio Telescope all-sky survey at 150 MHz (TGSS ADR1) found within the error ellipses of unassociated sources from the 3FGL catalog and a preliminary source list based on 7 years of LAT data. After follow-up interferometric observations to identify extended or variable sources, a list of 16 compact, steep-spectrum candidates is generated. An ongoing search for pulsations in these candidates, in gamma rays and radio, has found six millisecond pulsars and one normal pulsar. A comparison of this method with existing selection criteria based on gamma-ray spectral and variability properties suggests that the pulsar discovery space using Fermi may be larger than previously thought. Radio imaging is a hitherto underutilized source selection method that can be used, as with other multi-wavelength techniques, in the search for Fermi pulsars.

We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $\gamma$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$\sigma$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $\alpha\sim-3$ (where $S \propto \nu^{\alpha}$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607.

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.

Blazars are classically divided into the BL Lac (BLL) and Flat-Spectrum Radio Quasar (FSRQ) subclasses, corresponding to radiatively inefficient and efficient accretion regimes, respectively, largely based on the equivalent width (EW) of their optical broad emission lines (BEL). However, EW-based classification criteria are not physically motivated, and a few blazars have previously 'transitioned' from one subclass to the other. We present the first systematic search for these transition blazars in a sample of 602 unique pairs of repeat spectra of 354 blazars in SDSS, finding six clear cases. These transition blazars have bolometric Eddington ratios of ~0.3 and low-frequency synchrotron peaks, and are thus FSRQ-like. We show that the strong EW variability (up to an unprecedented factor of >60) is due to swamping of the BELs from variability in jet continuum emission, which is stronger in amplitude and shorter in timescale than typical blazars. Although these transition blazars appear to switch between FSRQ and BLL according to the phenomenologically-based EW scheme, we show that they are most likely rare cases of FSRQs with radiatively efficient accretion flows and especially strongly-beamed jets. These results have implications for the decrease of the apparent BLL population at high-redshifts, and may add credence to claims of a negative BLL redshift evolution.

We use Fermi Gamma-ray Space Telescope detections and upper limits on non-recycled pulsars obtained from the Large Area Telescope (LAT) to constrain how the gamma-ray luminosity L depends on the period P and the period derivative \dot{P}. We use a Bayesian analysis to calculate a best-fit luminosity law, or dependence of L on P and \dot{P}, including different methods for modeling the beaming factor. An outer gap (OG) magnetosphere geometry provides the best-fit model, which is L \propto P^{-a} \dot{P}^{b} where a=1.36\pm0.03 and b=0.44\pm0.02, similar to but not identical to the commonly assumed L \propto \sqrt{\dot{E}} \propto P^{-1.5} \dot{P}^{0.5}. Given upper limits on gamma-ray fluxes of currently known radio pulsars and using the OG model, we find that about 92% of the radio-detected pulsars have gamma-ray beams that intersect our line of sight. By modeling the misalignment of radio and gamma-ray beams of these pulsars, we find an average gamma-ray beaming solid angle of about 3.7{\pi} for the OG model, assuming a uniform beam. Using LAT-measured diffuse fluxes, we place a 2{\sigma} upper limit on the average braking index and a 2{\sigma} lower limit on the average surface magnetic field strength of the pulsar population of 3.8 and 3.2 X 10^{10} G, respectively. We then predict the number of non-recycled pulsars detectable by the LAT based on our population model. Using the two-year sensitivity, we find that the LAT is capable of detecting emission from about 380 non-recycled pulsars, including 150 currently identified radio pulsars. Using the expected five-year sensitivity, about 620 non-recycled pulsars are detectable, including about 220 currently identified radio pulsars. We note that these predictions significantly depend on our model assumptions.

We present a catalog of gamma-ray sources at energies above 10 GeV based on data from the Large Area Telescope (LAT) accumulated during the first 3 years of the Fermi Gamma-ray Space Telescope mission. This catalog complements the Second Fermi-LAT Catalog, which was based on 2 years of data extending down to 100 MeV and so included many sources with softer spectra below 10 GeV. The First Fermi-LAT Catalog of >10 GeV sources (1FHL) has 514 sources, and includes their locations, spectra, a measure of their variability, and associations with cataloged sources at other wavelengths. We found that 449 (87%) could be associated with known sources, of which 393 (76% of the 1FHL sources) are active galactic nuclei. We also highlight the subset of the 1FHL sources that are the best candidates for detection at energies above 50 GeV with ground-based gamma-ray observatories.

We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.

Spatial extension is an important characteristic for correctly associating gamma-ray-emitting sources with their counterparts at other wavelengths and for obtaining an unbiased model of their spectra. We present a new method for quantifying the spatial extension of sources detected by the Large Area Telescope (LAT), the primary science instrument on the Fermi Gamma-ray Space Telescope (Fermi). We perform a series of Monte Carlo simulations to validate this tool and calculate the LAT threshold for detecting the spatial extension of sources. We then test all sources in the second Fermi-LAT catalog (2FGL) for extension. We report the detection of seven new spatially extended sources.

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.

We report on a detailed investigation of the gamma-ray emission from 18 broad line radio galaxies (BLRGs) based on two years of Fermi Large Area Telescope (LAT) data. We confirm the previously reported detections of 3C 120 and 3C 111 in the GeV photon energy range; a detailed look at the temporal characteristics of the observed gamma-ray emission reveals in addition possible flux variability in both sources. No statistically significant gamma-ray detection of the other BLRGs was however found in the considered dataset. Though the sample size studied is small, what appears to differentiate 3C 111 and 3C 120 from the BLRGs not yet detected in gamma-rays is the particularly strong nuclear radio flux. This finding, together with the indications of the gamma-ray flux variability and a number of other arguments presented, indicate that the GeV emission of BLRGs is most likely dominated by the beamed radiation of relativistic jets observed at intermediate viewing angles. In this paper we also analyzed a comparison sample of high accretion-rate Seyfert 1 galaxies, which can be considered radio-quiet counterparts of BLRGs, and found none were detected in gamma-rays. A simple phenomenological hybrid model applied for the broad-band emission of the discussed radio-loud and radio-quiet type 1 active galaxies suggests that the relative contribution of the nuclear jets to the accreting matter is > 1 percent on average for BLRGs, whilst <0.1 percent for Seyfert 1 galaxies.

For more than 15 years, since the days of the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (CGRO; 1991-2000), it has remained an open question why the prominent blazar 3C 345 was not reliably detected at gamma-ray energies <=20 MeV. Recently a bright gamma-ray source (0FGL J1641.4+3939/1FGL J1642.5+3947), potentially associated with 3C 345, was detected by the Large Area Telescope (LAT) on Fermi. Multiwavelength observations from radio bands to X-rays (mainly GASP-WEBT and Swift) of possible counterparts (3C 345, NRAO 512, B3 1640+396) were combined with 20 months of Fermi-LAT monitoring data (August 2008 - April 2010) to associate and identify the dominating gamma-ray emitting counterpart of 1FGL J1642.5+3947. The source 3C 345 is identified as the main contributor for this gamma-ray emitting region. However, after November 2009 (15 months), a significant excess of photons from the nearby quasar NRAO 512 started to contribute and thereafter was detected with increasing gamma-ray activity, possibly adding flux to 1FGL J1642.5+3947. For the same time period and during the summer of 2010, an increase of radio, optical and X-ray activity of NRAO 512 was observed. No gamma-ray emission from B3 1640+396 was detected.

We report on a newly detected point-like source, HESS J1943+213 located in the Galactic plane. This source coincides with an unidentified hard X-ray source IGR J19443+2117, which was proposed to have radio and infrared counterparts. HESS J1943+213 is detected at the significance level of 7.9 \sigma (post-trials) at RA(J2000)=19h 43m 55s +- 1s (stat) +- 1s (sys), DEC(J2000) = +21deg 18' 8" +- 17" (stat) +- 20" (sys). The source has a soft spectrum with photon index Gamma = 3.1 +- 0.3 (stat) +- 0.2 (sys) and a flux above 470 GeV of 1.3 +- 0.2 (stat) +- 0.3 (sys) x 10^{-12} cm^{-2} s^{-1}. There is no Fermi/LAT counterpart down to a flux limit of 6 x 10^{-9} cm^{-2} s^{-1} in the 0.1-100 GeV energy range (95% confidence upper limit calculated for an assumed power-law model with a photon index Gamma=2.0). The data from radio to VHE gamma-rays do not show any significant variability. We combine new H.E.S.S., Fermi/LAT and Nancay Radio Telescope observations with pre-existing non-simultaneous multi-wavelength observations of IGR J19443+2117 and discuss the likely source associations as well as the interpretation as an active galactic nucleus, a gamma-ray binary or a pulsar wind nebula. The lack of a massive stellar counterpart disfavors the binary hypothesis, while the soft VHE spectrum would be very unusual in case of a pulsar wind nebula. In addition, the distance estimates for Galactic counterparts places them outside of the Milky Way. All available observations favor an interpretation as an extreme, high-frequency peaked BL Lac object with a redshift z>0.14. This would be the first time a blazar is detected serendipitously from ground-based VHE observations, and the first VHE AGN detected in the Galactic Plane.

Most extragalactic jets in radio-loud Active Galactic Nuclei are bright and variable gamma-ray sources, which are continuously monitored with Fermi/LAT. We present the gamma-ray properties of the MOJAVE and TANAMI AGN samples of radio-loud AGN. Both programs provide properties of the parsec-scale radio jets using Very Long Baseline Interferometry (VLBI) techniques. This information is important to understand the broad-band emission mechanism of these sources. In this work we compare the radio and gamma -ray properties of the two samples and present upper limits on the gamma-ray flux of the radio-brightest jet sources not yet detected by Fermi/LAT.

We report the discovery of eight gamma-ray pulsars in blind frequency searches using the LAT, onboard the Fermi Gamma-ray Space Telescope. Five of the eight pulsars are young (tau_c<100 kyr), energetic (Edot>10^36 erg/s), and located within the Galactic plane (|b|<3 deg). The remaining three are older, less energetic, and located off the plane. Five pulsars are associated with sources included in the LAT bright gamma-ray source list, but only one, PSR J1413-6205, is clearly associated with an EGRET source. PSR J1023-5746 has the smallest characteristic age (tau_c=4.6 kyr) and is the most energetic (Edot=1.1E37 erg/s) of all gamma-ray pulsars discovered so far in blind searches. PSRs J1957+5033 and J2055+25 have the largest characteristic ages (tau_c~1 Myr) and are the least energetic (Edot~5E33 erg/s) of the newly-discovered pulsars. We present the timing models, light curves, and detailed spectral parameters of the new pulsars. We used recent XMM observations to identify the counterpart of PSR J2055+25 as XMMU J205549.4+253959. In addition, publicly available archival Chandra X-ray data allowed us to identify the likely counterpart of PSR J1023-5746 as a faint, highly absorbed source, CXOU J102302.8-574606. The large X-ray absorption indicates that this could be among the most distant gamma-ray pulsars detected so far. PSR J1023-5746 is positionally coincident with the TeV source HESS J1023-575, located near the young stellar cluster Westerlund 2, while PSR J1954+2836 is coincident with a 4.3 sigma excess reported by Milagro at a median energy of 35 TeV. Deep radio follow-up observations of the eight pulsars resulted in no detections of pulsations and upper limits comparable to the faintest known radio pulsars, indicating that these can be included among the growing population of radio-quiet pulsars in our Galaxy being uncovered by the LAT, and currently numbering more than 20.

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)

The TANAMI program has been monitoring the parsec-scale radio jets of southern gamma-ray bright AGN with VLBI techniques simultaneously with Fermi/LAT monitoring of their gamma-ray emission. Here we present the gamma-ray properties of the TANAMI sources based on an analysis of the preliminary 1-year LAT source list. We present upper limits on the gamma-ray flux for TANAMI sources not detected by LAT.

We present results of a study of the localization capability of Fermi-LAT, using a large set of blazars with precise radio locations. Since the width of the PSF decreases with energy, the performance is typically dominated by a few high energy photons, so it is important to properly characterize the high-energy PSF. Using such data, we have found a need to modify the pre-launch high-energy (greater than a few GeV) PSF derived from extensive Monte Carlo simulations of particle interactions in the LAT; the resulting data-based PSF is shown

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

Pulsars are among the prime targets for the Large Area Telescope (LAT) aboard the recently launched Fermi observatory. The LAT will study the gamma-ray Universe between 20 MeV and 300 GeV with unprecedented detail. Increasing numbers of gamma-ray pulsars are being firmly identified, yet their emission mechanisms are far from being understood. To better investigate and exploit the LAT capabilities for pulsar science, a set of new detailed pulsar simulation tools have been developed within the LAT collaboration. The structure of the pulsar simulator package PulsarSpectrum is presented here. Starting from photon distributions in energy and phase obtained from theoretical calculations or phenomenological considerations, gamma rays are generated and their arrival times at the spacecraft are determined by taking into account effects such as barycentric effects and timing noise. Pulsars in binary systems also can be simulated given orbital parameters. We present how simulations can be used for generating a realistic set of gamma rays as observed by the LAT, focusing on some case studies that show the performance of the LAT for pulsar observations.

This paper presents the simulation of the GLAST high energy gamma-ray telescope. The simulation package, written in C++, is based on the Geant4 toolkit, and it is integrated into a general framework used to process events. A detailed simulation of the electronic signals inside Silicon detectors has been provided and it is used for the particle tracking, which is handled by a dedicated software. A unique repository for the geometrical description of the detector has been realized using the XML language and a C++ library to access this information has been designed and implemented.

The GLAST Balloon Flight Engineering Model (BFEM) represents one of 16 towers that constitute the Large Area Telescope (LAT), a high-energy (>20 MeV) gamma-ray pair-production telescope being built by an international partnership of astrophysicists and particle physicists for a satellite launch in 2006. The prototype tower consists of a Pb/Si pair-conversion tracker (TKR), a CsI hodoscopic calorimeter (CAL), an anti-coincidence detector (ACD) and an autonomous data acquisition system (DAQ). The self-triggering capabilities and performance of the detector elements have been previously characterized using positron, photon and hadron beams. External target scintillators were placed above the instrument to act as sources of hadronic showers. This paper provides a comprehensive description of the BFEM data-reduction process, from receipt of the flight data from telemetry through event reconstruction and background rejection cuts. The goals of the ground analysis presented here are to verify the functioning of the instrument and to validate the reconstruction software and the background-rejection scheme.

This paper reviews the field of gamma-ray astronomy and describes future experiments and prospects for advances in fundamental physics and high-energy astrophysics through gamma-ray measurements. We concentrate on recent progress in the understanding of active galaxies, and the use of these sources as probes of intergalactic space. We also describe prospects for future experiments in a number of areas of fundamental physics, including: searches for an annihilation line from neutralino dark matter, understanding the energetics of supermassive black holes, using AGNs as cosmological probes of the primordial radiation fields, constraints on quantum gravity, detection of a new spectral component from GRBs, and the prospects for detecting primordial black holes.