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.

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ~5.0 x 10^14 g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.

This paper reports the results from three targeted searches of Milagro TeV sky maps: two extragalactic point source lists and one pulsar source list. The first extragalactic candidate list consists of 709 candidates selected from the Fermi-LAT 2FGL catalog. The second extragalactic candidate list contains 31 candidates selected from the TeVCat source catalog that have been detected by imaging atmospheric Cherenkov telescopes (IACTs). In both extragalactic candidate lists Mkn 421 was the only source detected by Milagro. This paper presents the Milagro TeV flux for Mkn 421 and flux limits for the brighter Fermi-LAT extragalactic sources and for all TeVCat candidates. The pulsar list extends a previously published Milagro targeted search for Galactic sources. With the 32 new gamma-ray pulsars identified in 2FGL, the number of pulsars that are studied by both Fermi-LAT and Milagro is increased to 52. In this sample, we find that the probability of Milagro detecting a TeV emission coincident with a pulsar increases with the GeV flux observed by the Fermi-LAT in the energy range from 0.1 GeV to 100 GeV.

TeV flaring activity with time scales as short as tens of minutes and an orphan TeV flare have been observed from the blazar Markarian 421 (Mrk 421). The TeV emission from Mrk 421 is believed to be produced by leptonic synchrotron self-Compton (SSC) emission. In this scenario, correlations between the X-ray and the TeV fluxes are expected, TeV orphan flares are hardly explained and the activity (measured as duty cycle) of the source at TeV energies is expected to be equal or less than that observed in X-rays if only SSC is considered. To estimate the TeV duty cycle of Mrk 421 and to establish limits on its variability at different time scales, we continuously observed Mrk 421 with the Milagro observatory. Mrk 421 was detected by Milagro with a statistical significance of 7.1 standard deviations between 2005 September 21 and 2008 March 15. The observed spectrum is consistent with previous observations by VERITAS. We estimate the duty cycle of Mrk 421 for energies above 1 TeV for different hypothesis of the baseline flux and for different flare selections and we compare our results with the X-ray duty cycle estimated by Resconi et al. 2009. The robustness of the results is discussed.

We report on broad multi-wavelength observations of the 2010-2011 periastron passage of the gamma-ray loud binary system PSR B1259-63. High resolution interferometric radio observations establish extended radio emission trailing the position of the pulsar. Observations with the Fermi Gamma-ray Space Telescope reveal GeV gamma-ray flaring activity of the system, reaching the spin-down luminosity of the pulsar, around 30 days after periastron. There are no clear signatures of variability at radio, X-ray and TeV energies at the time of the GeV flare. Variability around periastron in the H$\alpha$ emission line, can be interpreted as the gravitational interaction between the pulsar and the circumstellar disk. The equivalent width of the H$\alpha$ grows from a few days before periastron until a few days later, and decreases again between 18 and 46 days after periastron. In near infrared we observe the similar decrease of the equivalent width of Br$\gamma$ line between the 40th and 117th day after the periastron. For the idealized disk, the variability of the H$\alpha$ line represents the variability of the mass and size of the disk. We discuss possible physical relations between the state of the disk and GeV emission under assumption that GeV flare is directly related to the decrease of the disk size.

The very high energy (VHE; >100 GeV) source HESS J0632+057 has been recently confirmed as a \gamma-ray binary, a subclass of the high mass X-ray binary (HMXB) population, through the detection of an orbital period of 321 days. We performed a deep search for the emission of HESS J0632+057 in the GeV energy range using data from the Fermi Large Area Telescope (LAT). The analysis was challenging due to the source being located in close proximity to the bright \gamma-ray pulsar PSR J0633+0632 and lying in a crowded region of the Galactic plane where there is prominent diffuse emission. We formulated a Bayesian block algorithm adapted to work with weighted photon counts, in order to define the off-pulse phases of PSR J0633+0632. A detailed spectral-spatial model of a 5 deg circular region centred on the known location of HESS J0632+057 was generated to accurately model the LAT data. No significant emission from the location of HESS J0632+057 was detected in the 0.1-100 GeV energy range integrating over ~3.5 years of data; with a 95% flux upper limit of F_{0.1-100 GeV} < 3 x 10-8 ph cm-2 s-1. A search for emission over different phases of the orbit also yielded no significant detection. A search for source emission on shorter timescales (days--months) did not yield any significant detections. We also report the results of a search for radio pulsations using the 100-m Green Bank Telescope (GBT). No periodic signals or individual dispersed bursts of a likely astronomical origin were detected. We estimated the flux density limit of < 90/40 \mu Jy at 2/9 GHz. The LAT flux upper limits combined with the detection of HESS J0632+057 in the 136-400 TeV energy band by the MAGIC collaboration imply that the VHE spectrum must turn over at energies <136 GeV placing constraints on any theoretical models invoked to explain the \gamma-ray emission.

We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1 to 200 TeV using the last 3 years of Milagro data (2005-2008), with the detector in its final configuration. MGRO J2019+37 is detected with a significance of 12.3 standard deviations ($\sigma$), and is better fit by a power law with an exponential cutoff than by a simple power law, with a probability $>98$% (F-test). The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of $7^{+5}_{-2}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$, a spectral index of $2.0^{+0.5}_{-1.0}$ and a cutoff energy of $29^{+50}_{-16}$ TeV. MGRO J2031+41 is detected with a significance of 7.3$\sigma$, with no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of $2.4^{+0.6}_{-0.5}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$ and a spectral index of $3.08^{+0.19}_{-0.17}$. The overall flux is subject to an $\sim$30% systematic uncertainty. The systematic uncertainty on the power law indices is $\sim$0.1. A comparison with previous results from TeV J2032+4130, MGRO J2031+41 and MGRO J2019+37 is also presented.

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.

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV - 100 TeV gamma ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between ~1 and ~100 TeV. The TeV emission is believed to be caused by inverse-Compton up-scattering scattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fit. The fit is consistent with values measured by IACTs between 1 and 20 TeV. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4 to 2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20 to 40 TeV.

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

One of the main results of the Fermi Gamma-Ray Space Telescope is the discovery of {\gamma}-ray selected pulsars. The high magnetic field pulsar, PSR J0007+7303 in CTA1, was the first ever to be discovered through its {\gamma}-ray pulsations. Based on analysis of 2 years of LAT survey data, we report on the discovery of {\gamma}-ray emission in the off-pulse phase interval at the ~ 6{\sigma} level. The flux from this emission in the energy range E \geq 100 MeV is F_100 = (1.73\pm0.40)\times10^(-8) photons/cm^2/s and is best fitted by a power law with a photon index of {\Gamma} = 2.54\pm0.14. The pulsed {\gamma}-ray flux in the same energy range is F_100 = (3.95\pm0.07)\times10^(-7) photons/cm^2/s and is best fitted by an exponentially-cutoff power-law spectrum with a photon index of {\Gamma} = 1.41 \pm 0.23 and a cutoff energy E_c = 4.04 \pm 0.20 GeV. We find no flux variability neither at the 2009 May glitch nor in the long term behavior. We model the {\gamma}-ray light curve with two high-altitude emission models, the outer gap and slot gap, and find that the model that best fits the data depends strongly on the assumed origin of the off-pulse emission. Both models favor a large angle between the magnetic axis and observer line of sight, consistent with the nondetection of radio emission being a geometrical effect. Finally we discuss how the LAT results bear on the understanding of the cooling of this neutron star.

We report on the discovery of \geq 100 MeV {\gamma} rays from the binary system PSR B1259-63/LS 2883 using the Large Area Telescope (LAT) on board Fermi. The system comprises a radio pulsar in orbit around a Be star. We report on LAT observations from near apastron to ~ 60 days after the time of periastron, tp, on 2010 December 15. No {\gamma}-ray emission was detected from this source when it was far from periastron. Faint {\gamma}-ray emission appeared as the pulsar approached periastron. At ~ tp + 30d, the \geq 100 MeV {\gamma}-ray flux increased over a period of a few days to a peak flux 20-30 times that seen during the pre-periastron period, but with a softer spectrum. For the following month, it was seen to be variable on daily time scales, but remained at ~ 1 - 4 \times 10^-6 cm^-2 s^-1 before starting to fade at ~ tp + 57d. The total {\gamma}-ray luminosity observed during this period is comparable to the spin-down power of the pulsar. Simultaneous radio and X-ray observations of the source showed no corresponding dramatic changes in radio and X-ray flux between the pre-periastron and post-periastron flares. We discuss possible explanations for the observed {\gamma}-ray-only flaring of the source.

We report the detection of high-energy gamma rays from the quiescent Sun with the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope (Fermi) during the first 18 months of the mission. These observations correspond to the recent period of low solar activity when the emission induced by cosmic rays is brightest. For the first time, the high statistical significance of the observations allows clear separation of the two components: the point-like emission from the solar disk due to cosmic ray cascades in the solar atmosphere, and extended emission from the inverse Compton scattering of cosmic ray electrons on solar photons in the heliosphere. The observed integral flux (>100 MeV) from the solar disk is (4.6 +/- 0.2 [statistical error] +1.0/-0.8 [systematic error]) x10^{-7} cm^{-2} s^{-1}, which is ~7 times higher than predicted by the "nominal" model of Seckel et al. (1991). In contrast, the observed integral flux (>100 MeV) of the extended emission from a region of 20 deg radius centered on the Sun, but excluding the disk itself, (6.8 +/-0.7 [stat.] +0.5/-0.4 [syst.]) x10^{-7} cm^{-2} s^{-1}, along with the observed spectrum and the angular profile, are in good agreement with the theoretical predictions for the inverse Compton emission.

We present observations of the young Supernova remnant (SNR) RX J1713.7-3946 with the Fermi Large Area Telescope (LAT). We clearly detect a source positionally coincident with the SNR. The source is extended with a best-fit extension of 0.55$^{\circ} \pm 0.04^{\circ}$ matching the size of the non-thermal X-ray and TeV gamma-ray emission from the remnant. The positional coincidence and the matching extended emission allows us to identify the LAT source with the supernova remnant RX J1713.7-3946. The spectrum of the source can be described by a very hard power-law with a photon index of $\Gamma = 1.5 \pm 0.1$ that coincides in normalization with the steeper H.E.S.S.-detected gamma-ray spectrum at higher energies. The broadband gamma-ray emission is consistent with a leptonic origin as the dominant mechanism for the gamma-ray emission.

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nan\c{c}ay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

Using the Parkes radio telescope we have carried out deep observations of eleven unassociated gamma-ray sources. Periodicity searches of these data have discovered two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 (1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). In addition we searched for but did not detect any radio pulsations from six gammaray pulsars discovered by the Fermi satellite to a level of - 0.04 mJy (for pulsars with a 10% duty cycle). Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 9' from the centroid of the gamma-ray source. Since these observations were carried out, independent evidence has shown that 1FGL J1103.9-5355 is associated with the flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not associated with the gamma-ray source, despite the seemingly low probability of a chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a chance discovery of a weak, long period pulsar and is unlikely to be associated with 1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very low mass companion with a 3.5 hour orbital period. The relatively high flux density and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for high precision timing experiments. The gamma-rays of 1FGL J2241.9-5236 have a spectrum that is well modelled by a power law with exponential cutoff, and phasebinning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. Observations with Chandra have identified a coincident X-ray source within 0.1" of the position of the pulsar obtained by radio timing

We report on observations of BL Lacertae during the first 18 months of Fermi-LAT science operations and present results from a 48-day multifrequency coordinated campaign from 2008 August 19 to 2008 October 7. The radio to gamma-ray behavior of BL Lac is unveiled during a low activity state thanks to the coordinated observations of radio-band (Metsahovi and VLBA), near-IR/optical (Tuorla, Steward, OAGH and MDM) and X-ray (RXTE and Swift) observatories. No variability was resolved in gamma-rays during the campaign, and the brightness level was 15 times lower than the level of the 1997 EGRET outburst. Moderate and uncorrelated variability has been detected in UV and X-rays. The X-ray spectrum is found to be concave indicating the transition region between the low and high energy component of the spectral energy distribution (SED). VLBA observation detected a synchrotron spectrum self-absorption turnover in the innermost part of the radio jet appearing to be elongated and inhomogeneous, and constrained the average magnetic field there to be less than 3 G. Over the following months BL Lac appeared variable in gamma-rays, showing flares (in 2009 April and 2010 January). There is no evidence for correlation of the gamma-rays with the optical flux monitored from the ground in 18 months. The SED may be described by a single zone or two zone synchrotron self-Compton (SSC) model, but a hybrid SSC plus external radiation Compton (ERC) model seems preferred based on the observed variability and the fact that it provides a fit closest to equipartition.

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.

The dramatic increase in the number of known gamma-ray pulsars since the launch of the Fermi Gamma-ray Space Telescope (formerly GLAST) offers the first opportunity to study a population of these high-energy objects. This catalog summarizes 46 high-confidence pulsed detections using the first six months of data taken by the Large Area Telescope (LAT), Fermi's main instrument. Sixteen previously unknown pulsars were discovered by searching for pulsed signals at the positions of bright gamma-ray sources seen with the LAT, or at the positions of objects suspected to be neutron stars based on observations at other wavelengths. Pulsed gamma-ray emission was discovered from twenty-four known pulsars by using ephemerides (timing solutions) derived from monitoring radio pulsars. Eight of these new gamma-ray pulsars are millisecond pulsars. The pulsed energy spectra can be described by a power law with an exponential cutoff, with cutoff energies in the range 1 to 5 GeV. The rotational energy loss rate (\dot{E}) of these neutron stars spans 5 decades, from ~3x10^{33} erg/s to 5x10^{38} erg/s, and the apparent efficiencies for conversion to gamma-ray emission range from ~0.1% to unity, although distance uncertainties complicate efficiency estimates. The pulse shapes show substantial diversity, but roughly 75% of the gamma-ray pulse profiles have two peaks, separated by >0.2 of rotational phase. For most of the pulsars, gamma-ray emission appears to come mainly from the outer magnetosphere, while polar-cap emission remains plausible for a remaining few. Finally, these discoveries suggest that gamma-ray-selected young pulsars are born at a rate comparable to that of their radio-selected cousins and that the birthrate of all young gamma-ray-detected pulsars is a substantial fraction of the expected Galactic supernova rate.

We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

The Extragalactic Background Light (EBL) includes photons with wavelengths from ultraviolet to infrared, which are effective at attenuating gamma rays with energy above ~10 GeV during propagation from sources at cosmological distances. This results in a redshift- and energy-dependent attenuation of the gamma-ray flux of extragalactic sources such as blazars and Gamma-Ray Bursts (GRBs). The Large Area Telescope onboard Fermi detects a sample of gamma-ray blazars with redshift up to z~3, and GRBs with redshift up to z~4.3. Using photons above 10 GeV collected by Fermi over more than one year of observations for these sources, we investigate the effect of gamma-ray flux attenuation by the EBL. We place upper limits on the gamma-ray opacity of the Universe at various energies and redshifts, and compare this with predictions from well-known EBL models. We find that an EBL intensity in the optical-ultraviolet wavelengths as great as predicted by the "baseline" model of Stecker et al. (2006) can be ruled out with high confidence.

The Fermi Large Area Telescope (LAT) data have confirmed the pulsed emission from all six high-confidence gamma-ray pulsars previously known from the EGRET observations. We report results obtained from the analysis of 13 months of LAT data for three of these pulsars (PSR J1057-5226, PSR J1709-4429, and PSR J1952+3252) each of which had some unique feature among the EGRET pulsars. The excellent sensitivity of LAT allows more detailed analysis of the evolution of the pulse profile with energy and also of the variation of the spectral shape with phase. We measure the cutoff energy of the pulsed emission from these pulsars for the first time and provide a more complete picture of the emission mechanism. The results confirm some, but not all, of the features seen in the EGRET data.

We report on the \textit{Fermi}-LAT observations of the Geminga pulsar, the second brightest non-variable GeV source in the $\gamma$-ray sky and the first example of a radio-quiet $\gamma$-ray pulsar. The observations cover one year, from the launch of the $Fermi$ satellite through 2009 June 15. A data sample of over 60,000 photons enabled us to build a timing solution based solely on $\gamma$ rays. Timing analysis shows two prominent peaks, separated by $\Delta \phi$ = 0.497 $\pm$ 0.004 in phase, which narrow with increasing energy. Pulsed $\gamma$ rays are observed beyond 18 GeV, precluding emission below 2.7 stellar radii because of magnetic absorption. The phase-averaged spectrum was fitted with a power law with exponential cut-off of spectral index $\Gamma$ = (1.30 $\pm$ 0.01 $\pm$ 0.04), cut-off energy $E_{0}$ = (2.46 $\pm$ 0.04 $\pm$ 0.17) GeV and an integral photon flux above 0.1 GeV of (4.14 $\pm$ 0.02 $\pm$ 0.32) $\times$ 10$^{-6}$ cm$^{-2}$ s$^{-1}$. The first uncertainties are statistical and the second are systematic. The phase-resolved spectroscopy shows a clear evolution of the spectral parameters, with the spectral index reaching a minimum value just before the leading peak and the cut-off energy having maxima around the peaks. Phase-resolved spectroscopy reveals that pulsar emission is present at all rotational phases. The spectral shape, broad pulse profile, and maximum photon energy favor the outer magnetospheric emission scenarios.

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.

We present detailed analysis of the two gamma-ray sources,1FGL J1801.3-2322c and 1FGL J1800.5-2359c,that have been found toward the supernova remnant(SNR) W28 with the Large Area Telescope(LAT) on board the Fermi Gamma-ray Space Telescope.1FGL J1801.3-2322c is found to be an extended source within the boundary of SNR W28,and to extensively overlap with the TeV gamma-ray source HESS J1801-233,which is associated with a dense molecular cloud interacting with the supernova remnant.The gamma-ray spectrum measured with LAT from 0.2--100 GeV can be described by a broken power-law function with a break of ~1GeV,and photon indices of 2.09$\pm$0.08(stat)$\pm$0.28(sys) below the break and 2.74$\pm$0.06(stat)$\pm$0.09(sys) above the break.Given the clear association between HESS J1801-233 and the shocked molecular cloud and a smoothly connected spectrum in the GeV--TeV band,we consider the origin of the gamma-ray emission in both GeV and TeV ranges to be the interaction between particles accelerated in the SNR and the molecular cloud.The decay of neutral pions produced in interactions between accelerated hadrons and dense molecular gas provide a reasonable explanation for the broadband gamma-ray spectrum. 1FGL J1800.5-2359c, located outside the southern boundary of SNR W28, cannot be resolved.An upper limit on the size of the gamma-ray emission was estimated to be ~16$'$ using events above ~2GeV under the assumption of a circular shape with uniform surface brightness. It appears to coincide with the TeV source HESS J1800-240B,which is considered to be associated with a dense molecular cloud that contains the ultra compact HII region W28A2(G5.89-0.39).We found no significant gamma-ray emission in the LAT energy band at the positions of TeV sources HESS J1800-230A and HESS J1800-230C.The LAT data for HESS J1800-230A combined with the TeV data points indicate a spectral break between 10GeV and 100GeV.

We present the results from the Suzaku X-ray observations of five flat-spectrum radio quasars (FSRQs), namely PKS0208-512, Q0827+243, PKS1127-145, PKS1510-089 and 3C 454.3. All these sources were additionally monitored simultaneously or quasi-simultaneously by the Fermi satellite in gamma-rays and the Swift UVOT in the UV and optical bands, respectively. We constructed their broad-band spectra covering the frequency range from 10^14 Hz up to 10^25 Hz, and those reveal the nature of high-energy emission of luminous blazars in their low-activity states. The analyzed X-ray spectra are well fitted by a power-law model with photoelectric absorption. In the case of PKS0208-512, PKS1127-145, and 3C 454.3, the X-ray continuum showed indication of hard-ening at low-energies. Moreover, when compared with the previous X-ray observations, we see a significantly increasing contribution of low-energy photons to the total X-ray fluxes when the sources are getting fainter. The same behavior can be noted in the Suzaku data alone. A likely explanation involves a variable, flat-spectrum component produced via inverse-Compton (IC) emission, plus an additional, possibly steady soft X-ray component prominent when the source gets fainter. This soft X-ray excess is represented either by a steep powerlaw (photon indices Gamma ~ 3 - 5) or a blackbody-type emission with temperatures kT ~ 0.1-0.2 keV. We model the broad-band spectra spectra of the five observed FSRQs using synchrotron self-Compton (SSC) and/or external-Compton radiation (ECR) models. Our modeling suggests that the difference between the low and high-activity states in luminous blazars is due to the different total kinetic power of the jet, most likely related to varying bulk Lorentz factor of the outflow within the blazar emission zone.

This paper presents light curves and the first systematic characterization of variability of the 106 objects in the Fermi Large Area Telescope (LAT) Bright AGN Sample (LBAS). Weekly light curves obtained during the first 11 months of survey (August 04, 2008 - July 04, 2009), are tested for variability, and their properties are quantified through autocorrelation and structure function analysis. For the brightest sources power density spectra (PDS) and fit of the temporal structure of major flares is performed. More than 50% of the sources are variable, where high states do not exceed 1/4 of the total observation range. Variation amplitudes are larger for FSRQs and low/intermediate synchrotron peaked (LSP/ISP) BL Lac objects. Autocorrelation time scales vary from 4 to a dozen of weeks. Variable sources of the sample have 1/(f^{a}) PDS and show two modes: (1) rather constant baseline with sporadic flaring activity characterized by flatter PDS slopes resembling flickering and red-noise with occasional intermittence, and (2) - measured for a few blazars showing strong activity - complex and structured temporal profiles characterized by longer-term memory and steeper PDS slopes typical of a random-walk underlying mechanism. The average PDS slope of the brightest 22 FSRQs and the 6 brightest BL Lacs is 1.5 and 1.7 respectively. The study of temporal profiles of well resolved flares shows that they generally have symmetric profiles and that their total duration vary between 10 and 100 days.

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.

We present the first catalog of active galactic nuclei (AGN) detected by the LAT, corresponding to 11 months of data collected in scientific operation mode. The First LAT AGN Catalog (1LAC) includes 671 gamma-ray sources located at high Galactic latitudes (|b| > 10 deg) that are detected with a test statistic (TS) greater than 25 and associated statistically with AGNs. Some LAT sources are associated with multiple AGNs, and consequently, the catalog includes 709 AGNs, comprising 300 BL Lacertae objects (BL Lacs), 296 flat-spectrum radio quasars (FSRQs), 41 AGNs of other types, and 72 AGNs of unknown type. We also classify the blazars based on their spectral energy distributions (SEDs) as archival radio, optical, and X-ray data permit. In addition to the format 1LAC sample, we provide AGN associations for 51 low-latitude LAT sources and AGN "affiliations" (unquantified counterpart candidates) for 104 high-latitude LAT sources without AGN associations. The overlap of the 1LAC with existing gamma-ray AGN catalogs (LBAS, EGRET, AGILE, Swift, INTEGRAL, TeVCat) is briefly discussed. Various properties--such as gamma-ray fluxes and photon power law spectral indices, redshifts, gamma-ray luminosities, variability, and archival radio luminosities--and their correlations are presented and discussed for the different blazar classes. We compare the 1LAC results with predictions regarding the gamma-ray AGN populations, and we comment on the power of the sample to address the question of the blazar sequence.

We report on analysis of timing and spectroscopy of the Vela pulsar using eleven months of observations with the Large Area Telescope on the Fermi Gamma-Ray Space Telescope. The intrinsic brightness of Vela at GeV energies combined with the angular resolution and sensitivity of the LAT allow us to make the most detailed study to date of the energy-dependent light curves and phase-resolved spectra, using a LAT-derived timing model. The light curve consists of two peaks (P1 and P2) connected by bridge emission containing a third peak (P3). We have confirmed the strong decrease of the P1/P2 ratio with increasing energy seen with EGRET and previous Fermi LAT data, and observe that P1 disappears above 20 GeV. The increase with energy of the mean phase of the P3 component can be followed with much greater detail, showing that P3 and P2 are present up to the highest energies of pulsation. We find significant pulsed emission at phases outside the main profile, indicating that magnetospheric emission exists over 80% of the pulsar period. With increased high-energy counts the phase-averaged spectrum is seen to depart from a power- law with simple exponential cutoff, and is better fit with a more gradual cutoff. The spectra in fixed-count phase bins are well fit with power-laws with exponential cutoffs, revealing a strong and complex phase dependence of the cutoff energy, especially in the peaks. By combining these results with predictions of the outer magnetosphere models that map emission characteristics to phase, it will be possible to probe the particle acceleration and the structure of the pulsar magnetosphere with unprecedented detail.

We report observation of the supernova remnant IC443 (G189.1+3.0) with the Fermi Gamma-ray Space Telescope Large Area Telescope (LAT) in the energy band between 200MeV and 50GeV. IC443 is a shell-type supernova remnant with mixed morphology located off the outer Galactic plane where high-energy emission has been detected in the X-ray, GeV and TeV gamma-ray bands. Past observations suggest IC443 has been interacting with surrounding interstellar matter. Proximity between dense shocked molecular clouds and GeV-TeV gamma-ray emission regions detected by EGRET, MAGIC and VERITAS suggests an interpretation that cosmic-ray (CR) particles are accelerated by the SNR. With the high gamma-ray statistics and broad energy coverage provided by the LAT, we accurately characterize the gamma-ray emission produced by the CRs accelerated at IC443. The emission region is extended in the energy band with theta_68 = 0.27 deg +/- 0.01 deg (stat) +/- 0.03 deg (sys) for an assumed 2-dimensional Gaussian profile and overlaps almost completely with the extended source region of VERITAS. Its centroid is displaced significantly from the known pulsar wind nebula (PWN) which suggests the PWN is not the major contributor in the present energy band. The observed spectrum changes its power-law slope continuously and continues smoothly to the MAGIC and VERITAS data points. The combined gamma-ray spectrum (200MeV <E< 2TeV) is reproduced well by decays of neutral pions produced by a broken power-law proton spectrum with a break around 70GeV.

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)

Context: The Large Magellanic Cloud (LMC) is to date the only normal external galaxy that has been detected in high-energy gamma rays. High-energy gamma rays trace particle acceleration processes and gamma-ray observations allow the nature and sites of acceleration to be studied. Aims: We characterise the distribution and sources of cosmic rays in the LMC from analysis of gamma-ray observations. Methods: We analyse 11 months of continuous sky-survey observations obtained with the Large Area Telescope aboard the Fermi Gamma-Ray Space Telescope and compare it to tracers of the interstellar medium and models of the gamma-ray sources in the LMC. Results: The LMC is detected at 33 sigma significance. The integrated >100 MeV photon flux of the LMC amounts to (2.6 +/- 0.2) * 10^-7 ph/cm2/s which corresponds to an energy flux of (1.6 +/- 0.1) * 10^-10 erg/cm2/s, with additional systematic uncertainties of ~16%. The analysis reveals the massive star forming region 30 Doradus as a bright source of gamma-ray emission in the LMC in addition to fainter emission regions found in the northern part of the galaxy. The gamma-ray emission from the LMC shows very little correlation with gas density and is rather correlated to tracers of massive star forming regions. The close confinement of gamma-ray emission to star forming regions suggests a relatively short GeV cosmic-ray proton diffusion length. Conclusions: The close correlation between cosmic-ray density and massive star tracers supports the idea that cosmic rays are accelerated in massive star forming regions as a result of the large amounts of kinetic energy that are input by the stellar winds and supernova explosions of massive stars into the interstellar medium.

We report on the first detection of GeV high-energy gamma-ray emission from a young supernova remnant with the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope. These observations reveal a source with no discernible spatial extension detected at a significance level of 12.2$\sigma$ above 500 MeV at a location that is consistent with the position of the remnant of the supernova explosion that occurred around 1680 in the Cassiopeia constellation - Cassiopeia A. The gamma-ray flux and spectral shape of the source are consistent with a scenario in which the gamma-ray emission originates from relativistic particles accelerated in the shell of this remnant. The total content of cosmic rays (electrons and protons) accelerated in Cas A can be estimated as $W_{\mathrm{CR}} \approx (1-4) \times 10^{49}$ erg thanks to the well-known density in the remnant assuming that the observed gamma-ray originates in the SNR shell(s). The magnetic field in the radio-emitting plasma can be robustly constrained as B $\gt 0.1$ mG, providing new evidence of the magnetic field amplification at the forward shock and the strong field in the shocked ejecta.

We present multiwavelength studies of the 106.6 ms gamma-ray pulsar PSR J1907+06 near the TeV source MGRO J1908+06. Timing observations with Fermi result in a precise position determination for the pulsar of R.A. = 19h07m547(2), decl. = +06:02:16(2) placing the pulsar firmly within the TeV source extent, suggesting the TeV source is the pulsar wind nebula of PSR J1907+0602. Pulsed gamma-ray emission is clearly visible at energies from 100 MeV to above 10 GeV. The phase-averaged power-law index in the energy range E > 0.1 GeV is = 1.76 \pm 0.05 with an exponential cutoff energy E_{c} = 3.6 \pm 0.5 GeV. We present the energy-dependent gamma-ray pulsed light curve as well as limits on off-pulse emission associated with the TeV source. We also report the detection of very faint (flux density of ~3.4 microJy) radio pulsations with the Arecibo telescope at 1.5 GHz having a dispersion measure DM = 82.1 \pm 1.1 cm^{-3}pc. This indicates a distance of 3.2 \pm 0.6 kpc and a pseudo-luminosity of L_{1400} ~ 0.035 mJy kpc^2. A Chandra ACIS observation revealed an absorbed, possibly extended, compact <(4 arcsec) X-ray source with significant non-thermal emission at R.A. = 19h07m54.76, decl. = +06:02:14.6 with a flux of 2.3^{+0.6}_{-1.4} X 10^{-14} erg cm^{-2} s^{-1}. From archival ASCA observations, we place upper limits on any arcminute scale 2--10 keV X-ray emission of ~ 1 X 10^{-13} erg cm^{-2} s^{-1}. The implied distance to the pulsar is compatible with that of the supernova remnant G40.5-0.5, located on the far side of the TeV nebula from PSR J1907+0602, and the S74 molecular cloud on the nearer side which we discuss as potential birth sites.

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

We report the first detection of very-high-energy (VHE) gamma-ray emission above 140 GeV from PKS 1424+240, a BL Lac object with an unknown redshift. The photon spectrum above 140 GeV measured by VERITAS is well described by a power law with a photon index of 3.8 +- 0.5_stat +- 0.3_syst and a flux normalization at 200 GeV of (5.1 +- 0.9_stat +- 0.5_syst) x 10^{-11} TeV^-1 cm^-2 s^-1, where stat and syst denote the statistical and systematical uncertainty, respectively. The VHE flux is steady over the observation period between MJD 54881 and 55003 (2009 February 19 to June 21). Flux variability is also not observed in contemporaneous high energy observations with the Fermi Large Area Telescope (LAT). Contemporaneous X-ray and optical data were also obtained from the Swift XRT and MDM observatory, respectively. The broadband spectral energy distribution (SED) is well described by a one-zone synchrotron self-Compton (SSC) model favoring a redshift of less than 0.1. Using the photon index measured with Fermi in combination with recent extragalactic background light (EBL) absorption models it can be concluded from the VERITAS data that the redshift of PKS 1424+240 is less than 0.66.

We report the detection of high-energy gamma-ray emission from two starburst galaxies using data obtained with the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. Steady point-like emission above 200 MeV has been detected at significance levels of 6.8 sigma and 4.8 sigma respectively, from sources positionally coincident with locations of the starburst galaxies M82 and NGC 253. The total fluxes of the sources are consistent with gamma-ray emission originating from the interaction of cosmic rays with local interstellar gas and radiation fields and constitute evidence for a link between massive star formation and gamma-ray emission in star-forming galaxies.

We report on gamma-ray observations of the Crab Pulsar and Nebula using 8 months of survey data with the Fermi Large Area Telescope (LAT). The high quality light curve obtained using the ephemeris provided by the Nancay and Jodrell Bank radio telescopes shows two main peaks stable in phase with energy. The first gamma-ray peak leads the radio main pulse by (281 \pm 12 \pm 21) mus, giving new constraints on the production site of non-thermal emission in pulsar magnetospheres. The improved sensitivity and the unprecedented statistics afforded by the LAT enable precise measurement of the Crab Pulsar spectral parameters: cut-off energy at E_c = (5.8 \pm 0.5 \pm 1.2) GeV, spectral index of Gamma = (1.97 \pm 0.02 \pm 0.06) and integral photon flux above 100 MeV of (2.09 \pm 0.03 \pm 0.18) x 10^{-6} cm^{-2} s^{-1}. The first errors represent the statistical error on the fit parameters, while the second ones are the systematic uncertainties. Pulsed gamma-ray photons are observed up to ~ 20 GeV which precludes emission near the stellar surface, below altitudes of around 4 to 5 stellar radii in phase intervals encompassing the two main peaks. The spectrum of the nebula in the energy range 100 MeV - 300 GeV is well described by the sum of two power-laws of indices Gamma_{sync} = (3.99 \pm 0.12 \pm 0.08) and Gamma_{IC} = (1.64 \pm 0.05 \pm 0.07), corresponding to the falling edge of the synchrotron and the rising edge of the inverse Compton components, respectively. This latter, which links up naturally with the spectral data points of Cherenkov experiments, is well reproduced via inverse Compton scattering from standard Magnetohydrodynamics (MHD) nebula models, and does not require any additional radiation mechanism.

We report on observations of TeV-selected AGN made during the first 5.5 months of observations with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope (Fermi). In total, 96 AGN were selected for study, each being either (i) a source detected at TeV energies (28 sources) or (ii) an object that has been studied with TeV instruments and for which an upper-limit has been reported (68 objects). The Fermi observations show clear detections of 38 of these TeV-selected objects, of which 21 are joint GeV-TeV sources and 29 were not in the third EGRET catalog. For each of the 38 Fermi-detected sources, spectra and light curves are presented. Most can be described with a power law of spectral index harder than 2.0, with a spectral break generally required to accommodate the TeV measurements. Based on an extrapolation of the Fermi spectrum, we identify sources, not previously detected at TeV energies, which are promising targets for TeV instruments. Evidence for systematic evolution of the $\gamma$-ray spectrum with redshift is presented and discussed in the context of interaction with the EBL.

We report the Fermi-LAT discovery of high-energy (MeV/GeV) gamma-ray emission positionally consistent with the center of the radio galaxy M87, at a source significance of over 10 sigma in ten-months of all-sky survey data. Following the detections of Cen A and Per A, this makes M87 the third radio galaxy seen with the LAT. The faint point-like gamma-ray source has a >100 MeV flux of 2.45 (+/- 0.63) x 10^-8 ph cm^-2 s^-1 (photon index = 2.26 +/- 0.13) with no significant variability detected within the LAT observation. This flux is comparable with the previous EGRET upper limit (< 2.18 x 10^-8 ph cm^-2 s^-1, 2 sigma), thus there is no evidence for a significant MeV/GeV flare on decade timescales. Contemporaneous Chandra and VLBA data indicate low activity in the unresolved X-ray and radio core relative to previous observations, suggesting M87 is in a quiescent overall level over the first year of Fermi-LAT observations. The LAT gamma-ray spectrum is modeled as synchrotron self-Compton (SSC) emission from the electron population producing the radio-to-X-ray emission in the core. The resultant SSC spectrum extrapolates smoothly from the LAT band to the historical-minimum TeV emission. Alternative models for the core and possible contributions from the kiloparsec-scale jet in M87 are considered, and can not be excluded.

We present the result of a search of the Milagro sky map for spatial correlations with sources from a subset of the recent Fermi Bright Source List (BSL). The BSL consists of the 205 most significant sources detected above 100 MeV by the Fermi Large Area Telescope. We select sources based on their categorization in the BSL, taking all confirmed or possible Galactic sources in the field of view of Milagro. Of the 34 Fermi sources selected, 14 are observed by Milagro at a significance of 3 standard deviations or more. We conduct this search with a new analysis which employs newly-optimized gamma-hadron separation and utilizes the full 8-year Milagro dataset. Milagro is sensitive to gamma rays with energy from 1 to 100 TeV with a peak sensitivity from 10-50 TeV depending on the source spectrum and declination. These results extend the observation of these sources far above the Fermi energy band. With the new analysis and additional data, multi-TeV emission is definitively observed associated with the Fermi pulsar, J2229.0+6114, in the Boomerang Pulsar Wind Nebula (PWN). Furthermore, an extended region of multi-TeV emission is associated with the Fermi pulsar, J0634.0+1745, the Geminga pulsar.

Observations by the Large Area Telescope (LAT) on the \textit{Fermi} mission of diffuse $\gamma$-rays in a mid-latitude region in the third quadrant (Galactic longitude $l$ from $200\arcdeg$ to $260\arcdeg$ and latitude $| b |$ from $22\arcdeg$ to $60\arcdeg$) are reported. The region contains no known large molecular cloud and most of the atomic hydrogen is within 1 kpc of the solar system. The contributions of $\gamma$-ray point sources and inverse Compton scattering are estimated and subtracted. The residual $\gamma$-ray intensity exhibits a linear correlation with the atomic gas column density in energy from 100 MeV to 10 GeV. The measured integrated $\gamma$-ray emissivity is $(1.63 \pm 0.05) \times 10^{-26} {\rm photons s^{-1} sr^{-1} H\mathchar`-atom^{-1}}$ and $(0.66 \pm 0.02) \times 10^{-26} {\rm photons s^{-1} sr^{-1} H\mathchar`-atom^{-1}}$ above 100 MeV and above 300 MeV, respectively, with additional systematic error of $\sim 10%$. The differential emissivity in 100 MeV--10 GeV agrees with calculations based on cosmic ray spectra consistent with those directly measured, at the 10% level. The results obtained indicate that cosmic ray nuclei spectra within 1 kpc from the solar system in regions studied are close to the local interstellar spectra inferred from direct measurements at the Earth within $\sim 10%$.

Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio pulsations from two of them. PSR J1741-2054, with period P=413ms, was detected in archival Parkes telescope data and subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to interstellar scintillation and it has a very small dispersion measure of DM=4.7pc/cc, implying a distance of ~0.4kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic emission, its gamma-ray luminosity above 0.1GeV corresponds to 25% of the spin-down luminosity of dE/dt=9.4e33erg/s. The gamma-ray profile occupies 1/3 of pulse phase and has three closely-spaced peaks with the first peak lagging the radio pulse by delta=0.29P. We have also identified a soft Swift source that is the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second source, PSR J2032+4127, was detected at the GBT. It has P=143ms, and its DM=115pc/cc suggests a distance of ~3.6kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by delta=0.15P the first of two narrow gamma-ray peaks that are separated by Delta=0.50P. Faint, diffuse X-ray emission in a Chandra image is possibly its pulsar wind nebula. PSR J2032+4127 likely accounts for the EGRET source 3EG J2033+4118, while its pulsar wind is responsible for the formerly unidentified HEGRA source TeV J2032+4130.

This Letter presents the first results from the observations of LSI +61 303 using Large Area Telescope data from the Fermi Gamma-Ray Space Telescope between 2008 August and 2009 March. Our results indicate variability that is consistent with the binary period, with the emission being modulated at 26.6 +/- 0.5 days. This constitutes the first detection of orbital periodicity in high-energy gamma rays (20 MeV-100 GeV, HE). The light curve is characterized by a broad peak after periastron, as well as a smaller peak just before apastron. The spectrum is best represented by a power law with an exponential cutoff, yielding an overall flux above 100 MeV of 0.82 +/- 0.03(stat) +/- 0.07(syst) 10^{-6} ph cm^{-2} s^{-1}, with a cutoff at 6.3 +/- 1.1(stat) +/- 0.4(syst) GeV and photon index Gamma = 2.21 +/- 0.04(stat) +/- 0.06(syst). There is no significant spectral change with orbital phase. The phase of maximum emission, close to periastron, hints at inverse Compton scattering as the main radiation mechanism. However, previous very high-energy gamma ray (>100 GeV, VHE) observations by MAGIC and VERITAS show peak emission close to apastron. This and the energy cutoff seen with Fermi suggest the link between HE and VHE gamma rays is nontrivial.

We report the discovery by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope of high-energy gamma-ray emission from the peculiar quasar PMN J0948+0022 (z=0.5846). The optical spectrum of this object exhibits rather narrow Hbeta (FWHM(Hbeta) ~ 1500 km s^-1), weak forbidden lines and is therefore classified as a narrow-line type I quasar. This class of objects is thought to have relatively small black hole mass and to accrete at high Eddington ratio. The radio loudness and variability of the compact radio core indicates the presence of a relativistic jet. Quasi simultaneous radio-optical-X-ray and gamma-ray observations are presented. Both radio and gamma-ray emission (observed over 5-months) are strongly variable. The simultaneous optical and X-ray data from Swift show a blue continuum attributed to the accretion disk and a hard X-ray spectrum attributed to the jet. The resulting broad band spectral energy distribution (SED) and, in particular, the gamma-ray spectrum measured by Fermi are similar to those of more powerful FSRQ. A comparison of the radio and gamma-ray characteristics of PMN J0948+0022 with the other blazars detected by LAT shows that this source has a relatively low radio and gamma-ray power, with respect to other FSRQ. The physical parameters obtained from modelling the SED also fall at the low power end of the FSRQ parameter region discussed in Celotti & Ghisellini (2008). We suggest that the similarity of the SED of PMN J0948+0022 to that of more massive and more powerful quasars can be understood in a scenario in which the SED properties depend on the Eddington ratio rather than on the absolute power.

We report the detection of pulsed gamma-rays from the young, spin-powered radio pulsar PSR J2021+3651 using data acquired with the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope (formerly GLAST). The light curve consists of two narrow peaks of similar amplitude separated by 0.468 +/- 0.002 in phase. The first peak lags the maximum of the 2 GHz radio pulse by 0.162 +/- 0.004 +/- 0.01 in phase. The integral gamma-ray photon flux above 100 MeV is (56 +/- 3 +/- 11) x 10^{-8} /cm2/s. The photon spectrum is well-described by an exponentially cut-off power law of the form dF/dE = kE^{-\Gamma} e^(-E/E_c) where the energy E is expressed in GeV. The photon index is \Gamma = 1.5 +/- 0.1 +/- 0.1 and the exponential cut-off is E_c = 2.4 +/- 0.3 +/- 0.5 GeV. The first uncertainty is statistical and the second is systematic. The integral photon flux of the bridge is approximately 10% of the pulsed emission, and the upper limit on off-pulse gamma-ray emission from a putative pulsar wind nebula is <10% of the pulsed emission at the 95% confidence level. Radio polarization measurements yield a rotation measure of RM = 524 +/- 4 rad/m^2 but a poorly constrained magnetic geometry. Re-analysis of Chandra data enhanced the significance of the weak X-ray pulsations, and the first peak is roughly phase-aligned with the first gamma-ray peak. We discuss the emission region and beaming geometry based on the shape and spectrum of the gamma-ray light curve combined with radio and X-ray measurements, and the implications for the pulsar distance. Gamma-ray emission from the polar cap region seems unlikely for this pulsar.

Following its launch in June 2008, the Fermi Gamma-ray Space Telescope (Fermi) began a sky survey in August. The Large Area Telescope (LAT) on Fermi in 3 months produced a deeper and better-resolved map of the gamma-ray sky than any previous space mission. We present here initial results for energies above 100 MeV for the 205 most significant (statistical significance greater than ~10-sigma) gamma-ray sources in these data. These are the best-characterized and best-localized point-like (i.e., spatially unresolved) gamma-ray sources in the early-mission data.

This is the first report of Fermi Gamma-ray Space Telescope observations of the quasar 3C 454.3, which has been undergoing pronounced long-term outbursts since 2000. The data from the Large Area Telescope (LAT), covering 2008 July 7 - October 6, indicate strong, highly variable gamma-ray emission with an average flux of ~3 x 10^{-6} photons cm^{-2} s^{-1}, for energies above 100 MeV. The gamma-ray flux is variable, with strong, distinct, symmetrically-shaped flares for which the flux increases by a factor of several on a time scale of about three days. This variability indicates a compact emission region, and the requirement that the source is optically thin to pair-production implies relativistic beaming with Doppler factor delta > 8, consistent with the values inferred from VLBI observations of superluminal expansion (delta ~ 25). The observed gamma-ray spectrum is not consistent with a simple power-law, but instead steepens strongly above ~2 GeV, and is well described by a broken power-law with photon indices of ~2.3 and ~3.5 below and above the break, respectively. This is the first direct observation of a break in the spectrum of a high luminosity blazar above 100 MeV, and it is likely direct evidence for an intrinsic break in the energy distribution of the radiating particles. Alternatively, the spectral softening above 2 GeV could be due to gamma-ray absorption via photon-photon pair production on the soft X-ray photon field of the host AGN, but such an interpretation would require the dissipation region to be located very close (less than 100 gravitational radii) to the black hole, which would be inconsistent with the X-ray spectrum of the source.

Results are presented of a harmonic analysis of the large scale cosmic-ray anisotropy as observed by the Milagro observatory. We show a two-dimensional display of the sidereal anisotropy pro jections in right ascension generated by the fitting of three harmonics to 18 separate declination bands. The Milagro observatory is a water Cherenkov detector located in the Jemez mountains near Los Alamos, New Mexico. With a high duty cycle and large field-of-view, Milagro is an excellent instrument for measuring this anisotropy with high sensitivity at TeV energies. The analysis is conducted using a seven year data sample consisting of more than 95 billion events, the largest such data set in existence. We observe an anisotropy with a magnitude around 0.1% for cosmic rays with a median energy of 6 TeV. The dominant feature is a deficit region of depth (2.49 +/- 0.02 stat. +/- 0.09 sys.)x10^(-3) in the direction of the Galactic North Pole centered at 189 degrees right ascension. We observe a steady increase in the magnitude of the signal over seven years.

We report the discovery of high-energy (E>100 MeV) gamma-ray emission from NGC 1275, a giant elliptical galaxy lying at the center of the Perseus cluster of galaxies, based on observations made with the Large Area Telescope (LAT) of the Fermi Gamma ray Space Telescope. The positional center of the gamma-ray source is only ~3' away from the NGC 1275 nucleus, well within the 95% LAT error circle of ~5'.The spatial distribution of gamma-ray photons is consistent with a point source. The average flux and power-law photon index measured with the LAT from 2008 August 4 to 2008 December 5 are F_gamma = (2.10+-0.23)x 10^{-7} ph (>100 MeV) cm^{-2} s^{-1} and Gamma = 2.17+-0.05, respectively. The measurements are statistically consistent with constant flux during the four-month LAT observing period.Previous EGRET observations gave an upper limit of F_gamma < 3.72x 10 ^{-8} ph (>100 MeV) cm^{-2} s^{-1} to the gamma-ray flux from NGC 1275. This indicates that the source is variable on timescales of years to decades, and therefore restricts the fraction of emission that can be produced in extended regions of the galaxy cluster. Contemporaneous and historical radio observations are also reported. The broadband spectrum of NGC 1275 is modeled with a simple one-zone synchrotron/synchrotron self-Compton model and a model with a decelerating jet flow.

We report the discovery by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope of high-energy gamma-ray (GeV) emission from the flat-spectrum radio quasar PKS 1454-354 (z=1.424). On 4 September 2008 the source rose to a peak flux of (3.5 +/- 0.7)x 10^-6 ph cm^-2 s^-1 (E > 100 MeV) on a time scale of hours and then slowly dropped over the following two days. No significant spectral changes occurred during the flare. Fermi/LAT observations also showed that PKS 1454-354 is the most probable counterpart of the unidentified EGRET source 3EG J1500-3509. Multiwavelength measurements performed during the following days (7 September with Swift; 6-7 September with the ground-based optical telescope ATOM; 13 September with the Australia Telescope Compact Array) resulted in radio, optical, UV and X-ray fluxes greater than archival data, confirming the activity of PKS 1454-354.

The first three months of sky-survey operation with the Fermi Gamma Ray Space Telescope (Fermi) Large Area Telescope (LAT) reveals 132 bright sources at |b|>10 deg with test statistic greater than 100 (corresponding to about 10 sigma). Two methods, based on the CGRaBS, CRATES and BZCat catalogs, indicate high-confidence associations of 106 of these sources with known AGNs. This sample is referred to as the LAT Bright AGN Sample (LBAS). It contains two radio galaxies, namely Centaurus A and NGC 1275, and 104 blazars consisting of 57 flat spectrum radio quasars (FSRQs), 42 BL Lac objects, and 5 blazars with uncertain classification. Four new blazars were discovered on the basis of the LAT detections. Remarkably, the LBAS includes 10 high-energy peaked BL Lacs (HBLs), sources which were so far hard to detect in the GeV range. Another 10 lower-confidence associations are found. Only thirty three of the sources, plus two at |b|>10 deg, were previously detected with EGRET, probably due to the variable nature of these sources. The analysis of the gamma-ray properties of the LBAS sources reveals that the average GeV spectra of BL Lac objects are significantly harder than the spectra of FSRQs. No significant correlation between radio and peak gamma-ray fluxes is observed. Blazar log N - log S and luminosity functions are constructed to investigate the evolution of the different blazar classes, with positive evolution indicated for FSRQs but none for BLLacs. The contribution of LAT-blazars to the total extragalactic gamma-ray intensity is estimated.

The Vela pulsar is the brightest persistent source in the GeV sky and thus is the traditional first target for new gamma-ray observatories. We report here on initial Fermi Large Area Telescope observations during verification phase pointed exposure and early sky survey scanning. We have used the Vela signal to verify Fermi timing and angular resolution. The high quality pulse profile, with some 32,400 pulsed photons at E>0.03 GeV, shows new features, including pulse structure as fine as 0.3ms and a distinct third peak, which shifts in phase with energy. We examine the high energy behavior of the pulsed emission; initial spectra suggest a phase-averaged power law index of Gamma=1.51{+0.05/-0.04} with an exponential cut-off at E_c=2.9+/-0.1 GeV. Spectral fits with generalized cut-offs of the form e^{-(E/E_c)^b} require b<1, which is inconsistent with magnetic pair attenuation, and thus favor outer magnetosphere emission models. Finally, we report on upper limits to any unpulsed component, as might be associated with a surrounding synchrotron wind nebula (PWN).

An analysis of 7 years of Milagro data performed on a 10-degree angular scale has found two localized regions of excess of unknown origin with greater than 12 sigma significance. Both regions are inconsistent with gamma-ray emission with high confidence. One of the regions has a different energy spectrum than the isotropic cosmic-ray flux at a level of 4.6 sigma, and it is consistent with hard spectrum protons with an exponential cutoff, with the most significant excess at ~10 TeV. Potential causes of these excesses are explored, but no compelling explanations are found.

This is a report on the findings of the SNR/cosmic-ray working group for the white paper on the status and future of ground-based gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

This is a report on the findings of the Galactic compact objects working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper is an APS commissioned document, and the full version has also been released and can be found on astro-ph. This section of the white paper discusses the potential of future ground-based gamma-ray experiments to advance our understanding of the physics of Galactic compact objects including pulsars, pulsar wind nebulae, and X-ray binaries.

Diffuse $\gamma$-ray emission produced by the interaction of cosmic-ray particles with matter and radiation in the Galaxy can be used to probe the distribution of cosmic rays and their sources in different regions of the Galaxy. With its large field of view and long observation time, the Milagro Gamma Ray Observatory is an ideal instrument for surveying large regions of the Northern Hemisphere sky and for detecting diffuse $\gamma$-ray emission at very high energies. Here, the spatial distribution and the flux of the diffuse $\gamma$-ray emission in the TeV energy range with a median energy of 15 TeV for Galactic longitudes between 30$^\circ$ and 110$^\circ$ and between 136$^\circ$ and 216$^\circ$ and for Galactic latitudes between -10$^\circ$ and 10$^\circ$ are determined. The measured fluxes are consistent with predictions of the GALPROP model everywhere except for the Cygnus region ($l\in[65^\circ,85^\circ]$). For the Cygnus region, the flux is twice the predicted value. This excess can be explained by the presence of active cosmic ray sources accelerating hadrons which interact with the local dense interstellar medium and produce gamma rays through pion decay.

Recent rapid localizations of short, hard gamma-ray bursts (GRBs) by the Swift and HETE satellites have led to the observation of the first afterglows and the measurement of the first redshifts from this type of burst. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Seventeen short duration (< 5 s) GRBs detected by satellites occurred within the field of view of the Milagro gamma-ray observatory between 2000 January and 2006 December. We have searched the Milagro data for >100 GeV counterparts to these GRBs and find no significant emission correlated with these bursts. Due to the absorption of high-energy gamma rays by the extragalactic background light (EBL), detections are only expected for redshifts less than ~0.5. While most long duration GRBs occur at redshifts higher than 0.5, the opposite is thought to be true of short GRBs. Lack of a detected VHE signal thus allows setting meaningful fluence limits. One GRB in the sample (050509b) has a likely association with a galaxy at a redshift of 0.225, while another (051103) has been tentatively linked to the nearby galaxy M81. Fluence limits are corrected for EBL absorption, either using the known measured redshift, or computing the corresponding absorption for a redshift of 0.1 and 0.5, as well as for the case of z=0.

A survey of Galactic gamma-ray sources at a median energy of ~20 TeV has been performed using the Milagro Gamma Ray Observatory. Eight candidate sources of TeV emission are detected with pre-trials significance $>4.5\sigma$ in the region of Galactic longitude $l\in[30^\circ,220^\circ]$ and latitude $b\in[-10^\circ,10^\circ]$. Four of these sources, including the Crab nebula and the recently published MGRO J2019+37, are observed with significances $>4\sigma$ after accounting for the trials involved in searching the 3800 square degree region. All four of these sources are also coincident with EGRET sources. Two of the lower significance sources are coincident with EGRET sources and one of these sources is Geminga. The other two candidates are in the Cygnus region of the Galaxy. Several of the sources appear to be spatially extended. The fluxes of the sources at 20 TeV range from ~25% of the Crab flux to nearly as bright as the Crab.

The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.